C FORTRAN RETRIEVAL CODE TO READ AND PRINT TTO DATA
REAL LATMIN,LONMIN,NIT,NITR
INTEGER ALK,TOTCO
CHARACTER*1 LATDIR,LONDIR
OPEN (unit=1, file='ndp004.rev')
OPEN (unit=2, file='ndp004.out')
C READ THE HEADER RECORD FOR EACH STATION
14 CONTINUE
READ(1,15) ILEG,ISTN,LATDEG,LATMIN,LATDIR,LONDEG,LONMIN,LONDIR,
1 IYEAR,IMTH,IDAY,IGMT,IDEPTH
15 FORMAT(1X,I1,1X,I3,1X,I2,F5.1,1X,A1,I3,F5.1,1X,A1,1X,I2,
1 1X,I2,1X,I2,1X,I4,1X,I4)
C WRITE THE HEADER RECORD INFORMATION
WRITE(2,16) ILEG,ISTN,LATDEG,LATMIN,LATDIR,LONDEG,LONMIN,
1 LONDIR,IYEAR,IMTH,IDAY,IGMT,IDEPTH
16 FORMAT(1X,'LEG=',I1,1X,'STATION # ',I3,1X,'LATITUDE ',I2,1X,F5.1,
1 1X,A1,1X,'LONGITUDE ',I3,F5.1,1X,A1,1X,'YEAR ',I2,1X,'MONTH ',I2,
2 1X,'DAY=',I2,1X,'GMT=',I4,1X,'BOTTOM DEPTH=',I4)
WRITE(2,18)
18 FORMAT(1X,'BOTTLE',1X,'PRESS.',1X,'DEPTH',1X,'TEMP',2X,'POT.',2X,
1 'SALINITY',1X,'THETA',1X,'SIGMA',2X,'SIGMA',2X,'O2',1X,
2 'SILICATE',1X,'PO4',1X,'NO3',1X,'NO2',1X,'ALKALINITY',1X,'TOTAL',
3 1X,'PCO2',3X,'TOTAL',/ 3X,'#',23X,'TEMP',17X,'2000',3X,'4000',
4 39X,'CO2',10X,'CO2')
1 CONTINUE
C READ THE TTO DATA FOR EACH STATION
READ(1,100,END=99) IBOT,IPRESS,IDEPTH,TEMP,TEMPPT,SALIN,THETA,
1 SIGMA,FOUR,IOXYG,SIL,PH,NIT,NITR,ALK,TOTCO,PCO,TOTCOR
100 FORMAT(2X,I3,1X,I4,1X,I4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
1 F6.3,1X,F6.3,1X,I3,1X,F5.1,1X,F4.2,1X,F4.1,1X,F4.2,1X,I4,1X,I4,
2 1X,F6.1,1X,F6.1)
C IF THE PROGRAM HAS REACHED THE BLANK LINE DELIMITER BETWEEN
C STATIONS, GO AND READ THE NEXT STATION HEADER RECORD
IF(IBOT .EQ. 0) GO TO 14
C WRITE THE TTO DATA FOR EACH STATION
WRITE(2,200) IBOT,IPRESS,IDEPTH,TEMP,TEMPPT,SALIN,THETA,
1 SIGMA,FOUR,IOXYG,SIL,PH,NIT,NITR,ALK,TOTCO,PCO,TOTCOR
200 FORMAT(2X,I3,3X,I4,2X,I4,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
1 F6.3,1X,F6.3,1X,I3,1X,F5.1,1X,F4.2,1X,F4.1,1X,F4.2,4X,I4,3X,I4,
2 1X,F6.1,1X,F6.1)
GO TO 1
99 CONTINUE
STOP
END
tbeaty 02/16/2001