为什么我用era5数据画斜剖图没有地形？

┏_(゜ω゜)=

照着别人的帖子画的。
https://bbs.06climate.com/forum. ... hlight=%D0%B1%C6%CA下载了0.25*0.25地形数据，修改了ctl的路径。
https://bbs.06climate.com/forum. ... hlight=%D0%B1%C6%CA

'reinit'
*斜剖
'open e:/era5/era.ctl'
'open e:/era5/orog0.25.ctl'
*注意修改ctl文件路径

'enable print e:/era5/5_21/xiepou.gmf'


*先插值地形到所需斜线

'set dfile 2'
'set t 1'
'set x 1'
'set y 1'
'set z 1 17'

*—————————————修改
lon1=110
lon2=130
lat1=20
lat2=24
*--------------------------
lon=lon1
'collect 3 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 3 gr2stn(orog,'lon','lat')'
lon=lon+1
endwhile


*再插值变量到所需斜线
'set dfile 1'
*________________________________根据ctl修改
'set x 1 361'
'set y 1 241'
'set z 1 21'

'set t 18'

*________________________________
*批量描述FNL，多时次
'define alfa=atan2('lat2-lat1','lon2-lon1')'
*atan2返回的结果单位是弧度
'define uv=Uprs*cos(alfa)+Vprs*sin(alfa)'

'set x 1'
'set y 1'
'set z 1 17'
*修改
*——————————————————
*修改
lon1=110
lon2=130
lat1=20
lat2=24
*——————————————————
lon=lon1

'collect 1 free'
'collect 2 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(-Wprs*100,'lon','lat')'
'collect 2 gr2stn(uv,'lon','lat')'
lon=lon+1
endwhile

'set grads off'
'set zlog on'
'set grid off'
*____________________________
'set lon 110 130'
*____________________________

'set z 1 17'

*——————————————
'set ylab off'
'set xlab off'
'set csmooth on'
*画地形
'set clevs 0'
'set ccols 0 1 0'
'set gxout shaded'
'd coll2gr(3,-u)'
*画变量
'set ylab on'
'set xlab on'
'set csmooth on'

*'set ylevs 1000 925 850 700 600 500 400 300 200 100'
*'set xlabs 35N,70E|37N,82E|39N,94E|41N,106E|43N,118E|45N,130E'
*'set xaxis 'lon1' 'lon2''

'set gxout stream'
'set clab on'
'd coll2gr(2,-u);coll2gr(1,-u)'
'draw title along (35N,70E) to (45N,130E)'


'print'
'disable print'

'reinit'

river

你这斜线在地图上画一下，看看它到底经不经过某些地形？

jeanne2006

俺也遇到了相同的问题，...

jeanne2006

循环中的lon=lon+1，要不要改为lon=lon+0.25?

jeanne2006

你用的哪个地图文件？NCEP0.25*0.25的地形太粗了，

sysue16

你好，请问解决了没？我也是剖面图画不出era5地形来

视线之外

你好，请问问题解决了吗