两层循环？nc提取赤道平流层纬向风平均值无法作图

qiaogj

想从nc提取赤道平流层纬向风平均值，做了两个嵌套，写完gs,写完ctl文件打卡画图，结果只有100hPa高度能出图，就是写的第一层有图
其它高度上画不出。求大神解答！附件gs和ctl


'reinit'
'sdfopen f:/qb/era/uwnd075.nc'
'set gxout fwrite'
'set fwrite  f:/qb/era/uwndave075.dat'
tt=1
zz=27
while (tt<=408)
'set lon 0'   
'set lat 0'
'set t ' tt
zz=27
while(zz<=37)
'set z 'zz  
'define ub=ave(u,lon=0,lon=359.25)'
'd ub'
say result
zz=zz+1
endwhile v
tt=tt+1
endwhile
'disable fwrite'
;

ctl文件如下
dset ^uwndave075.dat
title equator uwnd ave
undef  99999
xdef 1 linear 0 0.75
ydef 1  linear 0 0.75
zdef 11 levels  100 70 50 30 20 10 7 5 3 2 1
tdef 408 linear 00Z01SEP1979 1mo
vars 1
u 1 99  uwnd
endvars
;
如下图，画出100hPa层的风场正常，画70hPa就出问题了

做个霸气的木头

既然输出了11层，那么下图u后面不应该是11层吗，对不对

qiaogj

做个霸气的木头 发表于 2013-5-2 18:21

                               
登录/注册后可看大图

既然输出了11层，那么下图u后面不应该是11层吗，对不对

没错！就是这个错啊！找了几天不知道哪里错了……万分感谢

qiaogj

附图

loftydongshi

ctl里面经常犯这些小错误，哈哈，大家都仔细点~

loftydongshi

say result是在哪儿输出变量的结果么？

qiaogj

loftydongshi 发表于 2013-5-2 21:15

                               
登录/注册后可看大图

say result是在哪儿输出变量的结果么？

say result 可以在命令界面看到每个数据的输出结果，如果写的数据超出原来的维数或者文件出错会有显示

loftydongshi

qiaogj 发表于 2013-5-4 11:55

                               
登录/注册后可看大图

say result 可以在命令界面看到每个数据的输出结果，如果写的数据超出原来的维数或者文件出错会有显示

领教！！thx

赫小猫

做个霸气的木头 发表于 2013-5-2 18:21
既然输出了11层，那么下图u后面不应该是11层吗，对不对

一语点醒梦中人