cccccccccccccccccccccccccccccccccccccccccccc
c CFSv2 hindcast seasonal forecast SST Skill 
cccccccccccccccccccccccccccccccccccccccccccc
      parameter(im=360,jm=181)
      real*4 corl(im,jm)
      real*4 rmse(im,jm)
c
c corl==> correlatio, rmse==> root mean sqare error
c hindcasts from 1982-2009
c 7 leads, 12 initial months
c based on the 1st month of the target season
c 
c init=1,lead=1==>JFM fcst from DEC ic(lead=0-mon)
c init=1,lead=2==>FMA fcst from DEC ic(lead=1-mon)
c init=1,lead=3==>MAM fcst from DEC ic(lead=2-mon)
c init=1,lead=4==>AMJ fcst from DEC ic(lead=3-mon)
c init=1,lead=5==>MJJ fcst from DEC ic(lead=4-mon)
c init=1,lead=6==>JJA fcst from DEC ic(lead=5-mon)
c init=1,lead=7==>JAS fcst from DEC ic(lead=6-mon)
c 
c init=7,lead=1==>JAS fcst from JUN ic(lead=0-mon)
c init=7,lead=2==>ASO fcst from JUN ic(lead=1-mon)
c init=7,lead=3==>SON fcst from JUN ic(lead=2-mon)
c init=7,lead=4==>OND fcst from JUN ic(lead=3-mon)
c init=7,lead=5==>NDJ fcst from JUN ic(lead=4-mon)
c init=7,lead=6==>DJF fcst from JUN ic(lead=5-mon)
c init=7,lead=7==>JFM fcst from JUN ic(lead=6-mon)
c
c
      open(11,
     &file='corlRmseSSNxySST19822009.gr',
     &form='unformatted',
     &convert='big_endian',
     &access='direct',recl=im*jm*4)
c
      open(20,
     &file='corlRmseSSNxySST19822009.txt.sample',
     &form='formatted',access='sequential')
c
      write(20,*) '----------------------------------' 
      write(20,*) 'correlation/rmse at point longitude=180E,latitude=0'

      init=1        !initial time 1-->12 

      lead=1        !lead time 1-->7
      irec=(init-1)*2*7+lead
      read(11,rec=irec) corl
      irec=(init-1)*2*7+7+lead
      read(11,rec=irec) rmse
      write(20,*) 'JFM forecast from DEC initialCondition'
      write(20,'(2f8.4)') corl(181,91),rmse(181,91)

      lead=2        !lead time 1-->7
      irec=(init-1)*2*7+lead
      read(11,rec=irec) corl
      irec=(init-1)*2*7+7+lead
      read(11,rec=irec) rmse
      write(20,*) 'FMA forecast from DEC initialCondition'
      write(20,'(2f8.4)') corl(181,91),rmse(181,91)

      lead=3        !lead time 1-->7
      irec=(init-1)*2*7+lead
      read(11,rec=irec) corl
      irec=(init-1)*2*7+7+lead
      read(11,rec=irec) rmse
      write(20,*) 'MAM forecast from DEC initialCondition'
      write(20,'(2f8.4)') corl(181,91),rmse(181,91)

      lead=4        !lead time 1-->7
      irec=(init-1)*2*7+lead
      read(11,rec=irec) corl
      irec=(init-1)*2*7+7+lead
      read(11,rec=irec) rmse
      write(20,*) 'AMJ forecast from DEC initialCondition'
      write(20,'(2f8.4)') corl(181,91),rmse(181,91)
c
      stop
      end
cccccccccccccccccccccccccccccccc