;*************************************************
; zonal_2.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/contributed.ncl" 
;************************************************
begin
;************************************************
; variable and file handling
;************************************************
  in  = addfile("83.nc","r")
  ts  = in->T                                        ; select variable to ave
; use wrapper function so that meta data is preserved
  zave = dim_avg_Wrap(ts)                            ; calculate zonal ave
;************************************************
; plotting
;************************************************
  wks  = gsn_open_wks("ps","zonal")                ; Opens a ps file 

  res            = True                            ; plot mods desired
  res@trYReverse = True                            ; reverse Y axis
  res@trXReverse = True                            ; reverse X axis

  plot=gsn_csm_contour(wks,zave(0,:,:),res)             ; plot zonal ave
;***********************************************
end