Skip to main content
OCADSAccess DataNDP-004NDP-004r1 - Table 4

Table 4

Fortran code to read and print the revised TTO data file, ndp004.rev (ndp004.for).

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