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
Last modified: 2021-03-17T18:30:27Z