;*************************************************
; zonal_1.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" 
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
;************************************************
begin
;************************************************
; variable and file handling
;************************************************
  fn  = "83.nc"                                 ; define filename
  in  = addfile(fn,"r")                         ; open netcdf file
  ts = in->TS                                   ; select variable to ave
;************************************************
; data processing
;************************************************
  zave = dim_avg(ts)                            ; calculate zonal ave
;************************************************
; plotting
;************************************************
  wks  = gsn_open_wks("ps","zonal")                 ; Opens a ps file 
  
  plot=gsn_csm_xy(wks,zave(0,:),ts&lat,False)       ; plot zonal ave
;***********************************************
end