求一个计算某区域平均水汽通量的时间剖面图的gs程序

黑鹰

我要画一个某区域的从地面开始至高空200hpa的平均水汽通量随时间变化的剖面图，求高手给指导下啊！！！我有温度场、相对湿度场、比湿场以及风场资料！！非常感谢喽

云大学子

楼主，你弄出来了吗？能传授一下吗？我现在也是在弄这个问题

黑鹰

云大学子 发表于 2015-5-4 14:08
楼主，你弄出来了吗？能传授一下吗？我现在也是在弄这个问题

呵呵，不好意思，好久没有来了。出差去了，都没有弄。还没有整出来呢，你弄出来了可以分享下哈

云大学子

弄出来了，*水汽通量散度随时间变化图
'reinit'
'sdfopen d:\zpw\nc1.nc'
'set display color white'
'c'
'set grid off'
'set grads off'
'set csmooth on'
'set zlog on'
*等值线及坐标轴的设置
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set clopts 1 4 0.14'
'set lon 90 120'
'set lat 20 30'
'set z 1 11'
'set t 12 45'
*计算水汽通量散度的过程
*计算出饱和水汽压
'define e=(6.112*exp((17.67*(t-273.16))/(t-29.65)))'
*'define e=(6.1078*exp((17.13*t-273.16)/(t-38)))'
*计算出饱和比湿
'define qs=(0.622*e/(lev-0.378*e))'
*计算出比湿
'define q=qs*r/100'
'define qdiv=hdivg(q*u/9.8,q*v/9.8)*1000'

'set x 1'
'set y 1'
'define ww=aave(qdiv,lon=104,lon=109.5,lat=26,lat=28)'
'set ylevs 1000 950 900 850 800 750 700 650 600 550 500'
*'set xlabs 25|26|27|28|29|30'
'set xlabs |01-16|01-21|01-26|02-01|02-06|02-11|'
'set gxout shaded'
'set cmax 0'
'd ww'
'set ylpos 0 r'
'set gxout contour'
*'set cint 4'
'd ww'
'cbarn'
'printim d:\zpwsandushijian.jpg white'

黑鹰

云大学子 发表于 2015-5-17 18:57
弄出来了，*水汽通量散度随时间变化图
'reinit'
'sdfopen d:\zpw\nc1.nc'

太感谢了，呵呵，正在努力的发愁中呢，太感谢了

晓筱

非常感谢楼主！

wdd1992

学习了。。。。。