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/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin            
a = addfile("2020slp1.nc","r")      
b = addfile("2019w.nc","r")      
  wks   = gsn_open_wks ("eps", "slp")
  do it = 0,0,1
    msl1 = a->msl(it,:,:)     ; mean sea level pressure
    u101 = b->u10(it,:,:)     ; 10 metre U wind component
    v101 = b->v10(it,:,:)     ; 10 metre V wind component
    msl = short2flt(msl1)  
    msl = smth9(msl,0.50,0.25,False)
    u10 = short2flt(u101)
    v10 = short2flt(v101)
    msl = msl/100.
       msl@units="hPa"

;绘图参数控制
res                = True
res@gsnDraw        = False       ;为调整绘制顺序，让该图暂时不画出来
res@gsnFrame       = False
res@gsnLeftString = ""
res@gsnRightString = ""
  res@gsnAddCyclic = False
res@mpFillOn                = False
res@mpMinLatF               = 10.
res@mpMaxLatF               = 50.
res@mpMinLonF               = 75.
res@mpMaxLonF               = 135.
res@mpDataBaseVersion       = "MediumRes"
res@mpDataSetName           = "Earth..4"
res@mpOutlineSpecifiers     = (/"China:states"/)
res@pmTickMarkDisplayMode  = "Always"
res@tmXBMajorLengthF = 0.005

; Set some contouring resources.
res1   =  True
  res1 = res
  res1@cnSmoothingOn = True ;
  res1@cnSmoothingDistanceF   = 0.001

  res1@cnLevelSelectionMode  = "ManualLevels"
  res1@cnLevelSpacingF  = 2.5
  res1@cnInfoLabelOn = False
  res1@cnLineLabelsOn          =  True
  res1@cnLineLabelFontHeightF = 0.01
  res1@cnLineThicknessF = 2.0
  res1@cnLineLabelInterval = 1

  map = gsn_csm_contour_map(wks,msl, res1)
  draw(map)
  delete(res1)

;create vector plot
  res3 = True
  res3 = res
  res3@vcRefMagnitudeF= 20.0
  res3@vcRefLengthF= 0.05
  res3@vcMinDistanceF= 0.03
  res3@vcGlyphStyle= "CurlyVector"
  res3@vcRefAnnoOn = True   ;do not draw reference vector annotation
  res3@vcRefAnnoOrthogonalPosF = -0.15  ; vertical position
  res3@vcRefAnnoParallelPosF = 0.98
  res3@vcLineArrowThicknessF = 2.0
  res3@gsnDraw = False
  res3@gsnFrame = False
  res3@gsnLeftString = " "
  res3@gsnRightString = " "
  map2= gsn_csm_vector_map(wks,u10,v10,res3)
  draw(map2)
  delete(res3)
  frame(wks)
end do
end

1099221723

做一下平滑 smth9

赵狗剩

1099221723 发表于 2020-3-4 09:13
做一下平滑 smth9

现在已经是平滑后的结果了