a5dat.for (File 3)
This file contains a FORTRAN 77 program to read and print a5.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 a5.dat.
c****************************************************************
c* FORTRAN 77 data retrieval routine to read and print the
c* file named "a5.dat" (File 5).
c****************************************************************
CHARACTER qualt*12
INTEGER sta, cast, samp, bot
REAL pre, ctdtmp, ctdsal, ctdoxy, theta, sal, oxy, silca
REAL nitrat, nitrit, phspht, talk, tcarb, ph
OPEN (unit=1, file='a5.dat')
OPEN (unit=2, file='a5data.dat')
write (2, 5)
c*Writes out column labels*
5 format (2X,'STNNBR',2X,'CASTNO',2X,'SAMPNO',2X,'BTLNBR',2X,
1 'CTDPRS',2X,'CTDTMP',2X,'CTDSAL',3X,'CTDOXY',3X,'THETA',4X,
2 'SALNTY',2X,'OXYGEN',2X,'SILCAT',2X,'NITRAT',2X,'NITRIT',2X,
3 'PHSPHT',2X,'ALKALI',2X,'TCARBN',7X,'PH',7X,'QUALT1',/,
5 36X,'DBAR',2X,'ITS-90',2X,'PSS-78',2X,'µmol/KG',2X,'ITS-90',
6 4X,'PSS-78',1X,'µmol/KG',1X,6('µmol/KG',1X),1X,'@25_DEG',12X,
8 '*',/,25X,'*******',17X,2('*******',2X),9X,'*******',1X,
8 7('*******',1X),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, talk, tcarb, ph, qualt
10 format (5X, I3, 7X, I1, 6X, I2, 6X, I2, 1X, F7.1, 1X, F7.4,
1 1X, F7.4, 1X, F8.3, 1X, F7.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
2 1X, F7.2, 1X, F7.2, 1X, F7.2, 1X, F7.1, 1X, F7.1, 1X, F8.3,
3 1X, A12)
write (2, 20) sta, cast, samp, bot, pre, ctdtmp,
1 ctdsal, ctdoxy, theta, sal, oxy, silca, nitrat, nitrit,
2 phspht, talk, tcarb, ph, qualt
20 format (5X, I3, 7X, I1, 6X, I2, 6X, I2, 1X, F7.1, 1X, F7.4,
1 1X, F7.4, 1X, F8.3, 1X, F7.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
2 1X, F7.2, 1X, F7.2, 1X, F7.2, 1X, F7.1, 1X, F7.1, 1X, F8.3,
3 1X, A12)
GOTO 7
999 close(unit=1)
close(unit=2)
stop
end