请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4797|回复: 6

[图形美化] 为什么我用era5数据画斜剖图没有地形?

[复制链接]

新浪微博达人勋

发表于 2020-11-13 11:38:31 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
照着别人的帖子画的。
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'






360截图20201113113900536.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-11-28 20:12:26 | 显示全部楼层
你这斜线在地图上画一下,看看它到底经不经过某些地形?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-1-11 18:03:01 | 显示全部楼层
俺也遇到了相同的问题,...
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-1-11 18:06:36 | 显示全部楼层
循环中的lon=lon+1,要不要改为lon=lon+0.25?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-1-14 18:49:58 | 显示全部楼层
你用的哪个地图文件?NCEP0.25*0.25的地形太粗了,
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-2-9 17:52:41 | 显示全部楼层
你好,请问解决了没?我也是剖面图画不出era5地形来
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-3 11:51:54 | 显示全部楼层
你好,请问问题解决了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表