nc转格点多层数据时出错

我就是-猴子 · 发表于 2014-12-31 15:38:07

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

'reinit'
'sdfopen  e:/renyj/grads/nc2grd/shum.mon.mean.nc'
'set  gxout  fwrite'
'set  fwrite  e:/renyj/grads/nc2grd/shum10.jja.grd'

'set x 1 144'
'set y 1 73'

i=162
**162对应的是1961年6月*****
while(i<=786)
**786对应的是2013run年6月*****
'set t ' i

iz=1
while(iz<=8)
'set z ' iz
'define hgt1= ave(shum,t='i',t='i+2')'
'd hgt1'
iz=iz+1
endwhile

i=i+12
endwhile
'disable fwrite'
'reinit'

这是我的脚本，但画图显示每一层的结果都一样
都显示的是第一层的结果
麻烦高手给指导一下

我就是-猴子 · 发表于 2014-12-31 15:58:29

dset e:\renyj\grads\shum\shum.jja.grd
title grid data
undef -9.99e8
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 53 linear 01jan1961 1yr
vars 1
shum=>shum 0 t,y,x monthly means of hgt
endvars

这是ctl文件

鲸鱼便子 · 发表于 2014-12-31 16:27:10

ctl里面变量后面不应该是0吧？
写成 shum 8 99 ******试试

river · 发表于 2014-12-31 17:38:01

你确定是gs脚本有问题，而不是你的ctl 吗？

youyou · 发表于 2014-12-31 21:33:56

我也觉得的是CTL文件里的那个变量后面的0有问题。变量定义只有一层，你又如何要求人家出多层的图？

我就是-猴子 · 发表于 2015-1-4 16:59:54

谢谢各位的解答确实是ctl文件出错了

		自动登录	找回密码
密码			立即注册

[脚本编辑] nc转格点多层数据时出错

登录后查看更多精彩内容~