- 积分
- 808
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 wty118 于 2016-5-8 21:47 编辑
'reinit'
'open e:\grib\201604\fnl_201604.ctl'
'enable print e:\picture\thetase(huanghua).gmf'
tt=57
while(tt<=61)
'set t 'tt
'set lat 28.18'
'set lon 113.22'
'set lev 1000 200'
'set xaxis 315 355'
'define prs=lev'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set font 1'
'define es=(6.112*exp(17.67*(Tmpprs-273.15)/(Tmpprs-29.65)))'
'define q=RHprs*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(Tmpprs)-log(e)-4.805)'
'define theta=Tmpprs*pow((1000./prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set gxout contour'
'set ccolor 1'
'd eqt'
'q time'date=subwrd(result,3)
'draw title thetase(huanghua)at 'date''
'print'
'c'
tt=tt+1
endwhile
'disable print'
'reinit'
用set xaxis命令设置有坐标数据范围以后出图和实际数据不一致,只是坐标改变了,图并没有改变,这是为什么?
|
|