水汽通量散度计算gs文件

海贝

帮看下水汽通量散度的gs   不知道错哪了
'reinit'
'sdfopen d:/necp/uwnd/uwnd.2011.nc’
'sdfopen d:/necpwndwnd.2011.nc'
'sdfopen d:/necp/rhum/rhum.2011.nc'
'sdfopen d:/necp/air/air.2011.nc'
'set lon 105 120'
'set lat 20 35'
'set lev 850'
'set t 660'
'set grads off'
'set grid off'
'define tt=air.4+273.16'
'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=rhum.3/100.*qs'
'define xflux=qq*uwnd.1'
'define yflux=qq*vwnd.2'
'define sqdiv=hdivg(xflux,yflux)*1e8'
'd sqdiv'

传说中的谁

'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
似乎应该写成：'define es=6.11*exp(17.26*(tt-273.16))/(tt-35.86)'

KaKa

会不会是 'sdfopen d:/necpwndwnd.2011.nc' 文件名和路径有误，导致后面引用V风分量计算时vwnd.2无法识别？

飞雪残剑

'define es=6.11*exp((17.26*(tt-273.16))/(tt-35.86))'

这样就对了，但是出的图不能到图边，不知道怎么回事？

飞雪残剑

图形如下：

传说中的谁

飞雪残剑 发表于 2011-8-1 22:44

                               
登录/注册后可看大图

图形如下：

很简单，你定义变量的时候把范围设大一点，画图的时候再设置成你要的范围。

feiji158

飞雪残剑 发表于 2011-8-1 22:44

                               
登录/注册后可看大图

图形如下：

请问市县级底图 怎么画的？

feiji158

通量散度是X出来，的70度的 要68度资料才能插出来的吧

飞雪残剑

传说中的谁 发表于 2011-8-1 22:47

                               
登录/注册后可看大图

很简单，你定义变量的时候把范围设大一点，画图的时候再设置成你要的范围。

好的，谢谢传说！

飞雪残剑

feiji158 发表于 2011-8-1 23:15

                               
登录/注册后可看大图

通量散度是X出来，的70度的 要68度资料才能插出来的吧

有地图文件，你放到相应目录下，在grads安装软件下面的DAT目录下