求助。实习NCL制图显示如下

hhtmyc

好不容易做到这里，但是老师没具体讲画图。我们用的NCL也是直接拷贝的，这个NCL是从老师给的模板改的，但是我其实好多参数还没有弄懂。希望大神能有针对性的帮我分析一下。谢谢了~~


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/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "$NCL_FUNC/functions.ncl"

begin
;

   dirdta = "/opt/WRFV3/test/em_real/"

   nc_file = addfile(dirdta+"wrfout_d01_2009-08-09_12 00 00.nc","r")

   figname = "FILE"

   lats = wrf_user_getvar(nc_file,"lat",0)
   lons = wrf_user_getvar(nc_file,"lon",0)
   nn = dimsizes(lats)
   nx = nn(1)
   ny = nn(0)

   mdbz = wrf_user_getvar(nc_file,"mdbz",0)
   pre  = wrf_user_getvar(nc_file,"pressure",0)  ;hPa
   ua   = wrf_user_getvar(nc_file,"ua",0)
   va   = wrf_user_getvar(nc_file,"va",0)

   plev = 850.0
   uwnd = wrf_user_intrp3d(ua,pre,"h",plev,0,False)
   vwnd = wrf_user_intrp3d(va,pre,"h",plev,0,False)

   xstride = 10
   ystride = 10
   uwnds = new(dimsizes(uwnd),float)
   vwnds = new(dimsizes(uwnd),float)
   uwnds(::ystride,::xstride) = uwnd(::ystride,::xstride)
   vwnds(::ystride,::xstride) = vwnd(::ystride,::xstride)


   wks = gsn_open_wks("pdf",figname)

   cmap = (/(/1.0,   1.0,    1.0   /), (/   0.0, 0.0,    0.0/), \
           (/1.0,    1.0,    1.0   /), (/0.0039, 0.6172, 0.9570/), \
           (/0.0,    0.9219, 0.9219/), (/   0.0, 0.8438, 0.0/), \
           (/0.0039, 0.5625, 0.0   /), (/0.9922, 0.9922, 0.0/), \
           (/0.8984, 0.7500, 0.0   /), (/0.9922, 0.5625, 0.0/), \
           (/0.9922, 0.0,    0.0   /), (/0.8320, 0.0,    0.0/), \
           (/0.7422, 0.0,    0.0   /), (/0.9922, 0.0,    0.9375/), \
           (/0.5938, 0.0,    0.7   /), (/0.6758, 0.5625, 0.9375/)/)
   gsn_define_colormap(wks, cmap)

   mpres = True
;  mpres@gsnFrame = False
;  mpres@gsnDraw  = False
;   mpres@mpProjection = "LambertConformal"
;  mpres@mpLambertParallel1F = nc_file@TRUELAT1
;  mpres@mpLambertParallel2F = nc_file@TRUELAT2
;  mpres@mpLambertMeridianF  = nc_file@STAND_LON
;  mpres@mpLimitMode = "LatLon"
;  mpres@sfXArray  = lons
;  mpres@sfYArray  = lats
;  mpres@vfXArray  = lons
;  mpres@vfYArray  = lats
;  mpres@mpMinLatF = max(lats(0,:))
;  mpres@mpMaxLatF = min(lats(ny-1,:))
;  mpres@mpMinLonF = max(lons(:,0))
;  mpres@mpMaxLonF = min(lons(:,nx-1))
  mpres@mpDataBaseVersion="Ncarg4_1"
  mpres@mpDataSetName="Earth..4"
  mpres@mpOutlineSpecifiers=(/"China:states","Taiwan"/)
;  mpres@mpPerimOn = True
;  mpres@mpGridAndLimbOn      = False
;  mpres@pmTickMarkDisplayMode= "Always"
   mpres@tmXTOn = "False"
   mpres@tmYROn = "False"

   opts = mpres
   opts@cnFillOn = True
   opts@cnLinesOn = False
   opts@cnLineLabelsOn= False
   opts@cnInfoLabelOn = False
   opts@cnSpanFillPalette = True
   opts@cnLevelSelectionMode = "ExplicitLevels"
   opts@cnLevels = ispan(5,75,5)
   opts@lbLabelBarOn = True
   opts@gsnSpreadColorStart = 2        ; start at 2 (the default)
   opts@gsnSpreadColorEnd   = 15
   opts@lbTitleString = "~Z75~Composite reflectivity (dBZ)"
   opts@tiMainString = "~Z85~U-V wind at "+plev+"hPa"
;   ctr_ref = gsn_csm_contour_map(wks,mdbz,opts)
   ctr_ref = wrf_contour(nc_file,wks,mdbz,opts)
   delete(opts)

   opts = mpres
   opts@vcGlyphStyle = "LineArrow"
   opts@vcLineArrowThicknessF = 1.2
   opts@vcRefMagnitudeF = 30.0
   opts@vcRefLengthF = 0.04
   opts@vcMinFracLengthF = 0.3
   opts@vcRefAnnoZone = 0
   opts@vcRefAnnoOrthogonalPosF = -0.5
   opts@vcRefAnnoParallelPosF = 0.5
;   vtr_wnd = gsn_csm_vector(wks,uwnds,vwnds,opts)
   vtr_wnd = wrf_vector(nc_file,wks,uwnds,vwnds,opts)
   delete(opts)

;   overlay(ctr_ref,vtr_wnd)
   plot = wrf_map_overlays(nc_file,wks,(/ctr_ref,vtr_wnd/),True,mpres)


end

霹雳电光鼠

都有提示自己翻译一下就能看懂

sinagzy

同意楼上+1