求助用grads画湿Q矢量散度

心动勿行动

大家好，请问有人画过湿Q矢量吗？本人grads新手，想用湿Q矢量诊断暴雨。但是不知道怎么做，下面的程序是从一位同事那里获得的Q矢量散度的一个粗略的计算方法，但是500hpa画出来的都是0，700hpa和850hpa都是没有值，无法画，请众位大神指教，如果有低级错误，还请轻拍，楼主妹纸，打击太大就直接上天台去了。。。
'reinit'
'open d:\macc\NCEP\fnl_20070918_12_00.ctl'
'enable print d:\macc\NCEP\gsprint\070918_q_850hpa.gmf'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lev 850'
'f=2*7.292*sin(lat*3.14159/180)'
*'du=cdiff(ugrdprs*cos(lat*3.14159/180),y)'
*'dv=cdiff(vgrdprs,x)'
'du=ugrdprs(z-1)-ugrdprs(z+1)'
'dv=vgrdprs(z-1)-vgrdprs(z+1)'
'dp=100*(lev(z-1)-lev(z+1))'
'dx=cdiff(lon,x)*3.14159/180'
'dy=cdiff(lat,y)*3.14159/180'
'aa=du/dx'
'aa1=dv/dp'
'bb=dv/dx'
'bb1=du/dp'
'qx=f*(aa*aa1-bb*bb1)'
'aa2=du/dy'
'aa3=dv/dp'
'bb2=dv/dy'
'bb3=du/dp'
'qy=f*(aa2*aa3-bb2*bb3)'
'define qdiv=hdivg(qx,qy)*1e6'
'set gxout contour'
'set lon 115 125'
'set lat 22 36'
'run cbarn.gs'
'set ccolor 1'
'set cint 2'
'set clskip 1'
'set clopts -1 -1 0.11'
'print'
'disable print'

超台十八

你就跟着每一行都画一下吧。。。

d du, d dv。。。看哪里出问题

youyou

你去论坛搜索下相关帖子，你的这个GS看着有些乱。。别的不说，你要画的是什么的语句都没有。。。。。你去这个帖子看看是不是你想要的。。。http://bbs.06climate.com/forum.php?mod=viewthread&tid=32743

心动勿行动

youyou 发表于 2015-2-13 16:34
你去论坛搜索下相关帖子，你的这个GS看着有些乱。。别的不说，你要画的是什么的语句都没有。。。。。你去这 ...

谢谢提醒啊 我也才发现 怎么我贴的这个文件木有画图语句……

同路人..

给晨晨同学支持一下。。

心动勿行动

同路人.. 发表于 2015-2-13 19:43
给晨晨同学支持一下。。

难道不是应该帮我解决下吗……

同路人..

这个不是你的最终版吧。。 d 都没看到哦。。

str@nger

美女问题解决了吗？我是一名基层工作者，写论文想用Q矢量的东西，如果你搞懂了，能否将程序传给我一份，万谢！

yuhuiqing6

你好，我现在毕业论文也要分析湿Q矢量。请问你有完整的吗现在，求指导

d1268y

{:smile:}{:smile:}{:smile:}