*  THIS DATASET IS NAMED 'STORET.HELP.EXAMPLES.SAS.LIBRARY(THREED)'     00001000
*                                                                       00003000
*  THIS FILE CONTAINS A STORET/SAS RUN-STREAM WHICH READS LAT/LONG      00004000
*  AND PARAMETRIC DATA FROM A MORE=4 FCF, PLACING IT IN A PERMANENT     00005000
*  SAS DATA SET.                                                        00005100
*  ALSO GIVEN HERE ARE THE SAS INSTRUCTIONS NEEDED IN AN INTERACTIVE    00005200
*  TSO SESSION TO RE-ACQUIRE THE DATA FROM THIS PERMANENT SAS DATASET   00005300
*  AND TO CONVERT IT TO A SAS G3D PLOT.                                 00005400
*                                                                       00006000
** DELETE THIS LINE AND ALL PRECEEDING LINES BEFORE USING THIS FILE **  00007000
                                                                        00008000
PGM=RET,MORE=4,PRT=NO,                                                  00010000
NOECHO,                                                                 00020000
STC=51,CO=139,EXTRACT=AMBNT AND (WELL OR SPRING),                       00030000
P=400,P=410,                                                            00040000
##################################################                      00050000
SASPARMS=BEGIN,                                                         00060000
OPTIONS NOSOURCE;                                                       00070000
DATA FROMFCF;                                                           00080000
LENGTH STATION $ 15;                                                    00090000
RETAIN LAT LONG;                                                        00093000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00140000
FORMAT DATE YYMMDD8.;                                                   00150000
INFILE FCF LENGTH=L;                                                    00160000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00441000
INPUT @26 YYDLIM $2. @ ; DROP YYDLIM;                                   00442000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00443000
IF L=145 THEN DO;  *  L=145 SIGNALS POTENTIAL STATION HDR RECORD      ; 00444000
  INPUT @25 LINE $1. @; DROP LINE;                                      00445000
      * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00446000
       IF LINE='2' THEN DO;                                             00447000
                        INPUT @91  LATDEG  2.                           00448000
                              @94  LATMIN  2.                           00449000
                              @97  LATSEC  4.1                          00449100
                              @102 LONGDEG 3.                           00449200
                              @106 LONGMIN 2.                           00449300
                              @109 LONGSEC 4.1                          00449400
                              @;                                        00449500
                        LAT=LATDEG+LATMIN/60+LATSEC/3600;               00449600
                        LONG=LONGDEG+LONGMIN/60+LONGSEC/3600;           00449700
                        DROP LATDEG LATMIN LATSEC;                      00449800
                        DROP LONGDEG LONGMIN LONGSEC;                   00449900
                        END;                                            00450000
      * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00450100
  END;                                                                  00450800
IF L=305 | L=75; * OMITS PARAMETER HEADERS AND STATION HEADERS ;        00457000
IF YYDLIM^='99'; * OMITS DELIMITER RECORDS ;                            00460000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00470000
INPUT @9    STATION       $15.                                          00480000
      @26   DATE      YYMMDD6.                                          00490000
      @36   LAB_PH        RB4.                                          00500000
            TOT_ALK       RB4.                                          00510000
              ;                                                         00520000
    IF LAB_PH  > 0 AND LAB_PH  < 1.E-15 THEN LAB_PH  = . ;              00521000
    IF TOT_ALK > 0 AND TOT_ALK < 1.E-15 THEN TOT_ALK = . ;              00522000
    IF LAB_PH  > 14 THEN LAB_PH = . ;                                   00523000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00530000
PROC SORT DATA=FROMFCF; BY LAT LONG;                                    00530100
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00531000
PROC UNIVARIATE DATA=FROMFCF NOPRINT; BY LAT LONG;                      00540000
VAR LAB_PH TOT_ALK;                                                     00540100
OUTPUT OUT=KEEP.MAPPABLE MEAN=LAB_PH TOT_ALK;                           00541000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00542000
TITLE1 'GROUNDWATER DATA FROM PAGE COUNTY VIRGINIA';                    00543000
TITLE2 '- * -';                                                         00544000
PROC PRINT UNIFORM DATA=KEEP.MAPPABLE;                                  00566000
* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ; 00567000
STOPSAS,                                                                00570000
./III      JOB (AAAASTORP,MIII),'SAS-THREED',                           00580001
./             MSGLEVEL=(1,1),PRTY=4                                    00590000
**ROUTE  PRINT HOLD                                                     00600000
**JOBPARM LINES=10                                                      00610000
./KEEP DD DISP=(NEW,CATLG),UNIT=DISK,DSN=IIIAAAA.SAMPLE.SASKEEP,        00620001
./   SPACE=(TRK,(5,1))                                                  00630000
                                                                        00640001
                                                                        00650001
                                                                        00660001
                                                                        00670001
                                                                        00680001
                                                                        00690001
                                                                        00700001
                                                                        00710001
                                                                        00720001
                                                                        00730001
                                                                        00740001
   WHAT FOLLOWS IS A SAMPLE TSO DATASET WITH THE SAS INSTRUCTIONS       00750001
   TO READ DATA FROM THE TSO FILE SAMPLE.SASKEEP, WHICH CONTAINS        00760001
   A SAS MATRIX NAMED 'MAPPABLE'.  IN THE TSO ENVIRONMENT, THIS         00770001
   DATASET WILL CONTINUE TO BE REFERRED TO VIA DDNAME 'KEEP', SO AS     00780001
   TO BE CONSISTENT WITH THE BATCH RUN-STREAM ABOVE.                    00790001
                                                                        00800001
   THESE TSO INSTRUCTIONS MIGHT BE CREATED AND SAVED IN YOUR LIBRARY    00801001
   WITH DATA SET NAME 'SAS.INPUT'.                                      00802001
                                                                        00803001
   ASSUMING THEY WERE, YOU WOULD INVOKE THEM USING TSO COMMAND "SAS"    00804001
   AS FOLLOWS:                                                          00805001
            .                                                           00806001
            .                                                           00806101
            READY                                                       00807001
            SAS IN(SAS.INPUT)                                           00808001
            .                                                           00809001
            (SAS DISPLAYS FOLLOW)                                       00809101
            .                                                           00809201
            .                                                           00809301
            READY                                                       00809401
                                                                        00809501
                                                                        00809601
   CONTENTS OF DATASET 'SAS.INPUT' FOLLOW:                              00809701
                                                                        00809801
TSO ALLOC FI(KEEP) DA('IIIAAAA.SAMPLE.SASKEEP') SHR REU ;               00810001
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;        00820001
DATA FIXED; SET KEEP.MAPPABLE;                                          00830001
LENGTH COLORVAL $ 8;                                                    00840001
SIZEVAL=1.0; COLORVAL='RED';                                            00850001
IF _N_=1 THEN DO;                                                       00860001
OUTPUT;                                                                 00870001
*  THESE LINES GENERATE SOME "FAKE" DATA AT THE DESIRED        ;        00871001
*  CORNER-POINTS OF THE LAT/LONG WINDOW DESIRED.               ;        00872001
*  THIS METHOD OF CONTROLLING SCALE IS NOT REQUIRED.           ;        00873001
SIZEVAL=0; COLORVAL='CYAN';                                             00880001
LAT=38.4; LONG=78.3; TOT_ALK=0; OUTPUT;                                 00890001
LAT=38.4; LONG=78.7; TOT_ALK=0; OUTPUT;                                 00900001
LAT=38.8; LONG=78.3; TOT_ALK=0; OUTPUT;                                 00910001
LAT=38.8; LONG=78.7; TOT_ALK=0; OUTPUT;                                 00920001
END;                                                                    00930001
ELSE OUTPUT;                                                            00940001
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;        00950001
GOPTIONS DEVICE=TEK4107T VPOS=43 HPOS=85                                00960001
         CTITLE=WHITE                                                   00970001
         COLORS=(BLUE CYAN GREEN YELLOW RED PINK WHITE)                 00980001
          ;                                                             00990001
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;        01000001
PROC G3D DATA=FIXED;                                                    01010001
TITLE1 'ALKALINITY IN GROUNDWATER';                                     01020001
TITLE3 '- * -';                                                         01030001
TITLE4 'PAGE COUNTY            VIRGINIA';                               01040001
FOOTNOTE1 .J=R 'STORET SYSTEM';                                         01050001
SCATTER LONG*LAT=TOT_ALK/XTICKNUM=5 YTICKNUM=5 ZTICKNUM=7               01060001
COLOR=COLORVAL SIZE=SIZEVAL                                             01070001
ZMIN=0 ZMAX=300;                                                        01080001
