垂直水汽通量时间-纬度剖面图

...该...

各位达人好，我根据论坛里面的脚本文件画了个水汽通量的时间-剖面图，感觉很不对劲啊，跟人家做的相差太大了（见附件）。请各位高手给新手指正下。谢谢！本来是想贡献点的，一看自己也没多少啊，多体谅下。哈哈


脚本：
'reinit'
'sdfopen F:\data\201205-09\air.2012.nc'
'sdfopen F:\data\201205-09\uwnd.2012.nc'
'sdfopen F:\data\201205-09\vwnd.2012.nc'
'sdfopen F:\data\201205-09\rhum.2012.nc'
'sdfopen F:\data\201205-09\shum.2012.nc'
'set t 1 120'
'set lev 850'
'define prs=lev'
'define g=9.8'
'define es=6.112*exp(17.67*(air.1-273.16)/(air.1-29.65))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=qs*rhum.4/100'
'set lev 1000 300'
'define qdiv=hdivg(q*uwnd.2*1000/g,q*vwnd.3*1000/g)*100000'
'set lon 1'
'define qdiva=ave(qdiv,lon=100,lon=125)'
'define u1=ave(uwnd.2,lon=100,lon=125)'
'define v1=ave(vwnd.3,lon=100,lon=125)'

'set mpdset cnworld'
'set grads off'
'set lev 500'
'set xyrev on'
'set lat 16 29'
'set t 1 120'
'set gxout contour'
'set cint 0.4'
'd qdiva'
*'cbarn'
'd u1;v1'

'printim f:/mois.gif x2400 y1600 white'


附件

river

您这个gs仿照的是哪个帖子的啊，光set lev就有三个，每次都不一样，那你是想要话哪个高度的啊。显然楼主没有学到真正的方法，只是拿来用而已嘛

kongfeng0824

river 发表于 2013-3-19 07:24

                               
登录/注册后可看大图

您这个gs仿照的是哪个帖子的啊，光set lev就有三个，每次都不一样，那你是想要话哪个高度的啊。显然楼主没有 ...

看来图形的剖面图好像和看起来的不是很好。

river

kongfeng0824 发表于 2013-3-19 07:34

                               
登录/注册后可看大图

看来图形的剖面图好像和看起来的不是很好。

对啊，很多东西都是看起来容易，做起来就出各种问题

jacci_zou

感觉你的GS有点乱噶，发现有问题，坐标设置看不懂，'set t 1 120' 'set lev 850'，define hdivg之前， 你指固定了个层次，也就默认了是全球的经纬度范围，那么，lon lat t都是变量，计算出来肯定有问题，其次，是后面的'set lev 1000 300'这个是准备做垂直积分的么？前面在定义变量之前你已经固定好层次了，所以这里再来设置已经没有作用了， 'set lon 1' 也没怎么明白是做什么用的！

鳗鱼宝

学习了！

...该...

我 是从论坛里面找的，但是自己只是试了下，结果跟别人做的那图完全不一样了。。。该怎么改啊。。。5555，自己好好摸索下，再请教你们。。。

...该...

我把高度都作了调整，还是做不出那样的图啊。不知道该怎么改啊。望各位提点宝贵的建议
'reinit'
'sdfopen F:\data\201205-09\air.2012.nc'
'sdfopen F:\data\201205-09\uwnd.2012.nc'
'sdfopen F:\data\201205-09\vwnd.2012.nc'
'sdfopen F:\data\201205-09\rhum.2012.nc'
'sdfopen F:\data\201205-09\shum.2012.nc'
'set lev 850'
'define prs=lev'
'define g=9.8'
'define es=6.112*exp(17.67*(air.1-273.16)/(air.1-29.65))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=qs*rhum.4/100'
*'set lev 1000 300'
'define qdiv=hdivg(q*uwnd.2*1000/g,q*vwnd.3*1000/g)*100000'
'set lon 1'
'set t 1 120'
'define qdiva=ave(qdiv,lon=100,lon=125)'
'define u1=ave(uwnd.2,lon=100,lon=125)'
'define v1=ave(vwnd.3,lon=100,lon=125)'

'set mpdset cnworld'
'set grads off'
*'set lev 500'
'set xyrev on'
'set lat 16 29'
*'set t 1 120'
'set gxout shaded'
'set cint 0.4'
'd qdiva'
*'cbarn'
'd u1;v1'

'printim f:/mois.gif x2400 y1600 white'

jacci_zou

'reinit'
'open g:\zj\1118\1118.ctl'
'set mpdset cnworld fujian'
'set lon 109 126'
'set lat 14 31'
'set t 6 20'
'set lev 1000 100'
'set grads off'
'set grid off'
'define z=hgtprs'
'define t=tmpprs'
'define prs=1.333*760*pow((1-z/44300),5.256)'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define qs=0.62197*es/(prs-0.378*es)*1000'
'define q1=qs*ugrdprs/9.8'
'define q2=qs*vgrdprs/9.8'
'define qfdiv=hdivg(q1,q2)*1e5'
'set lon 118'
'set lat 24'
'set lev 1000 300'
'set cthick 5'
'set csmooth on'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set xlint 4'
'set ylint 100'
'set zlog on'
'set gxout shaded'
'set cmax 0'
'd qfdiv'
'd q1;q2'
'cbarn 0.7 1 10.62'
'draw title moisture flux and moisture flux divergence '
'set arrlab off'
'printim g:\zj\1118\tu\pl\sqtl.jpg white'
'enable print g:\zj\1118\tu\pl\sptl'
'print'
'disable print'
pull dummy
我画水汽通量和水汽通量散度的时间垂直剖面图，你可以这样试试噶

river

...该... 发表于 2013-3-19 21:26

                               
登录/注册后可看大图

我把高度都作了调整，还是做不出那样的图啊。不知道该怎么改啊。望各位提点宝贵的建议
'reinit'
'sdfopen ...

你要画水汽通量，就要先计算水汽通量，然后才固定经度，让时间和纬度变化再做剖面图。