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/"
  fili = "Tstorm.cdf"
  tf   = addfile(diri+fili,"r")                ; Open netCDF file.
  t    = (tf->t(0,:,:)  - 273.15)*9./5. + 32   ; Read in temp, 
                                               ; convert to C.

  print(tf)                               ; Print what's on the file.
  printVarSummary(t)                      ; Print information about "t".

  wks = gsn_open_wks("x11","contourmap2_ex")   ; Open X11 window
 
  res            = True                     ; Plot options desired.
  res@gsnMaximize = True                    ; Maximize plot in frame.

  plot = gsn_csm_contour(wks,t,res)
end