绘制海平面气压时在青藏高原上空不正常是怎么回事

ningning

用ERA5平均海平面气压画图，青藏高原地区等值线密集，很多小圈圈，不知是为啥？该如何解决呢？还有一个问题是底图怎么去除省界只保留国界呢？'reinit'
'sdfopen  D:\data\hgt\mean_sea_level_pressure.nc'
'set lat 20 50'
'set lon 80 130'
'set mpdset cnworld'
*'set poli off'
*'set map 1 1 4'
'draw map'
'set vpage 0 11 0 8.5'     
'set parea 1.2 10.3 1.0 8.0'  

i=1
while(i<=48)
'set t 'i''
'set grads off'
'set grid off'
'set xlpos -10'
'set ylpos -10'

'set gxout contour'
'set ccolor 1'
'set cthick 4'
'set cint 2.5'
'set clskip 2'
'd smth9(msl/100)'

'run D:\APPS\opengrads\Classic\scripts\axis.gs -type b  -position m -label on -interval 10  -lfont 5 -lsize 0.15 -suffix `3.'
'run D:\APPS\opengrads\Classic\scripts\axis.gs -type L  -position o -start 20 -end 60 -interval 10 -sinterval 1 -lfont 5 -lsize 0.15 -suffix `3.'
'gxprint D:\data\hgt\msl.png white'
'c'
i=i+1
endwhile
'reinit'
;

我画的图和micaps中同时次地面图对比下，可以看到区别

付亚男

这grads脚本看着无比的熟悉


再分析资料里海平面气压场在高原上是不准确的。
micaps里的线条是国家气象中心的预报员手动画的，按照天气分析中画等压线的原则，对高原进行了特殊的处理，一般会等压线尽量绕过高原，所以在高原上两个看起来不一样。

ningning

付亚男 发表于 2022-4-11 22:01
这grads脚本看着无比的熟悉

原来是这样！那应该怎么处理呢，用别的资料或者将青藏高原遮掩下？

栤唥de嶶笑

ningning 发表于 2022-4-13 08:53
原来是这样！那应该怎么处理呢，用别的资料或者将青藏高原遮掩下？

一、去除省界限，保留国界。只需要将'set mpdset cnworld'，改为'set mpdset lowres'。

二、可以用地形资料，或者西藏省界来遮掉青藏高原部分。

地形资料的话，函数这么写maskout(var1,var2-3000)，即把海拔3000米以上遮掉。var1指海平面气压；var2指地形高度。（要补充说明的是，地形资料的分辨率需要跟你气压资料分辨率一致，才能用；起止范围也建议保持一致）

西藏省界的话，函数这么写maskout(var,-mask)来实现。
用省界资料需要单独制作maskout文件，可用meteinfo软件来制作。
（maskout文件的分辨率需要跟你气压资料分辨率一致；起止范围也建议保持一致）

ningning

栤唥de嶶笑 发表于 2022-4-13 10:40
一、去除省界限，保留国界。只需要将'set mpdset cnworld'，改为'set mpdset lowres'。

二、可以用地 ...

非常感谢您的详细解答！
第一个问题试了下，‘set  mpdset lowres’国界也去掉了，保留了海岸线；改成‘set mpdset mres’可以保留国界；maskout我正在尝试

myj0219

111111111111111111111111111111