wrf_user_intrp3d插值位势高度到850hpa后出现数据空白，如何解决？

happymgj

如图，位势高度场有些区域是空白
脚本相关部分如下：
----------------------------------------------------------------------------------------------
f=addfile("E:/a/201311/em_real11/wrfout_d01_2013-11-25_00_00_00.nc","r")
hgt   = wrf_user_getvar(f,"z" ,0)
p  = wrf_user_getvar(f,"pressure",0) ; total pressure
pressure = 850.   ; 850 hPa
hgt850_wrf =   wrf_user_intrp3d(hgt(:,:,:),p(:,:,:),"h",pressure,0.,False)
if(any(ismissing(hgt850_wrf))) then
    print("Your data contains some missing values. Beware.")
  end if

-----------------------------------------------------------------------------------------------
ncl显示 (0)     Your data contains some missing values. Beware.

我认为是地形对插值造成了影响，又把地形画上去看了一下

好像地形确实影响了等高线插值。但是在fnl资料中，850hpa并没有缺测空白。
那我该如何把wrfout结果画出的有缺测的空白区域“订正”成有正常数据的情况呢？
ARWPOST会出现这种情况吗？

风之牧语

缺测试正常的，在气压层里面，wrfout数据是没有地面以下，NCL函数在插值的时候不会插值出地面以下的。你看看WRF模式的网格设定就知道了。

happymgj

风之牧语 发表于 2017-12-8 17:10
缺测试正常的，在气压层里面，wrfout数据是没有地面以下，NCL函数在插值的时候不会插值出地面以下的。你看 ...

谢谢。那像MICAPS资料，FNL资料中，850hpa高度场在我国区域内没有缺测值，连青藏高原也没有缺测。是不是因为：此“850hpa”并非实际气压，而是通过某种“订正”得到的。那么wrfout数据有没有这种订正呢？

风之牧语

happymgj 发表于 2017-12-10 09:36
谢谢。那像MICAPS资料，FNL资料中，850hpa高度场在我国区域内没有缺测值，连青藏高原也没有缺测。是不是 ...

这个和模式系统的网格设定有关系。再说青藏高原虽然没有缺测，但是结果是虚假的啊，有反而麻烦，作图时还要去屏蔽。

happymgj

风之牧语 发表于 2017-12-15 17:42
这个和模式系统的网格设定有关系。再说青藏高原虽然没有缺测，但是结果是虚假的啊，有反而麻烦，作图时还 ...

明白了 谢谢！

supersupergirl

请问楼主如何把地形也叠加高度场的呢

小张ccc

你好，请问可以分享一下完整代码吗？我在绘制图像时存在问题，想参考一下

Hannah_YCu

请问楼主缺测值如何订正