ncl怎么只画陆地的部分

chongzika

看了看官网MASKING的解释，;*************************************************; mask_1.ncl;; Concepts illustrated:;   - Using "mask" to set land or ocean values in your data to missing;   - Masking the ocean in a map plot;   - Masking land in a map plot;   - Spanning part of a color map for contour fill;;
************************************************l
oad "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
in = addfile("atmos.nc","r")  
ts = in->TS(0,:,:)
oro = in->ORO(0,:,:)
;************************************************
; use mask function to mask out land and then ocean data; ocean=0,land=1,sea_ice=2
;************************************************
land_only  = ts                              ; trick to keep cv's and atts
ocean_only = ts  
land_only  = mask(ts,oro,1)   ; Return the values of ts where oro=1  
ocean_only = mask(ts,oro,0)   ; Return the values of ts where oro=0
;************************************************; common resources
;************************************************
wks = gsn_open_wks("ps","mask")                  ; open a ps file  
gsn_define_colormap(wks,"BlAqGrYeOrRe")           ; choose colormap
res                      = True     ; plot mods desired
res@cnFillOn             = True     ; turn on color fill
res@cnLinesOn            = False    ; turn of contour lines  
res@gsnSpreadColors      = True     ; use full colormap
res@gsnSpreadColorStart  = 10       ; start at color 10
res@gsnSpreadColorEnd    = 96       ; end at color 96
;************************************************
; individual plots
;************************************************
res@cnLevelSpacingF         = 3                ; interval
  res@lbLabelStride           = 4                ; every 4th label  
res@tiMainString            = "Land Only"      ; title  
plot = gsn_csm_contour_map_ce(wks,land_only,res)
res@cnLevelSpacingF         = 2                ; interval
res@lbLabelStride           = 2                ; every 2nd label  
res@tiMainString            = "Ocean Only"     ; title
plot = gsn_csm_contour_map_ce(wks,ocean_only,res)
end
我的数据里面没有这TS，ORO怎么画呢

longlivehj

masking里面有很多例子，你提到的只是其中一种方式，看看其它例子，总有一个适合你！