求水平散度时遇到的奇葩问题

小傅之追梦

'reinit'
'sdfopen F:\bl\div\uwnd.2010.nc'
'sdfopen F:\bl\div\vwnd.2010.nc'
'set fwrite F:\bl\div\div.grd'
'set gxout fwrite'
'define div=hdivg(uwnd.1,vwnd.2)*1e5'
'd div'
'disable fwrite'
;

在grads里面直接输入下面的程序，出图是正常的，但是编成GS运行之后，配上CTL之后出图就错误了。两段程序是完全一样的，不知道这是神马原因，求大神帮忙看看啊~~~~~

正确的图

错误的图

小肥虾619

经纬度的问题吧？提取一个区域试试

lqouc

不管你cl写的对不对，输出文件的售后最好把经纬度和时间都定义一下，如果输出多时次，时间是要用循环的。好好改改你的gs吧。

小傅之追梦

小肥虾619 发表于 2013-4-28 12:02

                               
登录/注册后可看大图

经纬度的问题吧？提取一个区域试试

还是不行啊，在所提取的边界值异常大

小傅之追梦

lqouc 发表于 2013-4-28 12:07

                               
登录/注册后可看大图

不管你cl写的对不对，输出文件的售后最好把经纬度和时间都定义一下，如果输出多时次，时间是要用循环的。好 ...

这是我按照你的建议修改的GS和CTL文件，
'reinit'
'sdfopen F:\bl\div\uwnd.2010.nc'
'sdfopen F:\bl\div\vwnd.2010.nc'
'set fwrite F:\bl\div\div.grd'
'set gxout fwrite'
'set x 41 81'
'set y 40 60'
*'set lon 100 200'
*'set lat 7.5 57.5'
'set t 1'
'define div=hdivg(uwnd.1,vwnd.2)*1e5'
'd div'
'disable fwrite'
;
##############################
dset F:\bl\div\div.grd
title 4x daily NMC reanalysis(2010)
undef -9.99e+33
xdef 41 linear 100 2.5
ydef 21 linear 7.5 2.5
zdef 1 levels 1000
tdef 1 linear 00Z21JUL2010 360mn
vars 1
div 0 99
endvars
;可还是在图形的边界显示异常的大。。。

lqouc

小傅之追梦 发表于 2013-4-28 12:27

                               
登录/注册后可看大图

这是我按照你的建议修改的GS和CTL文件，
'reinit'
'sdfopen F:\bl\div\uwnd.2010.nc'

你的缺省值确定是这个么？反正我的grads如果计算之后在输出缺省值会变成-9.9900000e8。
你可以看看你的缺省值是什么。

小傅之追梦

lqouc 发表于 2013-4-28 12:32

                               
登录/注册后可看大图

你的缺省值确定是这个么？反正我的grads如果计算之后在输出缺省值会变成-9.9900000e8。
你可以看看你的缺 ...

可以怎么查看缺测啊？。。。。我的缺测值设为你给的这个后图形就正常多了。。。

小傅之追梦

lqouc 发表于 2013-4-28 12:32

                               
登录/注册后可看大图

你的缺省值确定是这个么？反正我的grads如果计算之后在输出缺省值会变成-9.9900000e8。
你可以看看你的缺 ...

真是神级的人啊，谢谢谢谢谢谢~~~~~我都纠结一个上午了这个问题。。。。现在就是不知道怎么确定这个缺测值了

mofangbao

小傅之追梦 发表于 2013-4-28 12:42

                               
登录/注册后可看大图

真是神级的人啊，谢谢谢谢谢谢~~~~~我都纠结一个上午了这个问题。。。。现在就是不知道怎么确定这个缺测值 ...

基础知识掌握不够 请自行补充

小傅之追梦

mofangbao 发表于 2013-4-28 14:22

                               
登录/注册后可看大图

基础知识掌握不够 请自行补充

清风大侠教训的是啊，好好学习，天天向上~~~~