c                                         Ishii M    May 9, 2009.

c These are text-formated files, storing historical ocean temperatures.
c The following shows a sample program for reading the data.
c In the this sample program,

c   dat(1,    1) contains S at 89.5S   0.5E,
c   dat(360,  1) contains S at 89.5S 359.5E,
c   dat(  1,180) contains S at 89.5N   0.5E, and
c   dat(360,180) contains S at 89.5N 359.5E,

c at k-th level. Variable "zlev" gives you depth of analyzed data.
c Variable "elm" indicates which the record stores temperature or
c analysis error.

      parameter (im=360,jm=180,km=28)
      dimension dat(im,jm)
      integer   idat(im)
      character*15 elm

      open(1,file='sal.1946.txt')
c
      do mon=1,12
         do ne=1,2
            do k=1,km
               read(1,'(a15,i5,i3,f6.1)') elm,iy,mo,zlev
               print '(a15,i5,i3,f6.1)',  elm,iy,mo,zlev
               do j=1,jm
                  do i=1,im,15
                     read(1,'(15i5)') (idat(ii+i),ii=0,14)
                  end do
                  do i=1,im
                     dat(i,j)=float(idat(i))*0.001
                  end do
               end do
            end do
         end do
      end do
c
      close(1)
C
      stop
      end