program readspi

c spi3:
c FMA 1981 - JFM 2010
C from CFS 0-month lead fcst:
C Initial conditions: Jan 1981 - December 2009
C 15 members

      PARAMETER ( NI = 192, NJ = 94)
      PARAMETER ( nyear=29)

      real spi(ni,nj,12*nyear)

      open(11,file='SPI3_cfs_0mleadfcst_15runs.gr',
     1 form='unformatted',access='direct',recl=ni*nj*12*nyear*4)

       DO id=1,15

c------------------------------------------------
c k=1: forecast made in 1981.01 for FMA 1981
c k=348: forecast made in 2009.12 for JFM 2010
c------------------------------------------------

       read(11,rec=id) (((spi(i,j,k),
     1              i=1,ni),j=1,nj),k=1,12*nyear)

       end do

      stop
      end