TSO ALLOC F(RAWDATA) DS('COINEWM.MIDAS.temp');
********************************************************************;
* THIS PROGRAM PRODUCES 1 CARTESIAN GRAPH USING THE SAS GPLOT      *;
* COMMAND. THIS IS AN INTERACTIVE PROGRAM WITH NO JCL. THE GRAPH   *;
* PLOTS THE VARIABLE AIRTEMP,                      FROM THE MIDAS  *;
* SEQUENTIAL DATA SET AGAINST THE NUMBER OF OBSERVATIONS.          *;
* NOTE THAT 7 VARIABLES HAVE BEEN OMITTED, THE GOPTIONS ON LINE 38 *;
* AND THE INFO ON LINE 1.                                          *;
* to concatonate files, type reuse on line 1 between ' and ;       *;
********************************************************************;
DATA P921;
   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-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 DIFVT1 366-372
         DIFVT2 374-380 SINGVT1 382-384 SINGVT2 386-392
         SINGVT3 394-400                                         ;
   observe = _N_;
GOPTIONS DEVICE=IBM3279 ROTATE=ROTATE;
TITLE1 J=C F=XSWISS 'water TEMPERATURE and salinity';
TITLE2 J=C F=XSWISS 'IN DEGREES CELCIUS';
symbol i=join;
PROC GPLOT;
   axis1 value = (f = xswiss h = 1)
         label = (f = xswiss h = 1 a=90 'water temperature')
         order = 20 to 25 by .5;
   axis2 value = (f = xswiss h = 1 )
         label = (f = xswiss h = 1 'salinity')
         order = 20 to 30 by .5;
   PLOT h2oTEMP * salin / FRAME
   VAXIS = axis1
   HAXIS = axis2;
RUN;
