//CR5BYSTA JOB ,'MASRIQUI',MSGCLASS=S,
//          CLASS=H,NOTIFY=COINEWM,TIME=24
/*JOBPARM SHIFT=D,LINES=990
/*ROUTE PRINT IOROOM
//*
//STEP1 EXEC SAS6,TIME=24
//RAWDATA DD DSN=COINEWM.CRUISE6.MIDAS,DISP=SHR
//CTDDATA DD DSN=COINEWM.CRUISE6.STA,DISP=SHR
//OUTDATA DD DSN=COINEWM.CRUISE6.BYSTA,DISP=SHR
//SYSIN DD *
********************************************************************;
*                                                                  *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
*                                                                  *;
********************************************************************;
*********                          *********************************;
*********                          *********************************;
********************************************************************;
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 ;
********************************************************************;
********************************************************************;
**************       CALCULATIONS     ******************************;
********************************************************************;
********************************************************************;
         FLAT = NDLATD  + ( NDLATM  + NDLATH/100  )/60 ;
         FLON = NDLONGD + ( NDLONGM + NDLONGH/100 )/60 ;
         RECNUM = _N_ ;
********************************************************************;
********************************************************************;
**************   DISCARD BAD DATA     ******************************;
********************************************************************;
********************************************************************;
*    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  <  5.0 OR H2OTEMP  > 40 ) THEN H2OTEMP  = . ;
     IF ( CONDUCT  < 00.0 OR CONDUCT  > 10 ) THEN CONDUCT  = . ;
     IF ( SALIN    <  5.0 OR SALIN    > 40 ) THEN SALIN    = . ;
     IF ( CHLFLUVT < 0.01 OR CHLFLUVT > 1  ) THEN CHLFLUVT = . ;
     IF ( LITTRAN  < 0.00 OR LITTRAN  > 99 ) THEN LITTRAN  = . ;
     IF ( CWNDSPD  < 0.00 OR CWNDSPD  > 40 ) THEN CWNDSPD  = . ;
     IF ( CMWDIR   < 0.00 OR CMWDIR  > 360 ) THEN CMWDIR   = . ;
     IF ( AIRTEMP  <  4.0 OR AIRTEMP  > 45 ) THEN AIRTEMP  = . ;
     IF ( BAROPRES < 10.0 OR BAROPRES > 50 ) THEN BAROPRES = . ;
     IF ( SOLRAD   < 0.00 OR SOLRAD   > 1  ) THEN SOLRAD   = . ;
C = ',' ;
********************************************************************;
********************************************************************;
********************************************************************;
************     READ  CTD STATION DATA    *************************;
********************************************************************;
********************************************************************;
********************************************************************;
DATA B ;
   INFILE CTDDATA DELIMITER = ',' MISSOVER ;
   INPUT STA$     LATD       LATM       LONGD       LONGM ;
         SLAT = LATD  + LATM/60  ;
         SLON = LONGD + LONGM/60 ;
********************************************************************;
D = '-' ;
C = ',' ;
********************************************************************;
E = ':' ;
********************************************************************;
************     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 NOPRINT ; 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 ;
C = ',' ;
FILE OUTDATA ;
PUT STA$ 1-8 C$ 9 M2OTEMP 12-19 .3 C$  20 MONDUCT 22-29 .3
C$ 30 MALIN 32-39 .3 C$ 40 MHLFLUVT 42-49 .3 C$ 50 MITTRAN
52-59 .3 C$ 60 MWNDSPD 62-69 .3 C$ 70 MMWDIR 72-79 .3 C$ 80
MIRTEMP 82-89 .3 C$ 90 MAROPRES 92-99 .3 C$ 100 MOLRAD 102-107 .3
C$ 108 ;
//
