//MASRIQUI JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=4
/*JOBPARM SHIFT=D,LINES=40
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=4
//RAWDATA DD DSN=COINEWM.MIDAS.SEQ,DISP=SHR
//STATION DD DSN=COINEWM.DANA.PGMII(CTDTIME),DISP=SHR
//SYSIN DD *
********************************************************************;
* THIS PROGRAM READS MIDAS DATA FROM : SEQUENTIAL DATA SET,        *;
* COINEWM.NIDAS.SEQ                                                *;
* THIS PROGRAM CALCULATES AND PRINTS MEAN OF SOME VARABLES OF      *;
* EVERY STATION PRESENT IN 'COINEWM.DANA.PGMII(CTDTIME).           *;
*                                                                  *;
* 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
         NDTD1 49-55 NDTD2 57-63 NDTD3 65-71 NDTD4 73-79
         NDTD5 81-87
         BLKSTAT $ 89 CYCSTAT $ 91 SNRSTAT $ 93
         GDOP 95-99 TRUCOG 101-107 MAGCOG 109-115 KNOTSOG 117-120

         ERRNM 144-147 ERRDIR $ 149
         WPLATD 151-152 WPLATM 155-156 WPLATH 158-159
         WPLONGD 164-166 WPLONGM 169-170 WPLONGH 172-173
         RNGWPT 178-182 BEARWPT 184-186
         LORLATD 188-189 LORLATM 192-193 LORLATH 195-196
         LORLONGD 201-203 LORLONGM 206-207 LORLONGH 209-210
         PRETD1 215-221 PRETD2 223-229
         GPSLATD 231-232 GPSLATM 235-236 GPSLATH 238-240
         GPSLONGD 245-247 GPSLONGM 250-251 GPSLONGH 253-255
         H2OTEMP 260-266 CONDUCT 268-274 SALIN 276-280
         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 DIFVT1 366-372
         DIFVT2 374-380 SINGVT1 382-384 SINGVT2 386-392
         SINGVT3 394-400 PULS1 402-405 PULS2 407-410 SERA 412-414;
         FLAT = NDLATD  + ( NDLATM  + NDLATH/100  )/60 ;
         FLON = NDLONGD + ( NDLONGM + NDLONGH/100 )/60 ;
DATA B ;
   INFILE STATION;
   INPUT  STA $ 1-7 LATD 9-10 LATM 12-16 LONGD 18-19 LONGM 21-25
   MMONTH 27 MDAY 29-30 MYEAR 32-35 MTIME 37-40 ;
        SLAT = LATD+(LATM)/60 ;
        SLON = LONGD+(LONGM)/60 ;
DATA C ; SET B ;
     DO NUMBER = 1 TO TOTAL   ;
        SET A POINT = NUMBER NOBS = TOTAL ;
        IF ( FLAT > (SLAT - 0.003) ) & ( FLAT < (SLAT+0.003) ) &
           ( FLON > (SLON - 0.003) ) & ( FLON < (SLON+0.003) ) THEN
             OUTPUT ;
     END ;
PROC SORT  DATA = C ; BY STA ;
PROC MEANS DATA = C ; BY STA ;
     VAR   H2OTEMP CONDUCT SALIN CHLFLUVT AUXFLUVT LITTRAN FLORATE
     CWNDSPD CMWDIR AIRTEMP BAROPRES SOLRAD ;
     OUTPUT OUT = RESULT  MEAN =
     MH2OTEM MCONDCT MSALIN MCHLFLVT MAUXFLVT MLITRAN MFLORAT
     MCWSPD MMWDIR MAIRTMP MBAROPRS MSOLRAD ;
PROC PRINT ;
//
