c     program readSST
      parameter(im=720,jm=360)
c CFSv2 operational daily time series SST from seasonal forecast 
c 1. 3 runs / day from (00z,run2,run3,run4)
c 2. 130 days time serires forecast starts from 
c    the 2nd day (24hr) of the initial date (00UTC) 
c 3. grid resolution 0.5deg longitude by 0.5deg latitude
c    starts from (89.75S,-359.75E),eastward and northward
c 4. undefined value=-99999.00 
c 5. format: binary in big_endian
c 6. file name: CFSv2.ocnsst.yyyymmdd.90dy.gr
c    yyyy=the year of initial date
c    mm=the month of initial date
c    dd=the day of initial date
c 7. unit=K
c
      real*4 sst2(im,jm)
      real*4 sst3(im,jm)
      real*4 sst4(im,jm)
c
      open(10,
     & file='CFSv2.ocnsst.20150407.90dy.gr',
     & form='unformatted',
     & convert='big_endian',
     & access='direct',recl=im*jm*4)
c  
      open(20,
     & file='CFSv2.ocnsst.20150407.90dy.txt.sample',
     & form='formatted',access='sequential')
c
      write(20,*) 'grid point (-179.75E,0.25N)'
      write(20,*) 
     &'fcst days: member1(K): member2(K): member3(K):'
      write(20,*) 
     &'fcst days:   00z-run2:   00z-run3:   00z-run4:'
      do ifcst=1,130
       irec=(ifcst-1)*3+1
       read(10,rec=irec) sst2
       irec=(ifcst-1)*3+2
       read(10,rec=irec) sst3
       irec=(ifcst-1)*3+3
       read(10,rec=irec) sst4
       write(20,'(i10,3f12.2)') ifcst,sst2(361,181),sst3(361,181),
     &     sst4(361,181)
      enddo
c
      stop
      end
cccccccccccccccccccccccccccccccc