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

begin
;
; Note that "t" on the file has no attributes like "long_name"
; or "units".
;
  tf = addfile("Tstorm.cdf","r")               ; Open netCDF file.
  T  = tf->t(0,:,:)                            ; Read in temperature data.

;  printVarSummary(T)

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

  res = True                          ; No plot options set.

  plot = gsn_csm_contour(wks,T,res)   ; Call the gsn_csm function for
                                      ; drawing the contours.
end