同样是grads提取数据问题，这次数据是fnl，里面自己算个物理量，维度有误

siljechen

具体是这样，我想画个涡度、散度的斜剖面，看了斜剖面的标准模板gs，多少还是不熟悉，决定先把要算的物理量单独先算出来。回头再套斜剖的模板，但遇到一系列问题。

先上gs:
'reinit'
'open G:\snowstorm\fnl_20130217_00_00_c.ctl'
'set gxout fwrite'
'set fwrite g:\snowstorm\gs\hdivg.grd'
i=1
while(i<=12)
'set t 'i''
j=1
while(j<=20)
'set z 'j''
'set y 1 181'
'set x 1 360'
'define div=hdivg(ugrdprs,vgrdprs)'  *这样写有问题么？
'd div'
j=j+1
endwhile
i=i+1
endwhile
'disable fwrite'
;

ctl如下：
dset g:\snowstorm\gs\hdivg.grd
title hdivg
undef -9.99e+33
ydef 179 linear -89.0 1
xdef 358 linear 2.0 1.0
zdef 20 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150
tdef 12 linear 00Z17feb2013 6hr
vars 1
div 0 99 hdivg
endvars

考虑到中央差分后，两头少了一圈，格点我就缩了两个。

根本错误在：z dimension上好像没有数据出来。检查了半天，改了n遍gs后，连数据本身也有点奇怪了，大神帮看看

siljechen

运行gs没有报错，打开ctl时候，感觉有问题，因为lev set 1000 1000这显然是z方向不对啊，但我明明循环写进去了啊

siljechen

运行gs没有报错，打开ctl时候，感觉有问题，因为lev set 1000 1000这显然是z方向不对啊，但我明明循环写进去了啊

兰北

我也研究研究，学习下

river

看gs是没有问题，但是要说明的是计算完散度确实会少一圈，但他还是有数据的，只是变成缺测了，所以你的ctl不能随便就减少经纬度范围。看错误提示说的是在维数Z上没有数据。你所用的那个fnl资料里面有那么多时次吗？

liutaoz21

前来围观~~~~~~~

cris

看了一下，UGRDprs 26 33,100,0 ** u wind [m/s]，VGRDprs 26 34,100,0 ** v wind [m/s]都是有26层的数据。

siljechen

gs不变，把ctl改回来，向river说的那样不要减少经纬度范围。结果仍然不变。
dset g:\snowstorm\gs\hdivg.grd
title hdivg
undef -9.99e+33
xdef 360 linear 0.0 1.0
ydef 181 linear -90.0 1.0
zdef 20 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150
tdef 12 linear 00Z17feb2013 6hr
vars 1
div 0 99 hdivg
endvars


先声明fnl资料里面ugrdprs,vgrdprs就是nc里面的u,v，是有很多层的。


fnl自己配套的ctl里面有描述。
UGRDprs 26 33,100,0 ** u wind [m/s]
VGRDprs 26 34,100,0 ** v wind [m/s]

先运行gs,没有错误提示。
运行ctl,，感觉lev set 1000 1000，就是有问题
。
画一下试试，数据好像是有的。
但是同一时次，数据就只有一层，无论怎么设定set lev 850也好，500也好，都是一样的图





查看下是不是z方向没有数据
发现就是没有，{:soso_e136:}{:soso_e136:}，到底出了什么事情啊！！！





无法理解。哭啊。大神们再给点建议呢？

斥鷃

ctl里面设置是
div 0 99 hdivg，
这里的0应该改成20那样子的层次数，要不然应该会认为div变量只有地面层有数据吧~这是我的想法~

siljechen

斥鷃 发表于 2013-7-3 22:20

                               
登录/注册后可看大图

ctl里面设置是
div 0 99 hdivg，
这里的0应该改成20那样子的层次数，要不然应该会认为div变量只有地面层有 ...

试了下，还是同样的问题呢。