单站水汽通量散度的问题

xiahou90

我想求取某个点的水汽通量散度随时间和高度的变化，但是散度是要求偏导数的，也就是说经纬度要在一定范围内变化，所以我用grads里的divg函数总是求不出，求高人指导
'open g:/date/d02_2012-06-04_15-00'
'define q1=qvapor*u/9.8'
'define q2=qvapor*v/9.8'
'define qfdiv=hdivg(q1,q2)*1e8'
'set lev 1000 100'
'set t 1 18'
'set lat 42'
'set lon 86'
'set gxout contour'
'set cstyle 1'
'set cmin 0'
'set black 0 0'
'd qfdiv'
'set cstyle 5'
'set cmax 0'
'set black 0 0'
'd qfdiv'
'printim g:/date/pic/zhandian/qfdiv.gif
white'

传说中的谁

根据散度的公式和定义可以知道，一个点当然没法求散度啊
如果你非要做的话，先插值到格点再求吧

river

你画等值线还设置两次线形，还设置了'set cmin 0'和'set cmax 0'，您明白什么意思不？
一个是表示不画大于0的，一个是不画小于0的，那您不就是让它别画么，所以它随您愿了图上什么都没有。
还有'set black x y '是说要给小于x和大于y的地方加上阴影，您直接放在等值线下面干嘛，这个命令不适合画等值线啊

xiahou90

river 发表于 2013-5-8 21:35

                               
登录/注册后可看大图

你画等值线还设置两次线形，还设置了'set cmin 0'和'set cmax 0'，您明白什么意思不？
一个是表示不画大于 ...

我的意思是0以上的值画为实线，0以下的值画为虚线 0线不画 你指出的这几个命令应该没错 因为我画其他物理量都没问题 只是画水汽通量散度的时候遇到问题了。

xiahou90

传说中的谁 发表于 2013-5-8 21:23

                               
登录/注册后可看大图

根据散度的公式和定义可以知道，一个点当然没法求散度啊
如果你非要做的话，先插值到格点再求吧

我的数据是格点数据 求指导哇 真的很想弄出来  因为论文要用到的缘故

river

xiahou90 发表于 2013-5-9 13:45

                               
登录/注册后可看大图

我的意思是0以上的值画为实线，0以下的值画为虚线 0线不画 你指出的这几个命令应该没错 因为我画其他物理 ...

我觉得吧你很幽默啊，grads本身默认就是正值画实线，负值画虚线，根本就不需要你自己设置

xiahou90

river 发表于 2013-5-9 22:22

                               
登录/注册后可看大图

我觉得吧你很幽默啊，grads本身默认就是正值画实线，负值画虚线，根本就不需要你自己设置

先不管实线虚线的问题了  水汽通量散度能画出来不 求指导

river

xiahou90 发表于 2013-5-9 22:54

                               
登录/注册后可看大图

先不管实线虚线的问题了  水汽通量散度能画出来不 求指导

我说过了，先把你那几个没用的命令set cmin、 set cmax 、set black都拿掉看能不能出图，能出的话你再逐句加上那几句命令来验证你说的“这几个命令应该没错”

xiahou90

'define qfdiv=hdivg(q1,q2)*1e8'
'set lev 1000 100'
'set t 1 18'
'set lat 42'
'set lon 86'
'set gxout contour'
*'set cstyle 1'
*'set cmin 0'
*'set black 0 0'
'd qfdiv'
*'set cstyle 5'
*'set cmax 0'
*'set black 0 0'
*'d qfdiv'

这样画出来 grads显示一个定值1e30  
加上我*掉的那些命令后 两次d qfdiv后 两次都是显示定值1e30

永远的王者

计算水汽通量散度的时候，设置经纬度在一个范围内变化；显示的时候再固定经纬度。