【完美】解决FNL资料任意两点间斜剖面叠加地形的问题

enso816

完美！赞楼主无私分享！

天边的那边

为什么画不出地形阴影呢，运行也没有报错

river

天边的那边 发表于 2024-8-24 11:17
为什么画不出地形阴影呢，运行也没有报错

你把脚本贴上来看看，或者发个求助帖，详细描述一下你的画图过程

天边的那边

river 发表于 2024-8-24 23:58
你把脚本贴上来看看，或者发个求助帖，详细描述一下你的画图过程

'reinit'
'open orog0.25.ctl'
'open grib20240625.ctl'
'enable print along98E32N-2024062512-2712.gmf'
#   r       g       b##91 light yellow
'set rgb 82   36       0     216'
'set rgb 83   24      28     247'
'set rgb 84   40      87     255'
'set rgb 85   61     135     255'
'set rgb 86   86     176     255'
'set rgb 87  117     211     255'
'set rgb 88  153     234     255'
'set rgb 89  188     249     255'
'set rgb 90  234     255     255'
'set rgb 91  255     255     234'      
'set rgb 92  255     241     188'
'set rgb 93  255     214     153'
'set rgb 94  255     172     117'
'set rgb 95  255     120      86'
'set rgb 96  255      61      61'
'set rgb 97  247      39      53'
'set rgb 98  216      21      47'
'set rgb 99  165       0      33'

i=13
while(i<=49)

*先插值地形到所需斜线
'set dfile 1'
'set t 1'
'set x 1'
'set y 1'
'set z 1 17'

lon1 = 98.0
lon2 = 102
lat1 = 32.0
lat2 = 31

lon=lon1
'collect 2 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 2 gr2stn(orog.1,'lon','lat')'
lon=lon+0.25
endwhile

*再插值变量到所需斜线
'set dfile 2'
'set x 1'
'set y 1'
'set t 'i
'q time'
a=subwrd(result,3)
'set lev 600 100'

lon1 = 98.0
lon2 = 102
lat1 = 32.0
lat2 = 31
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(wprs.2*10,'lon','lat')'
  lon = lon + 0.25
endwhile

'set vpage 0 11 0 8.5'
'set parea 1 10 1.5 8'
'set xlopts 1 6 0.22'
'set ylopts 1 6 0.22'
'set font 0'
'set grid off'
'set grads off'
*'set csmooth on'
'set zlog on'
'set xlab on'
'set lon 98 102'
'set lev 600 100'
'set xlabs 32N,98E|31.75N,99E|31.5N,100E|31.25N,101E|31N,102E'
*'set ylevs 700 650 600 550 500 450 400 350 300 250 200 150 100'
*'set ylab on'
*'set clab on'
*画地形
'set gxout shade2'
'set clevs 0'
'set ccols -1 1'
'd coll2gr(2,-u)'
*画变量
'set gxout shaded'
*'color  -12 0 2 -kind blue->lightblue->white'
*'set ccols 0 0 7 0'
'set clevs -14 -12 -10 -8 -6 -4 -2 0'
'set ccols 82 84 85 86 87 88 89 90 0'
'd coll2gr(1,-u)'
'cbar 1 1'
'set strsiz 0.22'
'set string 1 l 6'
'draw string 3 8.4 'a' vvelprs'

'print'
'c'
i=i+1
endwhile
'disable print'
'reinit'

天边的那边

river 发表于 2024-8-24 23:58
你把脚本贴上来看看，或者发个求助帖，详细描述一下你的画图过程

可以出垂直速度图，没有地形。单独画地形剖面没问题，是有的。

river

天边的那边 发表于 2024-8-25 17:56
可以出垂直速度图，没有地形。单独画地形剖面没问题，是有的。

你这个问题应该是变量的填色图覆盖了地形的填色图，画图顺序反过来，先画变量再画地形试试吧。

天边的那边

river 发表于 2024-8-25 20:28
你这个问题应该是变量的填色图覆盖了地形的填色图，画图顺序反过来，先画变量再画地形试试吧。

是的，如果如此，非常感谢。
我还以为地形一定要先画。

river

天边的那边 发表于 2024-8-25 22:15
是的，如果如此，非常感谢。
我还以为地形一定要先画。

帖子没看仔细，是地形插值过程一定要放在前面，不是先画地形