PROGRAM READSST
c
c     Program to read the VERSION 1 OI WEEKLY, MONTHLY compocp site, and CLIMATOLOGY grids.
c     All grids are 1x1, (360,180), (longitude,latitude).
c
c     Each global grid is preceeded by a single header record.  The table
c     below shows the meaning of the header record variables for each of
c     the three SST data products.
c
c              IYRST IMST IDST       IYREND IMEND IDEND     NDAYS    INDEX
c           -----------------------  -------------------    ------    ----
c    WEEKLY (start date, yr/mn/day) (end date, yr/mn/day)      7        ??
c    MONTHLY    (yr/mn/1)              (yr/mn/day)          days/mn     ??
c    CLIM.      (99/mn/1)              (99/mn/day)          days/mn      0
c           --------------------------------------------------------------
c           Note: numeric constants denote fixed place-holder values
c   
c
c  Following the header are integer SST values in degrees C time 100, with
c  format (20I4).
c  These data are read into integer array ISST, converted to degree C, and 
c  placed into real array SST.
c
c  The geo-location of the SST array elements are:
c    SST(1,1)     = 179.5W, 89.5S
c    SST(1,2)     = 179.5W, 88.5S
c    SST(2,1)     = 178.5W, 89.5S
c    SST(360,180) = 179.5E, 89.5N
c
c  NOTES: 
c     1) all values less than or equal to -1.78 C are ice
c     2) a land/sea mask should be used to mask out OI SST analyzed values
c        not located in the ocean, e.g. program ls.f and data file ls.dat

      DIMENSION SST(360,180)
      DIMENSION ISST(360,180)
      CHARACTER FLNM*64
      IGRD=0

      WRITE(*,1000)
 1000 FORMAT('NAME OF INPUT FILE : ',$)
      READ(*,'(A64)') FLNM

      OPEN(10,FILE=FLNM,form='formatted',status='old')


c     read the header and a grid
200   READ(10,6,END=100) IYRST,IMST,IDST,IYREND,IMEND,IDEND,NDAYS,INDEX
  6   FORMAT(8I5)
      READ(10,8,END=100) ISST
  8   FORMAT(20I4)


      DO 32 I=1,360
        DO 31 J=1,180
          SST(I,J) = 0.01*FLOAT(ISST(I,J))
  31    CONTINUE
  32  CONTINUE

c     Print some data from the first 10 grids
      IGRD = IGRD + 1
      IF (IGRD.LE.10) PRINT 7,IGRD,
     * IYRST,IMST,IDST,IYREND,IMEND,IDEND,SST(70,80)
    7  FORMAT ('IGRD =',I3,3X,'DATES =',3I3,' - ',3I3,3X,
     * 'SST (110.5W,10.5S) =',F6.2)
       GO TO 200

c     Print some data from the last grid           
  100 PRINT 7,IGRD,
     * IYRST,IMST,IDST,IYREND,IMEND,IDEND,SST(70,80)
      STOP                                                      
      END