- 积分
- 21490
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-5-16
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-11-10 19:18:52
|
显示全部楼层
不是很清楚错在哪里
'reinit'
'open G:\swv\thetase\thatase20151107\thetase_all.ctl'
'set grads off'
'set grid off'
*'set zlog on'
*'set csmooth on'
*'set font 0'
'set parea 1 9 1 8'
'set x 1'
'set y 1'
lon1=101.25
lon2=107.25
lat1=27.09
lat2=35.01
'set lev 1000 100'
'set t 22'
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(thetase,'lon','lat')'
lon=lon+0.3
endwhile
'set lon 101.25 107.25'
'set lev 1000 100'
'set ylevs 1000 950 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100'
'set xlabs 27.09N,101.25E|29.07N,102.75E|31.05N,104.25E|33.03N,105.75E|35.01N,107.25E'
'set t 22'
'set gxout contour'
'd smth9(coll2gr(1,-u))'
'printim G:\swv\thetase\thatase20151107\poumian(thetase_all).png white'
; |
|