ncl绘制wrfout数据错位

zyj1210

zyj1210

而且只在绘制wrfout数据时出现这样的问题，绘制其他数据可以随意变换绘图范围并且正确绘图

zyj1210

elev@description = "Terrain Height"
  elev@units = "m"
  wks = gsn_open_wks("png","dem")

;---Set some basic plot options
  ;mkres = True
  ;mkres@gsMarkerIndex = 16 ;标记的种类，16号为填色圆圈
  ;mk2res@gsMarkerSizeF = 5 ;设定标记的大小
  ;mk2res@gsMarkerColor = "red" ;设定标记颜色为colormap的第231个颜色

  res                   = True
  res@pmTickMarkDisplayMode = "Always"
  res@lbOrientation        = "Vertical"

  res@tiMainString         = "Southeast"
  ;res@mpOutlineSpecifiers = outline_areas
  res@gsnMaximize       = True   ; maximize plot in frame
  res@cnFillOn          = True
  res@cnLineLabelsOn     = False    ; turn off line labels
  res@cnInfoLabelOn      = False    ; turn off info label
  res@lbBoxLinesOn       = False    ; turn off labelbar box lines
  res@cnLevelSelectionMode = "ExplicitLevels"
  res@cnFillPalette        = "OceanLakeLandSnow"
  res@mpMinLatF            = 21         ; pamir limits
  res@mpMaxLatF            = 35
  res@mpMinLonF            = 97
  res@mpMaxLonF            = 113
  res@gsnAddCyclic        = False         ;
  res@mpOutlineOn         = True          ;
  res@mpGeophysicalLineColor = "black"    ;
  res@mpGeophysicalLineThicknessF = 2     ;
  res@mpOutlineBoundarySets  = "National" ;
  res@mpNationalLineColor = "black"        ;
  res@mpNationalLineThicknessF = 2.       ;
  res@mpDataBaseVersion = "MediumRes"
  res@mpOutlineOn       = True
  res@mpMaskAreaSpecifiers        = (/"Land"/)
  res@mpOutlineSpecifiers = outline_areas
  res@mpFillDrawOrder            = "PostDraw"
  res@cnConstFEnableFill  = True          ;
  ;res@cnLabelBarEndLabelsOn = True        ;
  ;res@cnConstFLabelOn     =  False        ;
  ;res@cnLabelBarEndStyle  = "ExcludeOuterBoxes" ;
  ;res@gsnSpreadColorStart = 2             ;
  ;res@gsnSpreadColorEnd   = 10            ;
  ;res@cnSpanFillPalette = True
  res@lbLabelAutoStride   = True          ;
  res@lbLabelStride       = 2             ;
  res@mpProvincialLineThicknessF    = 2   ;
  res@mpProvincialLineColor         = 1   ;
  res@tmYLLabelFontHeightF          = 0.02    ;set the size of y-Axis w
ords  ;res = wrf_map_resources(a,res)
;---Generate contour plot then overlay on map
  ;gsn_define_colormap(wks,"OceanLakeLandSnow");
  plot = gsn_csm_contour_map(wks,elev,res)
以上是代码主要部分

zyj1210

已解决：res@tfDoNDCOverlay      = False          ; Tell NCL you are doing a
native plot 这一句要设置为False