求助！！第一时间层可以画出图为什么之后的就画不出来了

lavenderSC361

求助各位大神，用GRADS转化nc资料为grd：uwnd.2102.nc和vwnd.2012.nc
程序如下：
'reinit'
'sdfopen g:\lunwen\uwnd.2012.nc'
'sdfopen g:\lunwen\vwnd.2012.nc'
'set fwrite g:\lunwen\u.grd'
'set gxout fwrite'
'set fwrite g:\lunwen\v.grd'
'set gxout fwrite'
'set lon 80 220'
'set lat -10 80'
'set z 3'
it=210
while(it<=217)
'set t 'it''
'd uwnd'
it=it+1
'd vwnd'
endwhile
'disable fwrite'
;
写是这么写的实际是分开来生成grd文件的，然后用FORTRAN编写程序计算涡度：
program main
implicit none
integer i,j,nx,ny,it,n
parameter(nx=57,ny=37,it=8)
real u(nx,ny,it),v(nx,ny,it),lat(ny)
real::a=6371000.0,pi=3.14
real vox1(nx,ny,it),vox2(nx,ny,it),vox(nx,ny,it)
real vox3(nx,ny,it),vox4(nx,ny,it),vox5(nx,ny,it)
open(1,file='g:\lunwen\vox.grd',form='binary')
open(2,file='g:\lunwen\u.grd',form='binary')
open(3,file='g:\lunwen\v.grd',form='binary')
do n=1,it
do j=1,ny
read(2)(u(i,j,n),i=1,nx)
read(3)(v(i,j,n),i=1,nx)
end do
end do
close(2);close(3)
do j=1,ny
lat(j)=10.0+float(j-1)*1.0
enddo
do n=1,it
do j=2,ny-1
  do i=2,nx-1
  vox1(i,j,n)=(v(i+1,j,n)-v(i-1,j,n))/((pi/180.0)*cosd(lat(j)))
  vox2(i,j,n)=(u(i+1,j,n)-u(i-1,j,n))/(pi/180.0)
  vox3(i,j,n)=(1.0/(2.0*a))*(vox1(i,j,n)-vox2(i,j,n)+2.0*u(i,j,n)*tand(lat(j)))
  vox4(i,j,n)=(-1.0)*(u(i,j,n)*(vox3(i+1,j,n)-vox3(i-1,j,n))/(2.0*a*(pi/180.0)*cosd(lat(j))))
  vox5(i,j,n)=(-1.0)*(v(i,j,n)*(vox3(i+1,j,n)-vox3(i-1,j,n))/(2.0*a*(pi/180.0)))
  vox(i,j,n)=(vox4(i,j,n)+vox5(i,j,n))*(10**10)
enddo
enddo
enddo
write(1)vox
      end
得出了vox.grd然后写ctl画图：
dset g:\lunwen\vox.grd
undef -9.99E+33
title NCEP/NCAR REANALYSIS PROJECT
xdef   57 linear   80.000  2.500
ydef   37 linear   -10.000  2.500
zdef   1 levels 850
tdef   8 linear   00:00z28JUL2012    24hr
vars  1
vox 11 99 temprature
endvars

然后画图就只能画出t=1的图。。。。。各位我是哪里出错了呢？帮个忙呀！

lavenderSC361

我知道哪里错了。。。。

nuister

觉悟挺快，刚发完帖子就知道错哪了，哈哈

做个霸气的木头

这个反映速度就像，上学去找老师问题，到了办公室门口就会做了。

lavenderSC361

nuister 发表于 2013-5-7 19:31

                               
登录/注册后可看大图

觉悟挺快，刚发完帖子就知道错哪了，哈哈

哈哈，犯了低级错误呢~

lavenderSC361

做个霸气的木头 发表于 2013-5-7 23:16

                               
登录/注册后可看大图

这个反映速度就像，上学去找老师问题，到了办公室门口就会做了。

被嘲笑了= =

river

lavenderSC361 发表于 2013-5-8 12:38

                               
登录/注册后可看大图

被嘲笑了= =

让大家笑一下，娱乐一下没什么不好的啊。不过不是嘲笑哦

做个霸气的木头

lavenderSC361 发表于 2013-5-8 12:38

                               
登录/注册后可看大图

被嘲笑了= =

额，真心不是嘲笑，都会遇到这样的问题，感慨一下下啦~

愉悦

楼主反应速度好快 我都没找出错误在哪

tcj-mlh

觉悟挺快，刚发完帖子就知道错哪了，哈哈