ncl在数据上叠加地形文件问题

chongzika

具体的脚本
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  filename = "2001-2010.nc"
  a        = addfile(filename,"r")
mbl  = a->MBL(0,:,:)   ; read in example data [2D only here]
  copy_VarCoords(dstsfmbl,mbl)
  type = "ps"    ;  "pdf" , "ncgm"
  wks = gsn_open_wks(type, "china")               ; open a  file
  gsn_define_colormap(wks,"MPL_YlOrBr")  ; choose colormap
; Set some basic resources
  res = True
  res@cnFillOn             = True                  ; turn on color fill
  res@mpLandFillColor        = "white"
  res@gsnSpreadColors      = True                  ; use full range of colors
  res@gsnRightString        = "t/km~S~2~N~";
  res@gsnLeftString        ="mbl"
; =================================================;
  res@mpMinLatF              =  17
  res@mpMaxLatF              =  55
  res@mpMinLonF              =  72
  res@mpMaxLonF              =  136
; =================================================;
  res@pmTickMarkDisplayMode = "Always"              ;
res@cnLevelSelectionMode = "ExplicitLevels"   ; set explicit contour levels
res@cnLevels            = (/10,20,50,100,200,500,1000/)
plot = gsn_csm_contour_map_ce（wks,mbl,res)   ; Create contour plot

;>============================================================<
;                      add China map
;>------------------------------------------------------------<
;---Attach the shapefile polylines using files read off gadm.org/country.
  china_shp_name=/data/map/bou1_4l.shp"
  province_shp_name=/data/map/bou2_4p.shp"
  river1_shp_name="/data/map/hyd1_4l.shp"

  lnres                  = True
  lnres@gsLineColor      = "gray25"
  lnres@gsLineThicknessF = 0.5   

  chi_id = gsn_add_shapefile_polylines(wks,plot,china_shp_name,lnres)
  pro_id = gsn_add_shapefile_polylines(wks,plot,province_shp_name,lnres)
  ri1_id = gsn_add_shapefile_polylines(wks,plot,river1_shp_name,lnres)

  draw(plot)       ; This will draw the map and the shapefile outlines.
  frame(wks)       ; Advance the frame

end
先画了一张没有地形文件的图