求气象大神助攻

surelee

这是fortran程序：
implicit none
   integer,parameter::nx=37,ny=17,nz=1,nt=48
   dimension u(nx,ny,nz,nt),v(nx,ny,nz,nt)
   real u,v
   integer it,i,j,iz
   open(1,file='d:\ll\u850.dat')
   open(2,file='d:\ll\v850.dat')
   do it=1,nt
      read(1,*) ((u(i,j,1,it),i=1,nx),j=1,ny)
          read(2,*) ((v(i,j,1,it),i=1,nx),j=1,ny)
   enddo
   open(12,file='d:\ll\uv.grd',form='binary')
   do it=1,nt
       do iz=1,nz
               write(12) ((u(i,j,iz,it),i=1,nx),j=1,ny)
           enddo
           do iz=1,nz
               write(12) ((v(i,j,iz,it),i=1,nx),j=1,ny)
           enddo
   enddo


这是.ctl文件
dset d:\ll\uv.grd
title NCEP/NCAR REANAL YSIS
undef -9.99E+33
xdef 37 linear 60.0 2.5
ydef 17 linear 0.0 2.5
zdef 1 levels 850
tdef 48 linear JAN1982 1mo
vars 2
u 1 99 u wind(m/s)
v 1 99 v wind(m/s)
endvars


为什么在grads中输入
open d:\data\uv.ctl         
set t 1                             
set z 1                           
d u;v   

以上命令画不出风场图呢