p13dat.for (File 3)
This file contains a FORTRAN 90 data-retrieval routine to read and print p13.dat (File 5). The following is a listing of this program. For additional information regarding variable definitions, variable lengths, variable types, units, and codes, please see the description for p13.dat.
c********************************************************************
c* FORTRAN 90 data retrieval routine to read and print the file
c* named "p13.dat" (File 5).
c********************************************************************
c*Defines variables*
INTEGER sta, cast, samp, bot
REAL pre, ctdtmp, ctdsal, ctdoxy, theta, sal, oxy, silca
REAL nitrat, nitrit, phspht, cfc11, cfc12, tcarb, talk
CHARACTER qualt*13
OPEN (unit=1, file='p13.dat')
OPEN (unit=2, file='p13.data')
write (2, 5)
c*Writes out column labels*
5 format (2X,'STNNBR',2X,'CASTNO',2X,'SAMPNO',2X,'BTLNBR',2X,
1 'CTDPRS',4X,'CTDTMP',4X,'CTDSAL',2X,'CTDOXY',5X,'THETA',4X,
2 'SALNTY',2X,'OXYGEN',2X,'SILCAT',2X,'NITRAT',2X,'NITRIT',2X,
3 'PHSPHT',3X,'CFC-11',3X,'CFC-12',2X,'TCARBN',2X,
4 'ALKALI',8X,'QUALT1',/,36X,'DBAR',4X,'ITS-90',4X,'PSS-78',
5 1X,'µmol/KG',4X,'ITS-90',4X,'PSS-78',1X,5('µmol/KG',1X),1X,
6 'PMOL/KG',2X,'PMOL/KG',1X,2('µmol/KG',1X),12X,'*',/,
7 25X,'*******',21X,2('*******',1X),12X,6('*******',1X),
8 1X,'*******',2X,3('*******',1X),12X,'*')
c*Sets up a loop to read and format all the data in the file*
read (1, 6)
6 format (//////////)
7 CONTINUE
read (1, 10, end=999) sta, cast, samp, bot, pre, ctdtmp,
1 ctdsal, ctdoxy, theta, sal, oxy, silca, nitrat, nitrit,
2 phspht, cfc11, cfc12, tcarb, talk, qualt
10 format (5X, I3, 7X, I1, 5X, I3, 4X, I4, 1X, F7.1, 1X, F9.4,
1 1X, F9.4, 1X, F7.1, 1X, F9.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
2 1X, F7.2, 1X, F7.2, 1X, F7.2, 1X, F8.3, 1X, F8.3, 1X, F7.1,
3 1X, F7.1, 1X, A13)
write (2, 20) sta, cast, samp, bot, pre, ctdtmp,
1 ctdsal, ctdoxy, theta, sal, oxy, silca, nitrat, nitrit,
2 phspht, cfc11, cfc12, tcarb, talk, qualt
20 format (5X, I3, 7X, I1, 5X, I3, 4X, I4, 1X, F7.1, 1X, F9.4,
1 1X, F9.4, 1X, F7.1, 1X, F9.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
2 1X, F7.2, 1X, F7.2, 1X, F7.2, 1X, F8.3, 1X, F8.3, 1X, F7.1,
3 1X, F7.1, 1X, A13)
GOTO 7
999 close(unit=1)
close(unit=2)
stop
end
