//CR3PTIME JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=125
/*JOBPARM SHIFT=D
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=125
//RAWDATA DD DSN=COINEWM.CRUISE5.MIDAS,DISP=SHR
//SYSIN DD *
GOPTIONS DEVICE=GDDMFAM4 GDDMNICKNAME=IBM3820 GDDMTOKEN=IMG240
   ROTATE ;
********************************************************************;
* THIS PROGRAM READS MIDAS DATA FROM : SEQUENTIAL DATA SET,        *;
* COINEWM.CRUISE3.MIDAS AND MAKES A GRAPH OF TIME VS A CHOSEN       *;
* PARAMETER (AIRTEMP, H2OTEMP, SOLRAD, ETC).                       *;
*                                                                  *;
*                                                                  *;
* BY  : H.S.MASHRIQUI.                                             *;
*                                                                  *;
********************************************************************;
DATA A ;
   INFILE RAWDATA;
   INPUT TDDATEM 1-2 TDDATED 4-5 TDDATEY 7-10
         TDTIMEH 12-13 TDTIMEM 15-16 TDTIMES 18-19
         NDLATD 21-22 NDLATM 25-26 NDLATH 28-29
         NDLONGD 34-36 NDLONGM 39-40 NDLONGH 42-43
         H2OTEMP 259-265 CONDUCT 267-273 SALIN 275-281
         CHLFLUVT 283-289 AUXFLUVT 291-297 FLOSWTCH 299-301
         LITTRAN 303-308 FLORATE 310-314 WNDSPD 316-320
         RWNDSPD 322-326 CWNDSPD 328-332 MWDIR 334-336
         DEVWDIR 338-340 CMWDIR 342-344 AIRTEMP 346-351
         BAROPRES 353-357 SOLRAD 359-365 ;
         OBS = _N_ ;
         FLAT = NDLATD +(NDLATM + NDLATS/100)/60 ;
         MTIME = TDDATED * 24 * 60 + TDTIMEH * 60 + TDTIMEM
                 + TDTIMES / 60 -18700 ;
********************************************************************;
********************************************************************;
********     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  < 15.0 OR H2OTEMP  > 23 ) THEN H2OTEMP  = . ;
     IF ( CONDUCT  < 00.0 OR CONDUCT  > 99 ) THEN CONDUCT  = . ;
     IF ( SALIN    < 15.0 OR SALIN    > 33 ) 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  > 20 ) THEN CWNDSPD  = . ;
     IF ( CMWDIR   < 0.00 OR CMWDIR  > 360 ) THEN CMWDIR   = . ;
     IF ( AIRTEMP  < 0.00 OR AIRTEMP  > 25 ) THEN AIRTEMP  = . ;
     IF ( BAROPRES < 29.0 OR BAROPRES > 31 ) THEN BAROPRES = . ;
     IF ( SOLRAD   < 0.00 OR SOLRAD   > 1  ) THEN SOLRAD   = . ;
********************************************************************;
********************************************************************;
 TITLE1 J=C F=SWISSB 'SALINITY' ;
 TITLE2 J=C F=SWISSB 'CRUISE 3 DATA  ' ;
 SYMBOL I= NEEDLE ;
 PROC GPLOT ;
 AXIS1 VALUE  = ( F = SWISSB H=1)
       LENGTH =   75 PCT
       LABEL  = ( F = SWISS H=1 A=90 'SALIN' ) ;
 AXIS2 VALUE  = ( F = SWISSB H = 1)
       LENGTH =   75 PCT
       LABEL  = ( F = SWISS H=1  'TIME MINUTES' );
       PLOT SALIN    * MTIME/FRAME
       HAXIS = AXIS2
       VAXIS = AXIS1 ;
 RUN ;
//
