;***********************************************
; CSM_Graphics: color_6.ncl
;***********************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" 

begin
;************************************************
; Read the file
;************************************************
  in    = addfile("80.nc","r")
  TS    = in->TS
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","color")              ; open a ps file
  gsn_define_colormap(wks,"gui_default")
  res                      = True               ; plot mods desired
  res@tiMainString         = "Explict Contour Levels"

  res@cnFillOn             = True               ; turn on color fill
  res@cnLinesOn            = False              ; turn off contour lines
  res@gsnSpreadColors      = True               ; use full colormap
  res@lbOrientation        = "vertical"         ; vertical label bars

  res@cnLevelSelectionMode = "ExplicitLevels"   ; set explicit contour levels
  res@cnLevels            = (/232.7,234.2,238,240,244,248.4,252,258.3,276,286.5,292.1,300,306/)

  res@mpFillOn             = False              ; turn off gray continents
  res@mpCenterLonF         = 180                ; Centers the plot at 180 

  plot = gsn_csm_contour_map_ce(wks,TS(0,:,:), res)   ; create plot
  
end