;
; Contours over a map exercise 3.
;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  diri = "$NCARG_ROOT/lib/ncarg/data/cdf/"
  tf = addfile(diri+"meccatemp.cdf","r")       ; Open netCDF file.
  t  = tf->t(0,:,:)                            ; Read in temperature data.

  wks = gsn_open_wks("x11","contourmap_ex03")  ; Open X11 window

  res = True                                  ; Plot options desired.

  res@gsnMaximize           = True            ; Maximize plot in frame.
  res@gsnAddCyclic          = False           ; Don't add a cyclic point.

  res@cnLevelSelectionMode = "ManualLevels"   ; Manually set contour levels
  res@cnMinLevelValF       =  195.0
  res@cnMaxLevelValF       =  328.0
  res@cnLevelSpacingF      =    5.0

  res@cnFillOn             = True            ; Turn on contour fill.
  res@cnLinesOn            = False           ; Turn off contour lines.
  res@cnFillPalette        = "rainbow"       ; Change the color map 

  plot = gsn_csm_contour_map(wks,t,res)
end