- 积分
- 117
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
上面是斜剖面,下面是39N剖面,看风向杆怎么杆的位置相反呢
lon1=116.0
lon2=121.0
lat1=41.0
lat2=35.0
'define alfa=atan2('lat1-lat2','lon2-lon1')'
'set x 1'
'set y 1'
'set z 1 20'
'set t 1'
lon=lon1
'collect 1 free'
'collect 2 free'
'collect 3 free'
while(lon<=lon2)
lat=lat2+(lat1-lat2)*(lon2-lon)/(lon2-lon1)
*将格点插值到站点(lon,lat)上
'collect 1 gr2stn(rhprs,'lon','lat')'
'collect 2 gr2stn(UGRDprs,'lon','lat')'
'collect 3 gr2stn(VGRDprs,'lon','lat')'
'collect 4 gr2stn(HGTsfc-hgtprs,'lon','lat')'
lon=lon+1
endwhile
'set lon 116 121'
'set z 1 13'
*'set ylevs 1000 925 850 700 600 500 400 300 200'
'set xlabs 41N,116E|39.8N,117E|38.6N,118E|37.4N,119E|36.2N,120E|35N,121E'
'set gxout shaded'
*'set csmooth on'
*'set clab on'
'd coll2gr(1,-u)'
'cbarn'
'set gxout barb'
'set cthick 6'
'd coll2gr(2,-u)*2.5;coll2gr(3,-u)*2.5'
沿任意两点斜剖面
沿39N剖面
|
|