wrfout中U、V、T、Q、PH变量插值到850hpa、500hpa和200hpa等压面出现缺测值

弟弟救我

各位大佬，
       wrfout中U、V、T、Q、PH变量插值到850hpa、500hpa和200hpa等压面，画出图来发现在850hpa上所有变量都出现了缺测，并且缺测的区域也是固定的，代码和图贴上来。有遇到过类似情况的么？
         diri = "/nuist/p/work/zhanghp/bgm/inputd02/"
   files = systemfunc("ls "+diri+"wrfout_d02_2015-07-26_00:00:00*")+ ".nc"
   n = addfiles(files,"r")
    xlat = wrf_user_getvar(n,"XLAT",0)
    xlong =wrf_user_getvar(n,"XLONG",0)
    p= wrf_user_getvar(n,"pressure",1) ; total pressure
    u = wrf_user_getvar(n,"U",-1)
    v = wrf_user_getvar(n,"V",-1)
    t = wrf_user_getvar(n,"T",-1)
    q = wrf_user_getvar(n,"QVAPOR",-1)
    ph = wrf_user_getvar(n,"PH",-1)
print(files)

u850=new((/15,south_north,west_east/),"float")
v850=new((/15,south_north,west_east/),"float")
t850=new((/15,south_north,west_east/),"float")
q850=new((/15,south_north,west_east/),"float")
ph850=new((/15,south_north,west_east/),"float")
do k=0,14
   u850(k,:,:) =wrf_user_intrp3d(u(2*k+1,:,:,0:197),p,"h",850,0.,False)
   v850(k,:,:) =wrf_user_intrp3d(v(2*k+1,:,0:167,:),p,"h",850,0.,False)   
   t850(k,:,:) =wrf_user_intrp3d(t(2*k+1,:,:,:),p,"h",850,0.,False)
   q850(k,:,:) =wrf_user_intrp3d(q(2*k+1,:,:,:),p,"h",850,0.,False)
   ph850(k,:,:) =wrf_user_intrp3d(ph(2*k+1,1:34,:,:),p,"h",850,0.,False)  
end do
u850@lat2d =xlat
u850@lon2d =xlong
uu=gsn_open_wks("ps","u850_do2")
         res=True
         res@mpInlandWaterFillColor = 0
         res@mpOutlineOn                 = True     
         res@mpDataBaseVersion       = "MediumRes"
         res@mpDataSetName           = "Earth..4"  
         res@mpOutlineBoundarySets = "AllBoundaries"        
         res@gsnMaximize          = True                             ; maximize plot in frame
         res@cnFillOn             = True                             ; turn on contour fill
         res@cnLinesOn            = False                             ;turn off contour lines
         res@cnLineLabelsOn       = False                            ; turn off line labels
        
          res@gsnFrame= True
         
          res@tfDoNDCOverlay       =True
          res@gsnAddCyclic          = False                           ; set to False if plotting regional data        
          res@mpMinLatF     = min(u850@lat2d)
         res@mpMaxLatF     = max(u850@lat2d)
          res@mpMinLonF     = min(u850@lon2d)
         res@mpMaxLonF     = max(u850@lon2d)
          res@mpProjection  = "CylindricalEquidistant"                ; The default
          res@pmTickMarkDisplayMode = "Always"                        ; nicer tickmarks
          res@mpGeophysicalLineThicknessF= 2.              
          res@mpNationalLineThicknessF= 2.                  
          res@lbOrientation="Vertical"
         
;         res = wrf_map_resources(nfile,res)
        plot = gsn_csm_contour_map(uu,u850(7,:,:),res)


file:///C:/Users/zhp/Documents/Tencent%20Files/1204036395/Image/C2C/GNS6)]3TDOO2%7DECYH(8%7B4K1.png

Liujiaxin

敢问：
u850=new((/15,south_north,west_east/),"float")
v850=new((/15,south_north,west_east/),"float")
t850=new((/15,south_north,west_east/),"float")
q850=new((/15,south_north,west_east/),"float")
ph850=new((/15,south_north,west_east/),"float")
中的south_north和west_east分别是什么啊？好像需要设定？

弟弟救我

还有就是在500hpa和200hpa上都是正常的

freekiller

固定有什么奇怪。你不觉得是地形吗

弟弟救我

freekiller 发表于 2018-5-4 13:46
固定有什么奇怪。你不觉得是地形吗

是因为地形海拔高然后不存在850hpa的等压面么？

astiny

九成是地形问题吧。

弟弟救我

astiny 发表于 2018-5-4 15:31
九成是地形问题吧。

十成是地形的问题，谢了

曼联现在没$

不是地形问题，我十一加班8天

Liujiaxin

p= wrf_user_getvar(n,"pressure",1) ; total pressure 似乎这里有什么问题，导致后面
u850(k,:,:) =wrf_user_intrp3d(u(2*k,:,:,0:99),p,"h",850,0.,False)找不到变量p

ZP盛夏光年

楼主的这个色标好好看，是哪个，可以分享一下吗