NCL绘制wrfout中qcloud的垂直剖面出错

小其其格

出现的错误是：fatal:Dimension (Vertical) of (qcloud) does not have an associated coordinate variable

代码如下：

begin
;---Open WRF output file
  filename = "wrfout_d03_2013-07-29_09_00_00"
  a        = addfile(filename,"r")

;---Read temperature at first time step
  time=0
  qcl = wrf_user_getvar(a,"QCLOUD",time)
  z  = wrf_user_getvar(a, "height",time)        ; grid point height
  p  = wrf_user_getvar(a, "pressure",time) ; total pressure
  dimsrh = dimsizes(qcl)
  plane = (/ dimsrh(2)/2, dimsrh(1)/2 /)    ; pivot point is center of domain

; Extract data on a W-E cross section that runs through "plane"
; Interpolate to height
  qcloud = wrf_user_intrp3d(qcl,z,"v",plane,90.0,False)

  wks = gsn_open_wks("pdf","wrf_gsn")

;---Resources for filled contour plot

  res                      = True                 ; plot mods desired

  res@cnLevelSelectionMode = "ManualLevels"       ; manually select levels
  res@cnLevelSpacingF      = 5.0                  ; contour spacing
  res@cnMinLevelValF       = -50.                 ; min level
  res@cnMaxLevelValF       =  50.                 ; max level
  res@cnLineLabelsOn       = True                 ; turn on line labels
  res@cnFillOn             = True                 ; turn on color fill
  res@cnFillPalette        = "BlWhRe"             ; choose colormap
  res@lbOrientation        = "Vertical"
  res@tiMainString         = filename

;---Necessary for contours to be overlaid correctly on WRF projection
  ;res@tfDoNDCOverlay       = True
  ;res@gsnAddCyclic         = False

;---Set map resources based on projection on WRF output file
  res = wrf_map_resources(a,res)

;---Call plotting routine which draws contours over a map
  ;nl   = 5    ; bottomost level
  ;plot = gsn_csm_contour_map(wks,qcloud(nl,:,:),res)
  plot = gsn_csm_pres_hgt(wks,qcloud({29.8},:),res)
end



请问大神们错在哪里了？

阳光不怕晒

你的plot = gsn_csm_pres_hgt(wks,qcloud({29.8},:),res),其中qcloud不是应该是三维变量吗，你少了一维垂直的

小其其格

阳光不怕晒 发表于 2017-10-30 11:57
你的plot = gsn_csm_pres_hgt(wks,qcloud({29.8},:),res),其中qcloud不是应该是三维变量吗，你少了一维垂直 ...

好的！我再研究研究，谢谢回复