Q矢量散度---修改时间或者高度出的图都是一样的T_T

Diary

850  15日20时

Q850 14081520

Q850 14081620

850 16日20时

在论坛上找到的程序  可是图真的都一样

Reinit

clear

sdfopen e:\nc\g\hgt.2014.nc

define hgtc=hgt.1

define ug= -9.8*cdiff(hgtc,y)/(111177*cdiff(lat,y))

define vg= 9.8*cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))

define dugdx=cdiff(ug,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))

define dvgdx=cdiff(vg,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))

define dtdx=cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))

define dugdy= cdiff(ug,y)/(111177*cdiff(lat,y))

define dvgdy= cdiff(vg,y)/(111177*cdiff(lat,y))

define dtdy= cdiff(hgtc,y)/(111177*cdiff(lat,y))

define qx=(-287.05/lev)*(dugdx*dtdx+dvgdx*dtdy)*100000000000000

define qy=(-287.05/lev)*(dugdy*dtdx+dvgdy*dtdy)*100000000000000

define div=hdivg(qx,qy)*100000

set mpdset cnworld

set map 15 4 15

set grads off

set grid off

set lev 850

set lon 90 120

set lat 30 45

set t 912

set gxout contour

set ccolor 2

set cthick 10

set clab forced

set clskip 1

d div

set mpdset cnriver

set map  3  1 8

draw map

printim f:\Q850hpa081620.png white x1200 y800

请大家帮忙~

clare

是不是define的那些需要放在set t 和 set lev的后面，你试试看？

river

计算之前加上高度和时间的变化设置，比如set lev 1000 10   set t 1 250
试试看

Diary

clare 发表于 2015-4-18 10:39
是不是define的那些需要放在set t 和 set lev的后面，你试试看？

difine 的部分是把Q矢量计算出来了吧？

clare

Diary 发表于 2015-4-18 14:20
difine 的部分是把Q矢量计算出来了吧？

楼上的解决方法你尝试了吗？我猜测可能原因也是你没有先set 范围，而grads的初始画图是默认set t 1等。所以楼上的方法应该正解。试试看啊

Diary

clare 发表于 2015-4-19 09:06
楼上的解决方法你尝试了吗？我猜测可能原因也是你没有先set 范围，而grads的初始画图是默认set t 1等。所 ...

嗯 修改后 出图了 但是 我看图反应出的情况和真实的暴雨情况不同 不知道怎么回事

Diary

clare 发表于 2015-4-19 09:06
楼上的解决方法你尝试了吗？我猜测可能原因也是你没有先set 范围，而grads的初始画图是默认set t 1等。所 ...

如果我要做 非地转湿Q矢量 应该就不是用这个了吧？

还有一个问题 111177是代表什么？

Diary

river 发表于 2015-4-18 11:53
计算之前加上高度和时间的变化设置，比如set lev 1000 10   set t 1 250
试试看

谢谢~已解决~

clare

Diary 发表于 2015-4-19 11:21
如果我要做 非地转湿Q矢量 应该就不是用这个了吧？

还有一个问题 111177是代表什么？

这个貌似是准地砖的吧？我没有计算过，所以也没有查阅，不好意思，帮不上什么

Diary

clare 发表于 2015-4-19 19:45
这个貌似是准地砖的吧？我没有计算过，所以也没有查阅，不好意思，帮不上什么

那还是感谢你了-----第一个在论坛上给我多次热心回复的~