//C2S5 JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=4
/*JOBPARM SHIFT=D,LINES=90
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=4
//RAWDATA DD DSN=COINEWM.CR2.SEQ,DISP=SHR
//OUTDATA DD DSN=COINEWM.C2S5.SEQ,DISP=SHR
//SYSIN DD *
GOPTIONS DEVICE=GDDMFAM4 GDDMNICKNAME=IBM3820 GDDMTOKEN=IMG240
   ROTATE ;
********************************************************************;
* THIS PROGRAM READS MIDAS DATA FROM : SEQUENTIAL DATA SET,        *;
* COINEWM.CR2.SEQ AND MAKES A GRAPH OF S1 LINES                    *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
* BY  : H.S.MASHRIQUI.                                             *;
*                                                                  *;
********************************************************************;
DATA A ;
   INFILE RAWDATA FIRSTOBS = 38972 OBS = 41791 ;
   INPUT TDDATEM 2-3 TDDATED 5-6 TDDATEY 8-11
         TDTIMEH 13-14 TDTIMEM 16-17 TDTIMES 19-20
         NDLATD 22-23 NDLATM 26-27 NDLATH 29-30
         NDLONGD 35-37 NDLONGM 40-41 NDLONGH 43-44
         H2OTEMP 260-266 CONDUCT 268-274 SALIN 276-282
         CHLFLUVT 284-290 AUXFLUVT 292-298 FLOSWTCH 300-302
         LITTRAN 304-309 FLORATE 311-315 WNDSPD 317-321
         RWNDSPD 323-327 CWNDSPD 329-333 MWDIR 335-337
         DEVWDIR 339-341 CMWDIR 343-345 AIRTEMP 347-351
         BAROPRES 353-357 SOLRAD 359-364 ;
         OBS1 = _N_ ;
         FLAT = NDLATD + ( NDLATM + NDLATH/100 )/60 ;
********************************************************************;
********************************************************************;
********     SCREEN DATA WITH GIVEN CONDITIONS    ******************;
********************************************************************;
********************************************************************;
     IF ( TDDATED = 15 AND TDTIMEH = 20 AND TDTIMEM > 34 )
          THEN DELETE ;
     IF ( NDLATD   < 25.0 OR NDLATD   > 35 ) THEN DELETE ;
********************************************************************;
     IF ( FLAT     < 28.0 OR FLAT   > 30.00 ) THEN DELETE ;
********************************************************************;
     IF ( NDLONGD  < 85.0 OR NDLONGD  > 99 ) THEN DELETE ;
     IF ( H2OTEMP  < 20.0 OR H2OTEMP  > 30 ) THEN H2OTEMP  = . ;
     IF ( CONDUCT  < 00.0 OR CONDUCT  > 10 ) THEN CONDUCT  = . ;
     IF ( SALIN    < 00.0 OR SALIN    > 36 ) THEN SALIN    = . ;
     IF ( CHLFLUVT < 0.05 OR CHLFLUVT > 1  ) THEN CHLFLUVT = . ;
     IF ( LITTRAN  < 0.00 OR LITTRAN  > 99 ) THEN LITTRAN  = . ;
     IF ( CWNDSPD  < 0.00 OR CWNDSPD  > 30 ) THEN CWNDSPD  = . ;
     IF ( CMWDIR   < 0.00 OR CMWDIR  > 360 ) THEN CMWDIR   = . ;
     IF ( AIRTEMP  < 0.00 OR AIRTEMP  > 39 ) THEN AIRTEMP  = . ;
     IF ( BAROPRES < 27.0 OR BAROPRES > 33 ) THEN BAROPRES = . ;
     IF ( SOLRAD   < 0.00 OR SOLRAD   > 1  ) THEN SOLRAD   = . ;
     D = '-' ;
     C = ',' ;
     E = ':' ;
     PROC SORT DATA = A ; BY FLAT ;
********************************************************************;
********************************************************************;
 DATA C ; SET A ;
 KEEP TDDATEM TDDATED TDDATEY TDTIMEH TDTIMEM TDTIMES FLAT
      H2OTEMP CONDUCT SALIN LITTRAN C D E ;
     PROC SORT  DATA = C ; BY FLAT ;
     PROC PRINT DATA = C ;
 DATA D ;
      SET C ;
      FILE OUTDATA ;
 PUT TDDATEM 1-2 D$3 TDDATED 4-5 D$6 TDDATEY 7-10 C$11
     TDTIMEH 13-14 E$15 TDTIMEM 16-17 E$18 TDTIMES 19-20 C$21
     FLAT 23-29 .3 C$30 H2OTEMP 32-38 .3 C$39 CONDUCT 41-47 .3
     C$48 SALIN 50-56 .3 C$57 LITTRAN 59-65 .3 C$66 ;
 RUN ;
//
