//MASRIQUI JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=25
/*JOBPARM SHIFT=D
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=25
//RAWDATA DD DSN=COINEWM.CR2.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 TIME VS H2OTEMP.            *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
* BY  : H.S.MASHRIQUI.                                             *;
*                                                                  *;
********************************************************************;
DATA A ;
   INFILE RAWDATA;
   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 ;
         OBS = _N_ ;
         MTIME = (TDDATED * 24 * 60 + TDTIMEH * 60 + TDTIMEM
                 + TDTIMES / 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 ( 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   = . ;
********************************************************************;
********************************************************************;
 TITLE1 J=C F=SWISSB 'CHLFLUVT AND SALINITY' ;
 TITLE2 J=C F=SWISSB 'CRUISE 2 DATA';
 SYMBOL H=0.25   ;
 PROC GPLOT ;
 AXIS1 VALUE  = ( F = SWISSB H=1)
       LENGTH =   75 PCT
       LABEL  = ( F = SWISS H=1 A=90 'CHLFLUVT' )
       ORDER  = 0.10 TO 1.0 BY 0.1 ;
 AXIS2 VALUE  = ( F = SWISSB H = 1)
       LENGTH =   75 PCT
       LABEL  = ( F = SWISS H=1  'SALINITY PARTS PER THOUSAND' )
       ORDER  = 15 TO 35 BY 5 ;
       PLOT CHLFLUVT * SALIN/FRAME
       HAXIS = AXIS2
       VAXIS = AXIS1 ;
 RUN ;
//
