//MASRIQUI JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=2
/*JOBPARM SHIFT=D,LINES=90
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=2
//RAWDATA DD DSN=COINEWM.CR2.SEQ,DISP=SHR
//SYSIN DD *
********************************************************************;
* THIS PROGRAM READS MIDAS DATA FROM : SEQUENTIAL DATA SET,        *;
* COINEWM.CR2.SEQ                                                  *;
* THIS PROGRAM CALCULATES AND PRINTS MEAN OF SOME VARABLES         *;
*                                                                  *;
*                                                                  *;
* 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 ;
         FLAT = NDLATD  + ( NDLATM  + NDLATH/100  )/60 ;
         FLON = NDLONGD + ( NDLONGM + NDLONGH/100 )/60 ;
     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.07 OR CHLFLUVT > 1  ) THEN CHLFLUVT = . ;
     IF ( AUXFLUVT < 9.00 OR AUXFLUVT > 10 ) THEN AUXFLUVT = . ;
     IF ( LITTRAN  < 0.00 OR LITTRAN  > 99 ) THEN LITTRAN  = . ;
     IF ( FLORATE  < 20.0 OR FLORATE  > 50 ) THEN FLORATE  = . ;
     IF ( CWNDSPD  < 0.00 OR CWNDSPD  > 30 ) THEN CWNDSPD  = . ;
     IF ( CMWDIR   < 0.00 OR CMWDIR  > 360 ) THEN CMWDIR   = . ;
     IF ( AIRTEMP  < 20.0 OR AIRTEMP  > 30 ) THEN AIRTEMP  = . ;
     IF ( BAROPRES < 28.0 OR BAROPRES > 31 ) THEN BAROPRES = . ;
     IF ( SOLRAD   < 0.00 OR SOLRAD   > 20 ) THEN SOLRAD   = . ;
PROC SORT  DATA = A ; BY TDDATED ;
PROC MEANS DATA = A ;
     VAR   H2OTEMP CONDUCT SALIN CHLFLUVT AUXFLUVT LITTRAN FLORATE
     CWNDSPD CMWDIR AIRTEMP BAROPRES SOLRAD ;
PROC CHART DATA = A ;
VBAR H2OTEMP / TYPE = FREQ MIDPOINTS = 22.0 22.5 23.0 23.5
         24.0 24.5 25.0 25.5 26.0 26.5 27.0 27.5 ;
PROC CHART DATA = A ;
VBAR CONDUCT / TYPE = FREQ MIDPOINTS = 3.0 3.25 3.50 3.75
         4.00 4.25 4.50 4.75 5.00 5.25 5.50 ;
PROC CHART DATA = A ;
VBAR SALIN   / TYPE = FREQ MIDPOINTS = 19.0 20.0 21.0 22.0
         23.0 24.0 25.0 26.0 27.0 28.0 29.0 30.0 31.0 32.0
         33.0 34.0 35.0 ;
PROC CHART DATA = A ;
VBAR CHLFLUVT / TYPE = FREQ MIDPOINTS = 0.0 0.25 0.50 0.75
         1.00 1.25 6.00 7.00 8.00  ;
PROC CHART DATA = A ;
VBAR AUXFLUVT / TYPE = FREQ MIDPOINTS = 3.00 4.00 5.00 6.00
         7.00 8.00 9.00 10.0  ;
PROC CHART DATA = A ;
VBAR LITTRAN / TYPE = FREQ MIDPOINTS = 10.0 20.0 30.0 40.0
         50.0 60.0 70.0 80.0 90.0 100.0 ;
PROC CHART DATA = A ;
VBAR FLORATE / TYPE = FREQ MIDPOINTS = 5.0 10.0 15.0 20.0 25.0
         30.0 35.0 40.0 45.0 50.0 60.0 65.0 70.0 ;
PROC CHART DATA = A ;
VBAR CWNDSPD  / TYPE = FREQ MIDPOINTS = 5.0 10.0 15.0 20.0
         25.0 30.0 35.0 40.0 45.0 50.0 55.0 60.0 ;
PROC CHART DATA = A ;
VBAR CMWDIR   / TYPE = FREQ MIDPOINTS = 0.00 100  200  300
         400 500 600 700 800 900 ;
PROC CHART DATA = A ;
VBAR AIRTEMP / TYPE = FREQ MIDPOINTS = 0.00 5.00 10.0 15.0
         20.0 25.0 30.0  ;
PROC CHART DATA = A ;
VBAR BAROPRES / TYPE = FREQ MIDPOINTS = 5.0 10.0 15.0 20.0 25.0
         30.0 35.0 ;
PROC CHART DATA = A ;
VBAR SOLRAD   / TYPE = FREQ MIDPOINTS = 0.0 0.25 0.50 0.75
         1.00 10.0 20.0 30.0 40.0 50.0 55.0 60.0 65.0 70.0 75.0
         80.0 85.0 ;
//
