c********************************************************************
c* FORTRAN 90 data retrieval routine to read and print the file     *           
c* named "i6s.dat"                                                  *
c********************************************************************
       CHARACTER samp*4, qualt*13
       INTEGER sta, cast, bot 
       REAL pre, ctdtmp, ctdsal, ctdoxy, theta, sal, oxy, silca
       REAL nitrat, phspht, cfc11, cfc12, tcarb, alkali 
       OPEN (unit=1, file='i6s.dat')
       OPEN (unit=2, file='i6s.data')
       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',2X,'CTDOXY',3X,'THETA',4X,
     2 'SALNTY',2X,'OXYGEN',2X,'SILCAT',2X,'NITRAT',2X,
     3 'PHSPHT',3X,'CFC-11',3X,'CFC-12',2X,'TCARBN',2X,'ALKALI',8X,
     4 'QUALT1',/,36X,'DBAR',2X,'ITS-90',2X,'PSS-78',1X,'UMOL/KG',
     5 2X,'ITS_90',4X,'PSS-78',1X,4('UMOL/KG',1X),1X,'PMOL/KG',2X,
     6 'PMOL/KG',1X,2('UMOL/KG',1X),12X,'*',/,25X,'*******',17X,
     7 2('*******',1X),10X,5('*******',1X),1X,'*******',2X,
     8 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, phspht, 
     2 cfc11, cfc12, tcarb, alkali, qualt

 10    format (5X, I3, 7X, I1, 4X, A4, 6X, I2, 1X, F7.1, 1X, F7.4,
     1 1X, F7.4, 1X, F7.1, 1X, F7.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
     2 1X, F7.2, 1X, F7.2, 1X, F8.3, 1X, F8.3, 1X, F7.1, 1X, F7.1, 
     3 1X, A13)

       write (2, 20) sta, cast, samp, bot, pre, ctdtmp,
     1 ctdsal, ctdoxy, theta, sal, oxy, silca, nitrat, 
     2 phspht, cfc11, cfc12, tcarb, alkali, qualt

 20    format (5X, I3, 7X, I1, 4X, A4, 6X, I2, 1X, F7.1, 1X, F7.4,
     1 1X, F7.4, 1X, F7.1, 1X, F7.4, 1X, F9.4, 1X, F7.1, 1X, F7.2,
     2 1X, F7.2, 1X, F7.2, 1X, F8.3, 1X, F8.3, 1X, F7.1, 1X, F7.1, 
     3 1X, A13)
     

       GOTO 7
 999   close(unit=1)    
       close(unit=2)
       stop
       end