用grads绘制沿任意折线的剖面

竹下客

学习http://bbs.06climate.com/forum.p ... rads%2B%C6%CA%C3%E6这个帖子，经过多次实验，成功绘出沿折线的剖面图。

river

直接贴上得了，下载好麻烦······

river

lon=lon1
while(lon<lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(-100*vvelprs.1,'lon','lat')'
'collect 2 gr2stn(UGRDprs.1*cos(alfa)+VGRDprs.1*sin(alfa),'lon','lat')'
'collect 3 gr2stn(hh.2(t=1),'lon','lat')'

lon=lon+1
endwhile

lon=lon2
while(lon>=lon3)
lat=lat2+(lat3-lat2)*(lon-lon2)/(lon3-lon2)
'collect 1 gr2stn(-100*vvelprs.1,'lon','lat')'
'collect 2 gr2stn(UGRDprs.1*cos(alfa2)+VGRDprs.1*sin(alfa2),'lon','lat')'
'collect 3 gr2stn(hh.2(t=1),'lon','lat')'

lon=lon-1
endwhile

我想请教一下，你是怎么确认这两个collect 1 和 collect 2 的数据是追加而不是覆盖的？我有点懵

竹下客

river 发表于 2019-9-6 22:27
lon=lon1
while(lon=lon3)
lat=lat2+(lat3-lat2)*(lon-lon2)/(lon3-lon2)

我实验的。先只画第一截，然后只画第二截，最后两截都画。发现它跟串珠一样的。collect几个站点数据，就依次串几颗珠子，最后coll2gr就是拎着绳子尾端一颗一颗按顺序释放出来。

river

竹下客 发表于 2019-9-9 09:08
我实验的。先只画第一截，然后只画第二截，最后两截都画。发现它跟串珠一样的。collect几个站点数据，就 ...

厉害了！回头有空我也试试