为什么用era5数据（已解决）画500hPa高度图会很怪？位势高度比较低，588线不见了

┏_(゜ω゜)=

'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）来计算地势高度。地势高度在天气学（天气模式分析）中起着重要作用。在恒定压力水平（例如300、500或850 hPa）下绘制的地势高度图可用于识别天气系统，例如旋风，反旋风，低谷和山脊。在地球表面，此参数显示表面的地势（高度）的变化，通常称为地形。 Zprs

┏_(゜ω゜)=

修改后

'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'

hbl1021

一楼正解，不是/100，应该是/98

t_imba

他高度场不是/10，是除9.8

/奋斗

t_imba 发表于 2020-11-2 22:31
他高度场不是/10，是除9.8

听说回复有金币，我试试

┏_(゜ω゜)=

t_imba 发表于 2020-11-2 22:31
他高度场不是/10，是除9.8

你说的对应该除以9.8再除以10

麻花

以前没发现，学到了

river

┏_(゜ω゜)= 发表于 2020-11-7 09:27
'reinit'
*中国500hPa 高度图
'open e:/era5/era.ctl'

再上个图，对比一下就完美了

┏_(゜ω゜)=

500HP高度图

zz460908990

我之前也遇到了这个问题 是位势高度算法不一样 楼上回答的都对