load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  tf = addfile("meccatemp.cdf","r")
  T  = tf->t(0,:,:)                      ; T will have metadata attached.

  printVarSummary(T)                     ; T has coord arrays and attributes
  printVarSummary(T&lon)                 ; T&lon has units attribute

  wks = gsn_open_wks("ps","contour2a")   ; Open PS file called "contour2a.ps"

  res = True                             ; No plot options set.

  plot = gsn_csm_contour_map(wks,T,res)  ; Call the gsn_csm function for
                                         ; drawing contours over a map.
end