//*  THIS DATASET IS NAMED 'STORET.HELP.FLOW.LIBRARY(SAS01)'            00010004
//*                                                                     00020104
/*CNTL CWT.TRY,EXC                                                      00020204
//FLOWSTEP  EXEC  FLOW                                                  00020306
//DISTR.INPUT DD *                                                      00020404
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # 00020800
#  THE FOLLOWING RETRIEVAL CONTROL CARDS REQUEST MEAN DAILY FLOW DATA # 00020900
#  FOR WATER YEARS 1964 THROUGH 1966 AT STATION 09268900.  DATA ARE   # 00021000
#  NOT DISTINGUISHED BY DEPTH OR CROSS SECTION LOCATOR.               # 00021100
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # 00021200
?RETRIEVE                                                               00021305
P60   3                                                                 00021405
D 631001660930                                                          00021500
S49 09268900                                                            00021600
?PASSFILE                                                               00021705
/*                                                                      00021800
//SASSTEP  EXEC  SAS                                                    00021906
//FLOWIN DD DSN=&DATA,DISP=(OLD,PASS)                                   00022000
//SYSIN DD *                                                            00022100
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00022200
*    SAS EXAMPLE SHOWING HOW TO PERFORM FLOW ANALYSIS INVOLVING   ;     00023004
*    SIX-DAY RUNNING AVERAGES OF DATA.                            ;     00023104
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00023200
DATA FLOWOUT(KEEP=STATION  WATYEAR  MONTH  DATE  FLOW  JDAY  RUN_AVE);  00023309
INFILE FLOWIN;                                                          00023400
LENGTH DEFAULT=4;                                                       00023500
FORMAT DATE YYMMDD8.;                                                   00023600
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00023700
INPUT  @10   GAGE    $CHAR15.                                           00023800
       @37   WYR         IB2.                                           00023900
       @41   NOVAL       RB4.                                           00024000
       @1611 BMONTH      IB2.                                           00024100
             @ ;                                                        00024200
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00024300
DO MO=BMONTH TO 12, 1 TO BMONTH-1;                                      00024400
   IF WYR < 1900 THEN DELETE;                    * SKIP 1800S DATA;     00024507
   IF MO >= BMONTH  THEN YEAR = WYR-1901;   ELSE YEAR = WYR-1900;       00024600
   OFFSET = (((MO-1)*31)*4)+1;                                          00024700
   DO DAY=1 TO 31;                                                      00024800
     INPUT  @OFFSET+44 VALUE RB4. @;                                    00024900
     IF VALUE <= 999990 THEN DO;                                        00025000
       FLOW = VALUE;                                                    00025100
       DATE = MDY(MO,DAY,YEAR);                                         00025200
       JDAY = JULDATE(DATE);                                            00025300
       JDAY = JDAY - 1000 * INT(JDAY/1000);                             00025400
       DELT = 273;                                                      00025500
       IF YEAR-4*INT(YEAR/4) = 0  THEN DELT = 274;                      00025600
       IF MO > 9                                                        00025700
         THEN JDAY = JDAY - DELT;                                       00025800
         ELSE JDAY = JDAY + 92;                                         00025900
       RUN_AVE = (FLOW + LAG1(FLOW) + LAG2(FLOW)                        00026000
                       + LAG3(FLOW) + LAG4(FLOW) + LAG5(FLOW))/6;       00026100
       OUTPUT FLOWOUT;                                                  00026200
       END;                                                             00026300
     OFFSET = OFFSET + 4;                                               00026400
     END;                                                               00026500
   END;                                                                 00026600
RENAME  GAGE=STATION  WYR=WATYEAR  MO=MONTH;                            00026700
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00026800
PROC SORT; BY STATION WATYEAR DATE;                                     00026900
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00027000
PROC PRINT; BY STATION WATYEAR;                                         00027100
ID DATE;                                                                00027200
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00027300
PROC UNIVARIATE PLOT NORMAL; BY STATION; ID DATE;                       00027400
VAR FLOW RUN_AVE;                                                       00027500
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00027600
PROC PLOT NOLEGEND; BY STATION;                                         00027700
TITLE1 'RAW FLOW & SIX DAY RUNNING AVERAGES';                           00027800
TITLE2 '-';                                                             00027900
PLOT FLOW*DATE="+" RUN_AVE*DATE="M" RUN_AVE*DATE='W' RUN_AVE*DATE='E' / 00028000
OVERLAY                                                                 00028100
HAXIS='1OCT63'D '1JAN64'D '1APR64'D '1JUL64'D '1OCT64'D                 00028200
'1JAN65'D '1APR65'D '1JUL65'D '1OCT65'D                                 00028300
HREF='1OCT63'D '1OCT64'D '1OCT65'D                                      00028400
VAXIS=0 TO 200 BY 50 VREF=0 TO 200 BY 25;                               00028500
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00028600
PROC SORT; BY STATION JDAY;                                             00028700
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00028800
PROC FORMAT DDNAME=LIBRARY;                                             00028900
 VALUE JMTH 1-31='OCT' 32-61='NOV' 62-92='DEC' 93-123='JAN'             00029000
  124-151='FEB' 152-181='MAR' 182-212='APR' 213-243='MAY'               00029100
  244-273='JUN' 274-304='JUL' 305-335='AUG' 336-365='SEP';              00029200
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - ;             00029300
PROC PLOT NOLEGEND; BY STATION;                                         00029400
LABEL JDAY='DAY OF WATER YEAR';                                         00029500
FORMAT JDAY JMTH.;                                                      00029600
TITLE1 'SIX-DAY RUNNING AVERAGES';                                      00029700
TITLE2 '-';                                                             00029800
PLOT RUN_AVE*JDAY='*' /                                                 00029900
HAXIS=0 15 45 76 107 135 165 197 227 258 288 320 340                    00030000
HREF=1 31 61 92 123 151 181 212 243 273 304 335 365                     00030100
VAXIS=0 TO 200 BY 50 VREF=0 TO 200 BY 25;                               00030200
/*                                                                      00030300
