- 积分
- 55
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-20
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-10-30 16:13:54
|
显示全部楼层
本帖最后由 may6117 于 2015-10-30 16:15 编辑
'reinit'
'open D:\NCEP\xm\20090728\fnl_20090727.ctl'
'set mpdset hires cnworld'
'set map 1 1 1'
*'set mpdset neimenggu neimengguq'
i=1
while (i<=12)
'set t ' i
'set lon 90 135'
'set lat 30 55'
'set lev 850'
'set xlopts 1 1 0.2'
'set ylopts 1 1 0.2'
'set grads off'
'set grid off'
*'set xlopts -10'
*'set ylopts -10'
'set cthick 6'
'set ccolor 1'
'set gxout contour'
'set cint 10'
'set ccolor 1'
'set cthick 6'
'set clopts 1 6 0.12'
'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 qx=q*UGRDprs/9.8'
'define qy=q*vGRDprs/9.8'
'define qdiv=hdivg(qx,qy)'
'set lon 98 118'
'set lat 35 48'
'd qdiv*1e6'
'q w2xy 111.7 40.8'
x=subwrd(result,3)
y=subwrd(result,6)
'draw string 'x' 'y' H'
'q w2xy 40.8 111.7'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.1'
'q time'
res = subwrd(result,3)
year= substr(res,9,4)
'q time'
res = subwrd(result,3)
hour = substr(res,1,2)
say hour
'q time'
res = subwrd(result,3)
day=substr(res,4,2)
say day
'q time'
res = subwrd(result,3)
month=substr(res,6,3)
say month
'set strsiz 0.3 0.3'
'set string 1 c 6'
'draw title 'year'.'month'.'day'.'hour' 850hPa qdiv'
'printim D:\NCEP\xm\result\'year'.'month'.'day'.'hour'.qdiv850.png white'
'c'
i=i+1
endwhile
;
|
|