Q矢量的计算

泡沫学天气

这是gs
'reinit'
'open F:\newwork\shiQ\wrfout_d02_2017-03-24.ctl'
'enable print F:\newwork\shiQ\shiQ.gmf'
'define hgtc=hgt.1'
*'define u=uwnd.2'
*'define v=vwnd.3'
'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 112 120'
'set lat 22 33'
'set t 3'
'set gxout contour'
'set ccolor 2'
'set cthick 2'
'set clab forced'
'set clskip 1'
'set cint -1e+09'
'd div'
'set gxout barb'
'set ccolor 1'
'd skip(u,16,16);skip(v,16,16)'
*'d u;v'
'set mpdset cnriver'
'set map  3  1 8'
'draw map'
'draw title 26MAY2014_700hPa_Q'
'printim  F:\newwork\shiQ\shiQ.png white x1200 y800'
'reinit'


gs是在论坛里找的
做出来的图片是这样的 ，不知道哪里错了，求大家帮忙

atsmd

盼望有大侠指导解惑，学习了

泡沫学天气

atsmd 发表于 2017-6-5 23:28
盼望有大侠指导解惑，学习了

{:eb303:}{:eb303:}

weiyh

楼主搞定了么？

大笨熊

楼主搞定了吗？同问

泡沫学天气

已经解决了，由于资料的精度太高了，所以平滑掉就好了

大笨熊

多谢楼主分享，学习了

slc2005011020

感谢楼主！！！！！！！！！！！！！！！！！！！！！！！！

小预报员

{:eb328:}{:eb334:}{:eb334:}{:eb334:}

tangjie199365

泡沫学天气 发表于 2017-7-5 12:22
已经解决了，由于资料的精度太高了，所以平滑掉就好了

具体怎么做 楼主 求解求解