- 积分
- 6328
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-26
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-5-18 16:49:09
|
显示全部楼层
完整gs如下:
'reinit'
'open d1.ctl'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.15'
'set zlog on'
*'set csmooth on'
'set t 17'
'set x 1'
'set y 1'
*'set lon 111 118'
*'set lat 30 33'
'set z 1 19'
lon1 = 101.5
lon2 = 102.2
lat1 = 36.65
lat2 = 35.25
lon = lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(dbz,'lon','lat')'
lon =lon+0.05
endwhile
'set cmin 0.01'
'set t 17'
'set lon 101.5 102.2'
'set xlint 0.1'
'set z 1 19'
*'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
*'set ccols 0 0 7 0'
'd coll2gr(1,-u)'
'cbarn_interp.gs'
'set cmin 0.01'
'set gxout contour'
'set cthick 1'
'd coll2gr(1,-u)'
'printim qxiepo.png white'
'reinit' |
|