grads画任意方向的垂直剖面问题

happymgj

下面是CTL：
dset F:/Paper0808/Radar_2_0808/Mref_Out/1708%d2%h2%n2.grd
options template
title radar data
undef 999.000
xdef  841   linear 113.470  0.01
ydef 757   linear  35.310  0.01
zdef 20  levels 500 1000 1500 2000 2500 3000 3500 4000 5000 6000 7000 8000 9000 10000 11000 12000 13000 14000 15000 16500
tdef 120   linear 12:06Z08AUG17 6mn
vars 1
dbz 21  99  dbz
endvars

下面是GS：
'reinit'
'open F:\Paper0808\Radar_2_0808\Mref_Pic\Mref.ctl'

'set t 30'
*剖面处理
'set lat 39'
'set lon 117'
'set z 1 18'

lon1=116
lon2=117
lat1=38
lat2=39
*剖面的角度
*'define alfa=atan2('lat2-lat1','lon2-lon1')'

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.01
endwhile
*以上是剖面插值

*页面设置
'set font 0'
'set parea 3 9 3 6'
'set grid off'
'set grads off'
'set xlopts 1 3 0.16'
'set ylopts 1 3 0.16'
'set lon 116 117.479'
'set xlint 0.2'

'set z 2 18'
'set ylint 50'
'set zlog on'

'd coll2gr(1,-u)'

'print'
'disable print'
画任意高度啥的都没问题，望各位老师指点迷津

river

你这是自己编的，还是按照别人的帖子编的？如果是按照别人的帖子编的，那就请仔细看看原贴

happymgj

river 发表于 2017-10-15 15:16
你这是自己编的，还是按照别人的帖子编的？如果是按照别人的帖子编的，那就请仔细看看原贴

这是我自己编的，数据是我自己做的，gs没有没问题，不过ctl改了很多东西，已经搞定了