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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 19271|回复: 25

[图形美化] 为什么用era5数据(已解决)画500hPa高度图会很怪?位势高度比较低,588线不见了

[复制链接]

新浪微博达人勋

发表于 2020-11-2 19:40:54 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 ┏_(゜ω゜)= 于 2020-11-5 10:23 编辑

'reinit'
*中国500hPa 高度图
'open e:/era5/era.ctl'
'enable print e:/era5/5_21/500hgt.gmf'
'set grads off'
'set grid off'

'set mpdset cnworld'

t=1
while(t<=24)
'set t 't''

'set lat 0 50'
'set lon 80 150'

'set ylint 10'
'set xlint 10'
'set ylopts 1 2 0.15'
'set xlopts 1 2 0.20'

'set lev 500'

'set gxout shaded'
'set rgb 20 255 255 255'
'set rgb 21 183 200 254'
'set rgb 22 127 150 254'
'set rgb 23 0 99 255'
'set rgb 24 0 198 51'
'set rgb 25 152 255 0'
'set rgb 26 255 255 1'
'set rgb 27 255 159 0'

'set rbcols 20 21 22 23 24 25 26 27 '
'set clevs 20 25 30 35 40 45 50'

'd mag(Uprs,Vprs)'
'run cbarn.gs'


*------------------------------------
****************500hpa高度场,单位:10gpm,间隔:2**********


'set gxout contour'
'set lev 500'
'set ccolor 1'
'set clopts -1 -1 0.08'
'set cthick 4'
'set cint 4'

'd Zprs/100'

'set clevs 588'
'set cthick 7'
'set ccolor 2'
'd Zprs/100'

'set clevs 586 '
'set cthick 6'
'set ccolor 1'
'd Zprs/100'

'set clevs 572 '
'set cthick 6'
'set ccolor 1'
'd Zprs/100'

'set clevs 556 '
'set cthick 7'
'set ccolor 1'
'd Zprs/100'

*************500hpa温度场,单位:度,间隔:4度*************
'set cthick 3'
'set cstyle 3'
'set ccolor 1'
'set cint 5'
'set clopts -1 -1 0.08'

'd Tprs-273.15'

*****************500hpa风场****************

'set gxout barb'

'set digsiz 0.05'
'set arrowhead 0.05'
'set ccolor 11'
'd skip(Uprs*2.5,8,8);skip(Vprs*2.5,8,8)'




*************************************
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
'draw title  'recs''

'print'

'c'
t=t+1
endwhile
*-------------------------------------------------------------------------------------
'disable print'
'reinit'


该参数是在特定位置相对于平均海平面的单位质量的重力势能。克服重力将单位质量从平均海平面提升到该位置也是要做的工作量。可以通过用地势除以地球的重力加速度g= 9.80665 m s-2)来计算地势高度。地势高度在天气学(天气模式分析)中起着重要作用。在恒定压力水平(例如300500850 hPa)下绘制的地势高度图可用于识别天气系统,例如旋风,反旋风,低谷和山脊。在地球表面,此参数显示表面的地势(高度)的变化,通常称为地形。 Zprs
500.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-11-9 19:50:53 | 显示全部楼层
一楼正解,不是/100,应该是/98
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2020-11-2 22:31:00 | 显示全部楼层
他高度场不是/10,是除9.8
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-11-3 15:28:26 | 显示全部楼层
t_imba 发表于 2020-11-2 22:31
他高度场不是/10,是除9.8

听说回复有金币,我试试
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-11-5 10:14:13 | 显示全部楼层
t_imba 发表于 2020-11-2 22:31
他高度场不是/10,是除9.8

你说的对应该除以9.8再除以10
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-11-6 16:41:51 | 显示全部楼层
以前没发现,学到了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-11-7 09:27:32 | 显示全部楼层
本帖最后由 ┏_(゜ω゜)= 于 2021-1-7 15:31 编辑

修改后

'reinit'
*中国500hPa 高度图
'open e:/era5/era.ctl'
'enable print e:/era5/5_21/500hgt.gmf'
'set grads off'
'set grid off'

'set mpdset cnworld'

t=1
while(t<=24)
'set t 't''

'set lat 0 50'
'set lon 80 150'

'set ylint 10'
'set xlint 10'
'set ylopts 1 2 0.15'
'set xlopts 1 2 0.20'

'set lev 500'

'set gxout shaded'
'set rgb 20 255 255 255'
'set rgb 21 183 200 254'
'set rgb 22 127 150 254'
'set rgb 23 0 99 255'
'set rgb 24 0 198 51'
'set rgb 25 152 255 0'
'set rgb 26 255 255 1'
'set rgb 27 255 159 0'

'set rbcols 20 21 22 23 24 25 26 27 '
'set clevs 20 25 30 35 40 45 50'

'd mag(Uprs,Vprs)'
'run cbarn.gs'


*------------------------------------
****************500hpa高度场,单位:10gpm,间隔:2**********


'set gxout contour'
'set lev 500'
'set ccolor 1'
'set clopts -1 -1 0.08'
'set cthick 4'
'set cint 4'

'd Zprs/9.8/10'

'set clevs 588'
'set cthick 7'
'set ccolor 2'
'd Zprs/9.8/10'

'set clevs 586 '
'set cthick 6'
'set ccolor 1'
'd Zprs/9.8/10'

'set clevs 572 '
'set cthick 6'
'set ccolor 1'
'd Zprs/9.8/10'

'set clevs 556 '
'set cthick 7'
'set ccolor 1'
'd Zprs/9.8/10'

*************500hpa温度场,单位:度,间隔:4度*************
'set cthick 3'
'set cstyle 3'
'set ccolor 1'
'set cint 4'
'set clopts -1 -1 0.08'

'd Tprs-273.15'

*****************500hpa风场****************

'set gxout barb'

'set digsiz 0.05'
'set arrowhead 0.05'
'set ccolor 11'
'd skip(Uprs*2.5,8,8);skip(Vprs*2.5,8,8)'




*************************************
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
'draw title  'recs''

'print'

'c'
t=t+1
endwhile
*-------------------------------------------------------------------------------------
'disable print'
'reinit'
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-11-7 12:39:28 | 显示全部楼层
┏_(゜ω゜)= 发表于 2020-11-7 09:27
'reinit'
*中国500hPa 高度图
'open e:/era5/era.ctl'

再上个图,对比一下就完美了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-11-7 16:49:28 | 显示全部楼层
500HP高度图
500.jpg
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-11-21 20:39:38 | 显示全部楼层
我之前也遇到了这个问题 是位势高度算法不一样 楼上回答的都对
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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