  This retrieval uses sas to compute and plot an aggressiveness         00000040
  index,  a function of pH, alkalinity, and Calcium hardness.           00000050
  The SAS part of the code contains a good example of the               00000051
  proper use of the MORE=4 output to pass Secondary Station             00000060
  Numbers and other information from the STORET station                 00000070
  header records of the FCF file into SAS.                              00000080
                                                                        00000090
PGM=MEAN,MORE=4,PRT=NO,PURP=TEST,A=1110NET,S=230040,S=230041,           00000100
P=400,PM=403,P=410,P=900,PM=916*2.5,                                    00000110
BD=61,ED=62,R=*,ANC=A,DSROC=1,BOT=W,                                    00000120
DG=Q,SF=MN,PCTL=50,BASE=0101,                                           00000130
SASPARMS=BEGIN,                                                         00000140
DATA AI;                                                                00000150
INFILE FCF LENGTH=L;                                                    00000160
FORMAT DATE YYMMDD8.;                                                   00000170
FORMAT TIME HHMM5.;                                                     00000180
FORMAT DEPTH $6.;                                                       00000190
LENGTH SEC_STN $ 12; RETAIN SEC_STN ' ';                                00000200
LENGTH LAT_LONG $ 25; RETAIN LAT_LONG ' ';                              00000210
LENGTH LOCATION $ 48; RETAIN LOCATION ' ';                              00000220
LENGTH AGENCY $ 8;                                                      00000230
LENGTH STATION $ 15;                                                    00000240
ARRAY P(I) P1-P50;                                                      00000250
ARRAY R(I) $1 R1-R50;                                                   00000260
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000270
INPUT @26 YYDLIM $2. @ ;                                                00000280
IF L=145 THEN DO;                                                       00000290
  INPUT @25 LINE $1. @; IF LINE='1' THEN INPUT @108 SEC_STN $12. @;     00000300
                        IF LINE='2' THEN INPUT @91 LAT_LONG $22. @;     00000310
                        IF LINE='3' THEN INPUT @91 LOCATION $48. @;     00000320
  INPUT;                                                                00000330
  END;                                                                  00000340
IF L=305 | L=75; * OMITS PARAMETER HEADERS AND STATION HEADERS ;        00000350
IF YYDLIM^='99'; * OMITS DELIMITER RECORDS ;                            00000360
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000370
INPUT @1 AGENCY $8. @9 STATION $15. @26 DATE YYMMDD6.                   00000380
@32 KHR $2. KMN $2.                                                     00000390
@36 (P1-P50)(RB4.)                                                      00000400
@236 (R1-R50)($1.) @300 DEPTH $6.;                                      00000410
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000420
IF KHR<'25' THEN TIME=HMS(KHR,KMN,0);                                   00000430
DO OVER P;  IF P > 0. & P < 1.E-15 THEN P=. ;                           00000440
 END;   * END OF DO OVER P LOOP ;                                       00000450
IF P2*P3>0 THEN P4=P1+LOG10(P2*P3);                                     00000460
IF SEC_STN=' ' THEN SEC_STN='NONE SPEC.';                               00000470
RENAME P1=P_H P2=TOT_ALK P3=CA_HARD P4=AGGR_IND;                        00000480
IF P4^=.;                                                               00000490
DROP KHR KMN YYDLIM I P5-P50 R1-R50 LINE;                               00000500
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000510
PROC SORT; BY AGENCY STATION DATE TIME;                                 00000520
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000530
PROC PRINT PAGE=2;ID DATE;                                              00000540
BY AGENCY STATION SEC_STN LAT_LONG LOCATION;                            00000550
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000560
PROC PLOT;                                                              00000570
BY AGENCY STATION SEC_STN LAT_LONG LOCATION;                            00000580
TITLE 'AGGRESSIVE INDEX AND P_H';                                       00000590
PLOT AGGR_IND*DATE='#' P_H*DATE='*'/OVERLAY                             00000600
HAXIS=366 456 547 639 731 821 912 1004 1096                             00000610
HREF=456 547 639 731 821 912 1004 1096                                  00000620
VAXIS=6 TO 16 BY 2                                                      00000630
VREF=8 TO 16 BY 2;                                                      00000640
PLOT AGGR_IND*P_H='*';                                                  00000650
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - -;            00000660
STOPSAS,                                                                00000670
./iii      JOB (aaaaSTORPuuu,Miii),STORET,NOTIFY=iii,TIME=(,30),        00000680
./             MSGLEVEL=(1,1),PRTY=4                                    00000690
**ROUTE  PRINT HOLD                                                     00000700
**ROUTE  PRINT RMT13                                                    00000710
**JOBPARM LINES=10                                                      00000720
