c     program readSST
      parameter(im=720,jm=360)
c CFSv2 operational daily time series SST from seasonal forecast 
c 1. 9 runs / day from
c    06z,run2,run3,run4;12z,run2,run3,run4;18z,run2,run3,run4;
c 2. 45 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.45dy.gr',
     & form='unformatted',
     & convert='big_endian',
     & access='direct',recl=im*jm*4)
c  
      open(20,
     & file='CFSv2.ocnsst.20150407.45dy.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:   06z-run2:   06z-run3:   06z-run4:'
      do ifcst=1,45
       irec=(ifcst-1)*9+1
       read(10,rec=irec) sst2
       irec=(ifcst-1)*9+2
       read(10,rec=irec) sst3
       irec=(ifcst-1)*9+3
       read(10,rec=irec) sst4
       write(20,'(i10,3f12.2)') ifcst,sst2(361,181),sst3(361,181),
     &     sst4(361,181)
      enddo
c
      write(20,*) 
     &'fcst days: member4(K): member5(K): member6(K):'
      write(20,*) 
     &'fcst days:   12z-run2:   12z-run3:   12z-run4:'
      do ifcst=1,45
       irec=(ifcst-1)*9+4
       read(10,rec=irec) sst2
       irec=(ifcst-1)*9+5
       read(10,rec=irec) sst3
       irec=(ifcst-1)*9+6
       read(10,rec=irec) sst4
       write(20,'(i10,3f12.2)') ifcst,sst2(361,181),sst3(361,181),
     &     sst4(361,181)
      enddo
c
      write(20,*) 
     &'fcst days: member7(K): member8(K): member9(K):'
      write(20,*) 
     &'fcst days:   18z-run2:   18z-run3:   18z-run4:'
      do ifcst=1,45
       irec=(ifcst-1)*9+7
       read(10,rec=irec) sst2
       irec=(ifcst-1)*9+8
       read(10,rec=irec) sst3
       irec=(ifcst-1)*9+9
       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