求助，grads对于格点数据的要求，如何修改FORTRAN文件

水木缘

FORTAN程序

!!var()各点数据,varcli()月平均气候场!!!
      program main
      parameter(ii=37,jj=17,mo=12,yr=4)
      real var(ii,jj,mo,yr)
        real varcli(ii,jj,mo)
        real iy,m,i,j,ny,irec
      open(7,file='F:\QiXianTonJi\h500.dat')
ccccccccccccccccccc 读数据资料 h500  ccccccccccccccccccccccccccc
       do  iy=1,yr
       do  m=1,mo
       read(7,1000) !读 1982 1等类似的两个整形数据
       read(7,3000) ((var(i,j,m,iy),i=1,ii),j=1,jj)
       enddo
       enddo
1000   format(2i7)
3000   format(37f8.1)
       close(7)
cccccccccccccccccccc以下自己编程计算cccccccccccccccccccccccccccc  
!!!!!!!!!!!!!!!!!!!!!!!!气候场!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        do m=1,mo
      do j=1,jj
        do i=1,ii
          temp=0.0
          ny=0.0
        do iy=1,yr
           temp=temp+var(i,j,m,iy)
           ny=ny+1.0
          end do
          varcli(i,j,m)=temp/ny
        end do
        end do
        end do
             open(12,file='F:\QiXianTonJi\varcli\varcli.grd',
     & form='unformatted',access='direct',recl=ii*jj*4 )
              irec=0
         do it=1,mo
              irec=irec+1
            write(12,rec=irec)((varcli(i,j,it),i=1,ii),j=1,jj)
         enddo
        close(12)
        end

CTL程序

dset  F:\QiXianTonJi\varcli\varcli.grd
undef -9.99E+33
title NCEP/NCAR REANALYSIS PROJECT
xdef   37 linear   60.000  2.500
ydef   17 linear   0.000  2.500
zdef    1 levels   500
tdef 12 linear   JAN1982   1mo
vars  1
H  1 99 H500
endvars

gs文件
'reinit'
'open F:\QiXianTonJi\varcli\varcli.ctl'
'set lat 0 40'
'set lon 60 150'
'set lev 500'
mo=1
while(mo<13)
i=mo
'set t 'i''
'set cint 15'
'd H'
'set annot 1 3'
'draw title  'i'month'
'printim F:\QiXianTonJi\varcli\picture\varcli'mo'.png white'
'c'
mo=mo+1
endwhile
;