c    67  1         2         3         4         5         6         7 2
c* Writes out the GPCP monthly 2.5-d data and covar monthly EOFs for
c* the 384 months, 1979-2010.  The grid is centered on 1.25E-358.75E by
c* 88.75S-88.75N.  Data and EOFs are global.
c
      parameter(im=144,jm=72,ntm=384,nmo=30)
      real a2(im,jm),ts(nmo)
      character fm*11,st*3
      data fm/'unformatted'/, st/'old'/
c
c* Input files:
      open(51,file='data/ascii.gpcp2.5.dat')
      open(52,file='data/ascii.gpcp2.5.anom.dat')
      open(53,file='data/ascii.gpcp2.5.anom.eofev.dat')
      open(54,file='data/ascii.gpcp2.5.anom.eofts.dat')
c
      open(21,form=fm,status=st,file='data/gpcp2.5.1979.2010.dat')
      open(22,form=fm,status=st,file='data/gpcp2.5.anom.1979.2010.dat')
      open(23,form=fm,status=st,file='data/gpcp.mon.eof2.5.ev')
      open(24,form=fm,status=st,file='data/gpcp.mon.eof2.5.ts')
c
c* Get the data.
      do iu=21,22
       iuw=iu+30
       do nt=1,ntm
        read(iu) a2
        do i=1,im
         write(iuw,81) (a2(i,j),j=1,jm)
        end do
       end do
       close(iu)
       close(iuw)
      end do
c
c* Get the spatial EOF modes.
      do nm=1,nmo
       read(23) a2
       do i=1,im
        write(53,81) (a2(i,j),j=1,jm)
       end do
      end do
      close(23)
      close(53)
c
c* Get the EOF time series.
      do nt=1,ntm
       read(24) (ts(nm),nm=1,nmo)
       write(54,81) (ts(nm),nm=1,nmo)
      end do
      stop
 81   format(200f7.2)
      end