作物理量诊断分析时部分物理量gs文件

yanyan330

垂直速度：
'open d:\fnl\fnl.ctl'
'set mpdset cnworld cnriver cnbase'
"set lon 105 120"
"set lat 25 35"
"set lev 850"
"set t 1"
'd VVELprs'
"printim D:\fnl\082008.gif white"
'disable print'
'reinit'
散度：
'open d:\fnl\fnl.ctl'
'set mpdset cnworld cnriver cnbase'
"set lon 105 120"
"set lat 25 35"
"set lev 850"
"set t 1"
'd hdivg(UGRDprs,VGRDprs)*100000'
"printim D:\fnl\082008.gif white"
'disable print'
'reinit'
涡度：
'open d:\fnl\fnl.ctl'
'set mpdset cnworld cnriver cnbase'
"set lon 105 120"
"set lat 25 35"
"set lev 850"
"set t 1"
'd hcurl(UGRDprs,VGRDprs)*100000'
"printim D:\fnl\082008.gif white"
'disable print'
'reinit'
水汽通量散度：
'open d:\fnl\fnl.ctl'
'set mpdset cnworld cnriver cnbase'
'set lon 1 360'
'set lat 1 90'
'set lev 1000 100'
*'define pres=1.333*760*pow((1-hgtprs/44300),5.256)'
'define es=6.112*exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))'
'define qs=0.62197*es/(lev-0.378*es)*1000'
'define q=qs*RHprs/100.0'
'define q1=q*ugrdprs/9.8'
'define q2=q*vgrdprs/9.8'
'define qfdiv=hdivg(q1,q2)*1e6'
'set lat 25 40'
'set lon 115'
'set lev 1000 100'
'set t 4'
'set clopts -1 -1 0.15'
'set xlopts 1 4 0.20'
'set ylopts 1 4 0.20'
'set ccolor 1'
'set cmin 0'
'set cstyle 1'
'set csmooth on'
'set cint 4'
'd qfdiv'
'set cmax 0'
'set cstyle 4'
'set csmooth on'
'set cint 8'
'd qfdiv'
'printim d:\fnl\qfdiv.gif white'

水汽通量：
'open d:\fnl\fnl.ctl'
'set mpdset cnworld cnriver cnbase hub'
'set lon 70 130'
'set lat 15 60'
'set lev 700'
'define pres=1.333*760*pow((1-hgtprs/44300),5.256)'
'define es=6.112*exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))'
'define qs=0.62197*es/(pres-0.378*es)*1000'
'define q=qs*RHprs/100.0'
'define q1=q*ugrdprs/9.8'
'define q2=q*vgrdprs/9.8'
'define qflux=mag(q1,q2)'
'set grid off'
'set gxout contour'
'set ccolor 1'
'set cint 2'
'd qflux'
'set gxout barb'
*'d ugrdprs;vgrdprs'
'printim d:\fnl\qflux700hpa.gif white'

topmad

感谢楼主的分享，要是有附图就更漂亮了

善人/jw

感谢LZ，可以弄个集嘛~~

meehooqq

楼主用的是不是ncep的6h的fnl数据？

Donald

很好，学习啦！

luoziwuhui

非常不错的东西，对我帮助很大，感谢楼主大多无私分享！

小肥虾619

楼主好人，但是我想问一下水汽通量散度那儿计算q时为什么用lev呢？求解释

疯一样的男人

谢谢分享！

Forcast

小肥虾619 发表于 2012-6-7 18:06

                               
登录/注册后可看大图

楼主好人，但是我想问一下水汽通量散度那儿计算q时为什么用lev呢？求解释

LEV就是气压，应该这样理解吧

小肥虾619

Forcast 发表于 2012-6-8 07:57

                               
登录/注册后可看大图

LEV就是气压，应该这样理解吧

好的好的，谢谢楼主了