Example of Fortran program to read CTD data from the R/V Le Noroit

 ####  1rst record : Cruise header
       Number of stations recorded in the file, 
     Pressure resolution (dbar), Text of 50 characters

      [nbsta,dz,text] in the fortran program.

    Example Cruise header red in ctdpo01.a file :
   80   2.00  COARE-POI, 156E, 2/12/92 au 28/02/93
   with
      80   = Number of stations in this file
      2.00 = Pressure resolution. 
      COARE-POI, 156E, 2/12/92 au 28/02/93 = Cruise comment text     

 #### Followed by one record per station :
       (i.e. 80 station records stored in ctdpo01.a
             85 station records in ctdpo01.b)

       Station number, Number of physical parameters,
      Levels number,Day, Month, Year, Hour(GMT), Minutes,
      Latitude,Longitude,(VALUE(k,j),k=1,nblev),j=1,npara)

      [ is,npar,nbniv,day,month,year,hour,imin,xlat,xlon,
        ((val(k,j),k=1,nbniv),j=1,npar)] in the fortran program

    Example station #1 (the first station of the ctdpo01.a file)
     
     is =1   station number
     npar =  6  parameters are measured at each level
     nbniv = 504  levels
     day,month,year,hour,imin = 2/12 /92   21h32 
     Latitude = -20.00  (20S)
     Longitude= 163.23  (163.23E)
     val = matix of the parameters values ( 504 lines*6 columns)

     Parameters occurrence (for i=1 to 504 levels) :
     val(i,1) = Pressure (dbar)
     val(i,2) = Temperature (deg.C)
     val(i,3) = Salinity (PSU)
     val(i,4) = Sigmateta (kg/m3)
     val(i,5) = Dynamic Height / 0 Dbar (Dyn. m.)
     val(i,6) = Oxygen (ml/l)
    missing data to 9999.

     The 10 first levels of the station #1 (val(k,j),k=1,10),j=1,6)
  Pressure   Temp.    Sal.   Sigmateta   Dynamic   Oxygen
                                          Height
     0.0    25.877    35.673  23.58127     0.000  9999.000
     4.0    25.877    35.673  23.58154     0.017  9999.000
     6.0    25.770    35.663  23.60804     0.026  9999.000
     8.0    25.610    35.667  23.66033     0.034  9999.000
    10.0    25.581    35.662  23.66576     0.043  9999.000
    12.0    25.510    35.662  23.68800     0.051  9999.000
    14.0    25.503    35.663  23.69123     0.060  9999.000
    16.0    25.498    35.664  23.69305     0.068  9999.000
    18.0    25.498    35.663  23.69254     0.076  9999.000
    20.0    25.482    35.662  23.69688     0.085  9999.000

#### The following fortran program read the data of ctdpo1.a file
 (or ctdpo01.b) and print them on standard output.
      

	parameter(NPARA=6,NBLEV=1500)
	real val(NBLEV,NPARA),xlat,xlon,dz
        integer*4 is,npara,nbniv,day,month,year,hour,imin
        integer*4 nbsta
	character*50 text
        data name/'ctdpo01.a'/
c       data name/'ctdpo01.b'/
c
	open(unit=22,file=name,form='un',status='old')
c
c................ read the cruise header
c
	read(22)nbsta,dz,text
        write(6,*)nbsta,dz,text

c
c.................read the data of each station
c
	do ista=1,nbsta
	   read(22)is,npar,nbniv,day,month,year,hour,imin,
     &     xlat,xlon,((val(k,j),k=1,nbniv),j=1,npar)
           write(6,*)is,npara,nbniv,day,month,year,hour,imin,
     &     xlat,xlon
c
           do k=1,nbniv
c.................Print the kth level parameters
             write(6,*)(val(k,j),j=1,npar)
           enddo
c
	enddo
        close(22)
	stop 000
	end