如何将地图外等值线去除

英英茵

想要将地图外多余的等值线去除  试了res@cnLineDrawOrder = "PreDraw"  但是还是不行  有没有大神能够指点一下

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.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/wrf/WRF_contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"

begin
fname = "/home/WuQiuuu/quyu/esccr132.txt"
lines = asciiread(fname,-1,"string")

;delim = str_get_tab()
lon = stringtofloat(str_get_field(lines(1:),2," "))
lat = stringtofloat(str_get_field(lines(1:),3," "))
prec = stringtofloat(str_get_field(lines(1:),4," "))
prex = stringtofloat(str_get_field(lines(1:),5," "))
preq = stringtofloat(str_get_field(lines(1:),6," "))
pred = stringtofloat(str_get_field(lines(1:),7," "))


olon = fspan(95,115,100)
olat = fspan(20,35,30)


olon!0          = "lon"
olon@long_name  = "lon"
olon@units      = "degrees-east"
olon&lon        = olon
olat!0          = "lat"
olat@long_name  = "lat"
olat@units      = "degrees_north"
olat&lat        = olat

prec@_FillValu4e = -999.00
rscan = (/5,3,1/)  ;连续的有效半径大小，最大为10，依次递减
;data1 = obj_anal_ic_deprecated(lon,lat,prec,olon,olat,rscan, False)  ;Creanm插值
data1 = obj_anal_ic_Wrap(lon,lat,prec,olon,olat,rscan, False)  ;Creanm插值


wks = gsn_open_wks("png","/home/WuQiuuu/quyu/prechunxishu")

;-----------------map
mpres= True

mpres@gsnDraw   = False   ; don't draw yet
mpres@gsnFrame  = False
mpres@gsnMaximize   = True
;mpres@gsnAddCyclic        = False
mpres@mpMinLonF            = 96.               ; select a subregion
mpres@mpMaxLonF            = 113.
mpres@mpMinLatF            = 20.
mpres@mpMaxLatF            = 35.

mpres@mpPerimOn   = True
mpres@mpOutlineOn = True
mpres@mpFillOn    = False
mpres@vpWidthF    = 0.68
mpres@vpHeightF   = 0.6
mpres@mpDataBaseVersion = "Ncarg4_1"
mpres@mpDataSetName = "Earth..4"
mpres@mpOutlineBoundarySets     = "NoBoundaries"     
;mpres@mpOutlineSpecifiers      = (/"Continents","China"/)  ; continents only
mpres@mpNationalLineThicknessF= 2
mpres@mpGeophysicalLineThicknessF = 2
mpres@mpNationalLineColor    = "gray52"

map = gsn_csm_map(wks, mpres)

ShpDir ="D:/CYGwin/usr/local/lib/ncarg/nclscripts/cnmap/xinan4.shp"
resmp = True
resmp@gsLineThicknessF= 2
resmp@gsLineColor    = "black"

resmp@mpAreaMaskingOn = True
resmp@mpMaskAreaSpecifiers = (/"xinan4.shp"/)
resmp@mpLandFillColor = 0
resmp@mpOceanFillColor=0
resmp@mpOutlineOn=True
resmp@mpOutlineSpecifiers="xinan4.shp"

poly = gsn_add_shapefile_polylines(wks,map,ShpDir,resmp)

res = True

res@vpWidthF    = 0.68
res@vpHeightF   = 0.6
res@gsnAddCyclic        = False
res@gsnSpreadColors     = True     ; spread out color table
res@cnFillOn = False                ; turn on color
res@cnLinesOn = True
res@cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
res@cnMinLevelValF       = -1              ; set min contour level
res@cnMaxLevelValF       =  1              ; set max contour level
res@cnLevelSpacingF      =  0.05              ; set contour spacing
res@cnLineDrawOrder = "PreDraw"

res@sfXArray    = olon
res@sfYArray    = olat

res@cnLineLabelsOn       = True
res@cnLabelDrawOrder = "PreDraw"
res@gsnContourZeroLineThicknessF = 2. res@gsnContourNegLineDashPattern = 1         ; sets negative contours to dash pattern 1

plot=gsn_csm_contour(wks,data1,res)

overlay(map, plot)

maximize_output(wks,False)        


end

月满西楼

在地图部分加上这三句试试：
mpres@mpLandFillColor                = "white"
mpres@mpInlandWaterFillColor     = "white"
mpres@mpOceanFillColor              = "white"

英英茵

月满西楼 发表于 2018-12-17 09:27
在地图部分加上这三句试试：
mpres@mpLandFillColor                = "white"
mpres@mpInlandWaterFillC ...

设置了的  并没有用