- 积分
- 201
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
来自:http://hzysara.blog.163.com/blog/static/142648576201061111856281/
这个是用NCEP资料画温度平流和水汽通量叠加的,不知道有用没
'reinit'
'sdfopen e:\surface\uwnd.2008.nc'
'sdfopen e:\surface\vwnd.2008.nc'
'sdfopen e:\surface\rhum.2008.nc'
'sdfopen e:\surface\air.2008.nc'
'set map 1 1 1'
'set mpdraw on'
'set poli on'
'set mproj latlon'
'set clab forced'
'set mpdset cnworld cnriver'
'set cthick 6'
'set lon 95 130'
'set lat 25 50'
'set lev 850'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.15'
'set clopts 1 1 0.15'
'set xlint 5'
'set ylint 5'
'set t 905'
'set grads off'
'set grid off'
'set gxout shaded'
'set clevs 1 2 4 6 8 10 12 14 16 18'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'd mag(q*uwnd.1/g,q*vwnd.2/g)*1000'
*'set gxout vector'
*'d q*uwnd.1/g;q*vwnd.2/g'
'set cthick 6'
'set ccolor 1'
'set gxout barb'
'd uwnd.1;vwnd.2'
'run cbarn'
*'d uwnd.1/g;vwnd.2/g;mag(uwnd.1,vwnd.2)'
'set gxout contour'
'define dtx=cdiff(air.4-273.15,x)'
'define dty=cdiff(air.4-273.15,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'set ccolor 2'
'set cstyle 2'
*'set cint 1'
'set csmooth on'
'd -1*((uwnd.1*dtx)/(cos(lat*3.1416/180)*dx)+vwnd.2*dty/dy)/6.37e6*100000'
'printim e:\surface\theta\a\sdtl2008081408-850hPa-3.png white'
;
|
评分
-
查看全部评分
|