//C3STATST JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=5
/*JOBPARM SHIFT=D,LINES=90
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=5
//RAWDATA DD DSN=COINEWM.CR3CLEAN.SEQ,DISP=SHR
//CTDDATA DD DSN=COINEWM.CR3.LONGLAT(C3GOODLL),DISP=SHR
//OUTDATA DD DSN=COINEWM.C3STAVG.SEQ,DISP=SHR
//SYSIN DD *
********************************************************************;
* THIS PROGRAM READS MIDAS DATA FROM : SEQUENTIAL DATA SET,        *;
* COINEWM.CR2.SEQ                                                  *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
* BY  : H.S.MASHRIQUI.                                             *;
*                                                                  *;
********************************************************************;
*********       READ THE DATA FILE *********************************;
*********        COINEWM.CR2.SEQ   *********************************;
********************************************************************;
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   24-25    NDLATH   27-28
         NDLONGD  30-32   NDLONGM  34-35    NDLONGH  37-38
         H2OTEMP   40-47  CONDUCT   49-56   SALIN     58-65
         CHLFLUVT  67-74  AUXFLUVT  76-83   FLOSWTCH  85-92
         LITTRAN   94-101 FLORATE  103-110  WNDSPD   112-119
         RWNDSPD  121-128 CWNDSPD  130-137
                          CMWDIR   139-146  AIRTEMP  148-155
         BAROPRES 157-164 SOLRAD   166-173 ;
********************************************************************;
********************************************************************;
**************       CALCULATIONS     ******************************;
********************************************************************;
********************************************************************;
         FLAT = NDLATD  + ( NDLATM  + NDLATH/100  )/60 ;
         FLON = NDLONGD + ( NDLONGM + NDLONGH/100 )/60 ;
         RECNUM = _N_ ;
********************************************************************;
********************************************************************;
************     READ  CTD STATION DATA    *************************;
********************************************************************;
********************************************************************;
********************************************************************;
DATA B ;
   INFILE CTDDATA;
   INPUT STA$ 1-8 LATD 10-11 LATM 13-14 LATH 16-17 LONGD 19-20
         LONGM 22-23 LONGH 25-26 ;
         SLAT = LATD  + (LATM+LATH/100)/60 ;
         SLON = LONGD + (LONGM+LONGH/100)/60 ;
********************************************************************;
********************************************************************;
********************************************************************;
************     DATA SEARCH BY STATION    *************************;
********************************************************************;
********************************************************************;
********************************************************************;
DATA C ; SET B ;
      DO NUMBER = 1 TO TOTAL ;
         SET A POINT = NUMBER NOBS = TOTAL ;
         IF ( FLAT > (SLAT - 0.0055)) AND (FLAT < (SLAT + 0.0055))
            AND (FLON > (SLON - 0.0055)) AND (FLON < (SLON + 0.0055))
         THEN OUTPUT ;
      END ;
PROC SORT DATA = C ; BY STA ;
PROC MEANS DATA = C ; BY STA ;
         VAR H2OTEMP CONDUCT SALIN   CHLFLUVT LITTRAN
             CWNDSPD CMWDIR  AIRTEMP BAROPRES SOLRAD ;
         OUTPUT OUT = AA MEAN= M2OTEMP MONDUCT MALIN MHLFLUVT MITTRAN
             MWNDSPD MMWDIR  MIRTEMP MAROPRES MOLRAD ;
             PROC PRINT ;
DATA D ;
     SET AA ;
     E=',' ;
     FILE OUTDATA ;
PUT STA 1-8 E$9 M2OTEMP 10 - 16 .4 E$17 MONDUCT 18-24 .4 E$25
    MALIN 26-32 .4 E$33 MHLFLUVT 34-40 .4 E$41 MITTRAN 42-48 .4 E$49
    MWNDSPD 50-56 .4 E$57 MMWDIR 58-64 .4 E$65 MIRTEMP 66-72 .4 E$73
    MAROPRES 74-80 .4 E$81 MOLRAD 82-88 .4 ;
********************************************************************;
//
