（已解决）为什么用era5数据自己算的散度和数据自带的散度画出来的图不一样。

┏_(゜ω゜)=

'reinit'
*200hPa散度图
'open e:/era5/era.ctl'
'enable print e:/era5/5_21/div200_1.gmf'


'set mpdset cnworld'

t=1
while(t<=24)
'set t 't''
'set lev 200'
'set lon 90 130'
'set lat 0 40'  


*-------------------------------------------------------------
'set parea 1 9.4 1 8'
*设置画布尺寸
'set xlopts 1 4 0.14'
*设置X轴标注颜色 粗细 大小
'set ylopts 1 4 0.14'

'set grid off'
'set grads off'
'set xlint 5'
'set ylint 5'
'set gxout shaded'
'set frame on'

'set lon 100 125'
'set lat 15 30'  

'set rgb 20 255 255 255'
'set rgb 21 181 202 255'
'set rgb 22 143 179 255'
'set rgb 23 171 207 99'
'set rgb 24 255 250 20'
'set rgb 25 255 209 33'
'set rgb 26 255 76 0'

'set rbcols 20 21 22 23 24 25 26 '
'set clevs 0 2 4 6 8 10 '


'd VOprs*1e5'
*数据本身的变量
'run cbarn.gs'


*-------------------------------------------------------------------------------------

'set gxout stream'
'set strmden 3'
'd  Uprs;Vprs'

*--------------------------------------------------------------------------------------
'q dims'

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

'print'
'c'
t=t+1
endwhile
*-------------------------------------------------------------------------------------
'disable print'
'reinit'
*————————————————————————————————————————————————————————————————————
'reinit'
*200hPa散度图
'open e:/era5/era.ctl'
'enable print e:/era5/5_21/div200.gmf'


'set mpdset cnworld'

t=1
while(t<=24)
'set t 't''
'set lev 200'
'set lon 90 130'
'set lat 0 40'  

'define vor=hdivg(uprs,vprs)'

*-------------------------------------------------------------
'set parea 1 9.4 1 8'
*设置画布尺寸
'set xlopts 1 4 0.14'
*设置X轴标注颜色 粗细 大小
'set ylopts 1 4 0.14'

'set grid off'
'set grads off'
'set xlint 5'
'set ylint 5'
'set gxout shaded'
'set frame on'

'set lon 100 125'
'set lat 15 30'  

'set rgb 20 255 255 255'
'set rgb 21 181 202 255'
'set rgb 22 143 179 255'
'set rgb 23 171 207 99'
'set rgb 24 255 250 20'
'set rgb 25 255 209 33'
'set rgb 26 255 76 0'

'set rbcols 20 21 22 23 24 25 26 '
'set clevs 0 2 4 6 8 10 '


'd vor*1e5'
*自己算的散度
'run cbarn.gs'


*-------------------------------------------------------------------------------------

'set gxout stream'
'set strmden 3'
'd  Uprs;Vprs'

*--------------------------------------------------------------------------------------
'q dims'

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

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

┏_(゜ω゜)=

d voprs  画的图

┏_(゜ω゜)=

自己定义散度，用速度画出来的图

平流层的萝卜

hdfivg是散度吧。。。。

┏_(゜ω゜)=

平流层的萝卜 发表于 2020-10-31 14:25
hdfivg是散度吧。。。。

你说的对，照搬程序出错了，应该求的是散度。用公式画和用era5的变量画是没有区别的

┏_(゜ω゜)=

'reinit'
*200hPa散度图
'open e:/era5/era.ctl'
'enable print e:/era5/5_21/div200_1.gmf'


'set mpdset cnworld'

t=1
while(t<=24)
'set t 't''
'set lev 200'
'set lon 90 130'
'set lat 0 40'  


*-------------------------------------------------------------
'set parea 1 9.4 1 8'
*设置画布尺寸
'set xlopts 1 4 0.14'
*设置X轴标注颜色 粗细 大小
'set ylopts 1 4 0.14'

'set grid off'
'set grads off'
'set xlint 5'
'set ylint 5'
'set gxout shaded'
'set frame on'

'set lon 100 125'
'set lat 15 30'  

'set rgb 20 255 255 255'
'set rgb 21 181 202 255'
'set rgb 22 143 179 255'
'set rgb 23 171 207 99'
'set rgb 24 255 250 20'
'set rgb 25 255 209 33'
'set rgb 26 255 76 0'

'set rbcols 20 21 22 23 24 25 26 '
'set clevs 0 2 4 6 8 10 '


'd Dprs*1e5'
'run cbarn.gs'


*-------------------------------------------------------------------------------------

'set gxout stream'
'set strmden 3'
'd  Uprs;Vprs'

*--------------------------------------------------------------------------------------
'q dims'

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

'print'
'c'
t=t+1
endwhile
*-------------------------------------------------------------------------------------
'disable print'
'reinit'
正确画法

clarawho

┏_(゜ω゜)= 发表于 2020-11-02 20:01
你说的对，照搬程序出错了，应该求的是散度。用公式画和用era5的变量画是没有区别的

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

来自苹果APP客户端