使用wrfout数据绘制CIN失败

SG哭晕在厕所

按照官网说明使用wrf_cape_2d这个函数，出来的结果只有第三维的LCL和第四维的LFC，而第一维和第二维的CAPE和CIN全部显示为缺测值，非常费解，不知道怎么解决，跪求坛子里各位大神帮忙看看
以下是脚本，核心部分都是直接copy官网的例子
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"`
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin

  a = addfile("wrfchem1/WRFV3/test/em_real/wrfout_d02_2016-11-03_06:00:00.nc","r")
  
type = "png"
wks = gsn_open_wks(type,"cin")

  T  = a->T
  P  = a->P
  PB = a->PB
  QV = a->QVAPOR
  PH  = a->PH
  PHB = a->PHB
  HGT = a->HGT
  PSFC = a->PSFC

  T = T+300.
  P  = P+PB
  tk = wrf_tk(P,T)
  PH =  PH + PHB
  z = wrf_user_unstagger(PH,PH@stagger)
  z = z/9.81
; print(P)
  cinfo = wrf_cape_2d(P,tk,QV,z,HGT,PSFC,True)
  print(cinfo(1,:,:,:))
  cin   = cinfo(1,:,:,:)
   lat = a->XLAT(0,:,50)
   lat@unit = "degree_north"
   lon = a->XLONG(0,50,:)
   lon@unit = "degree_east"
   cin!1 = "lat"
   cin!2 = "lon"
   cin&lat = lat
   cin&lon = lon
  printVarSummary(cin)


  res               = True              ; Variable to hold plot options
  res@gsnDraw       = False
  res@gsnFrame      = False
  res@gsnAddCyclic  = False
  res@mpOutlineOn   = True
  res@cnFillOn      = True              ; Turn on contour fill.
  res@cnLinesOn     = False
  res@cnLineLabelsOn = False
;  res@cnLevelSelectionMode = "ManualLevels"       ; manually select levels
;  res@cnLevelSpacingF      = 1000                  ; contour spacing
;  res@cnMinLevelValF       = 95000                 ; min level
;  res@cnMaxLevelValF       = 103000                 ; max leve
;  res@lbOrientation = "Vertical"        ; Move labelbar to side.
  res@gsnMajorLatSpacing = 3
  res@gsnMajorLonSpacing = 3
  res@mpGeophysicalLineThicknessF = 3.
  res@mpProvincialLineThicknessF = 3.
  res@pmTickMarkDisplayMode = "Always"
  res@mpDataBaseVersion = "MediumRes"
  res@mpDataSetName     = "Earth..4"
  res@mpOutlineSpecifiers = (/"China:states"/)
  res@mpLandFillColor = 0
  res@lbLabelBarOn  = True
  res@mpLimitMode = "Latlon"
  res@mpMinLatF = 34
  res@mpMaxLatF = 43
  res@mpMinLonF = 112
  res@mpMaxLonF = 122
;  res@cnFillDrawOrder = "PreDraw"
  res@mpOceanFillColor = 0
  res@mpInlandWaterFillColor = 0
  res@gsnLeftString = ""
  res@gsnRightString = ""  

plot = gsn_csm_contour_map(wks, cin(0,:,:), res)
  draw(plot)
end

SG哭晕在厕所

自顶{:5_275:}

SG哭晕在厕所

真的没人吗。。。求救啊

南南偶卖糕Oo

请问楼主有解决吗哭唧唧