;*************************************************
; mask_2.ncl
;************************************************
load "$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,:,:)
;************************************************
; 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@lbLabelStride        = 4        ; every 4th label

  res@gsnSpreadColors      = True     ; use full colormap
  res@gsnSpreadColorStart  = 10       ; start at color 10
  res@gsnSpreadColorEnd    = 96       ; end at color 96

  res@cnLevelSpacingF      = 3        ; interval

  res@tiMainString         = "Draw land ON TOP of contours"      ; title

  res@cnFillDrawOrder      = "Predraw"  ; draw contours first
  plot = gsn_csm_contour_map_ce(wks,ts,res)

end