//C4TSTPGM 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.CRUISE4.CLEANMID,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 C$ 20 NDLATD 21-22 E$ 23 NDLATM 24-25
         F$ 26 NDLATH 27-28 C$29 NDLONGD 30-32 E$ 33
         NDLONGM 34-35 F$ 36 NDLONGH 37-38 C$ 39 H2OTEMP 40-47
         .4 C$ 48 CONDUCT 49-56 .4 C$ 57 SALIN 58-65 .4 C$ 66
         CHLFLUVT 67-74 .4 C$ 75 AUXFLUVT 76-83 .4 C$ 84
         FLOSWTCH 85-92 .4 C$93 LITTRAN 94-101 .4 C$ 102 FLORATE
         103-110 .4 C$ 111 WNDSPD 112-119 .4 C$ 120 RWNDSPD 121-128
         .4 C$ 129 CWNDSPD 130-137 .R C$ 138 CMWDIR 139-146 .4 C$
         147 AIRTEMP 148-155 .4 C$ 156 BAROPRES 157-164 .4 C$ 165
         SOLRAD 166-173 .R C$174;
         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  > 34 ) THEN H2OTEMP  = . ;
     IF ( CONDUCT  < 00.0 OR CONDUCT  > 99 ) THEN CONDUCT  = . ;
     IF ( SALIN    < 15.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  > 18 ) THEN CWNDSPD  = . ;
     IF ( CMWDIR   < 0.00 OR CMWDIR  > 360 ) THEN CMWDIR   = . ;
     IF ( AIRTEMP  < 0.00 OR AIRTEMP  > 34 ) THEN AIRTEMP  = . ;
     IF ( BAROPRES < 29.0 OR BAROPRES > 30 ) THEN BAROPRES = . ;
     IF ( SOLRAD   < 0.00 OR SOLRAD   > 1  ) THEN SOLRAD   = . ;
********************************************************************;
********************************************************************;
 TITLE1 J=C F=SWISSB 'SALINITY' ;
 TITLE2 J=C F=SWISSB 'CRUISE 4 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 ;
//
