program reaa500yr

      PARAMETER ( NI = 144, NJ = 73, ny=500)
      PARAMETER ( mi = 180, mj = 139)
      PARAMETER ( ix = 192, il = 94)
      real a(ni,nj),b(mi,mj),c(ix,il)
      real wind(ni,nj,12,ny)
      real sst(mi,mj,12,ny)
      real prate(ix,il,12,ny)

c-----------------------------------------------------------
c All monthly mean data from NORTH Pole to SOUTH Pole in lat
c (big endian)

C For Wind (m/s)
c  grid: 144x73, 2.5x2.5 (lon x lat)
c  lat  90.000000 to -90.000000 by 2.500000
c  lon 0.000000 to 357.5 by 2.500000

c For SST (K)
c  grid: 180x139,  2x1 (lon x lat)
c  lat  64.0 to -74.0 by 1.0
c  lon  1.0 to 359.0 by 2.0
c  Value on LAND  9.999E+20

c For prate (kg/m^2/s, Precipitation rate)
c  grid: 192x94,  T62
c  lon  0.0 to 358.125 by 1.875
c
c  lat (T62, from SOUTH to NORTH):
c -88.542 -86.653 -84.753 -82.851 -80.947 -79.043 -77.139 -75.235 -73.331 -71.426
c -69.522 -67.617 -65.713 -63.808 -61.903 -59.999 -58.094 -56.189 -54.285 -52.380
c -50.475 -48.571 -46.666 -44.761 -42.856 -40.952 -39.047 -37.142 -35.238 -33.333
c -31.428 -29.523 -27.619 -25.714 -23.809 -21.904 -20.000 -18.095 -16.190 -14.286
c -12.381 -10.476  -8.571  -6.667  -4.762  -2.857  -0.952   0.952   2.857   4.762
c   6.667   8.571  10.476  12.381  14.286  16.190  18.095  20.000  21.904  23.809
c  25.714  27.619  29.523  31.428  33.333  35.238  37.142  39.047  40.952  42.856
c  44.761  46.666  48.571  50.475  52.380  54.285  56.189  58.094  59.999  61.903
c  63.808  65.713  67.617  69.522  71.426  73.331  75.235  77.139  79.043  80.947
c  82.851  84.753  86.653  88.542

c-----------------------------------------------------------
      
      open(11,file='u200_2.5x2.5_freerun500yr',
     1 form='unformatted',status='old')

      open(12,file='sst_2x1_freerun500yr',
     1 form='unformatted',status='old')

      open(13,file='prate_T62_freerun500yr',
     1 form='unformatted',status='old')

      do iyear=1,ny
      do imonth=1,12

      read(11) a

        do i=1,ni
        do j=1,nj
         wind(i,j,imonth,iyear)=a(i,j)
        end do
        end do

      read(12) b

        do i=1,mi
        do j=1,mj
         sst(i,j,imonth,iyear)=b(i,j)
        end do
        end do

      read(13) c

        do i=1,ix
        do j=1,il
         prate(i,j,imonth,iyear)=c(i,j)
        end do
        end do


      end do
      end do

      stop
      end