ncl 画青藏高原边界内图问题

雾撒微晨

begin

;************************************************
;************************************************
; create pointer to file and read in data
;************************************************
shp_fname = "/cygdrive/c/Users/Administrator/Desktop/NN/tibetan/tibet.shp"

fils = systemfunc("ls /cygdrive/c/Users/Administrator/Desktop/NN/GLDAS_CLM10SUBP_3H.A2010001.0000.001.grb")
printVarSummary(fils)   ; full path to files
f = addfiles(fils, "r")
ListSetType(f,"join")
printVarSummary(f)

; ListSetType(f,"join")
sw = f[:]->SOIL_M_GDS0_DBLY
printVarSummary(sw)
printVarSummary(swy)
copy_VarCoords(sw(0,:,:), swy)
udims = dimsizes(swy)
     opt             = True
     opt@return_mask = True
     TP_mask        = shapefile_mask_data(swy,shp_fname,opt)
   res                       = True   
   wks = gsn_open_wks("png", "sm11")
   gsn_define_colormap(wks,"MPL_RdYlBu")     

   ;res@gsnMaximize           = True   
   res@gsnDraw               = False
   res@gsnFrame              = False
   res@mpOutlineOn           = False
   res@cnFillOn              = True
   res@cnLinesOn             = False  
   res@cnLineLabelsOn        = False   
   res@lbLabelBarOn          = True   
   res@lbBoxLinesOn          = True
   res@cnSmoothingOn           = True     
   ;res@cnSmoothingDistanceF  =   True
   res@cnSmoothingTensionF   =   True
   res@gsnSpreadColors     = True      
   res@gsnSpreadColorStart  = 3            
  res@gsnSpreadColorEnd   = 9               
  res@cnLevelSelectionMode    = "ExplicitLevels"

   res@gsnAddCyclic        = False
  res@cnFillDrawOrder     = "Draw"
  ;axis and ticks set
   res@mpMinLonF            = 70               ; set min lon
   res@mpMaxLonF            = 110               ; set max lon
   res@mpMinLatF            = 24               ; set min lat
   res@mpMaxLatF            = 45               ; set max lat
      res@tmXBValues           = (/70.,80.,90.,100.,110./);,120.,130.,140/)
   res@tmXBLabels           = (/"70E","80E","90E","100E","110E"/);,"120E","130E","140E"/)
   res@tmYLValues           = (/20.,25.,30.,35.,40.,45./)        ;,50.,55/)     ;(/15.,
   res@tmYLLabels           = (/"20N","25N","30N","35N","40N","45N"/)       ;,"50N","55N"/)  ;(/"15N",  
   res@tmXTOn               =False           ;turns off the top tick marks
   res@tmYLMinorOn          =False           ;Turns off minor tick marks

  res@mpDataBaseVersion     = "MediumRes"
   res@mpLandFillColor         = "white"
   res@mpInlandWaterFillColor  = "white"
   res@mpOceanFillColor        = "white"
   res@gsnLeftString        ="soil moisture"
   res@gsnRightString        ="m~S~3~N~/m~S~3~N~"

plots1 = gsn_csm_contour_map(wks,swy,res)

   lnres  = True
   poly_mask0 = gsn_add_shapefile_polylines(wks,plots1,shp_fname,lnres)

   draw(plots1)
   frame(wks)

end