;*************************************************
; zonal_3.ncl
;
; Concepts illustrated:
;   - Attaching a zonal means plot to a cylindrical equidistant contour plot
;************************************************
;
; These files are loaded by default in NCL V6.2.0 and newer
; 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
;************************************************
  in  = addfile("83.nc","r") 
  ts = in->TS                                   ; select variable to ave
;************************************************
; plotting
;************************************************
  wks  = gsn_open_wks("png","zonal")            ; send graphics to PNG file
  
  res                 = True		        ; make plot mods
  res@tiMainString    = "Zonal Average"         ; Title for the plot
  res@gsnZonalMean    = True                    ; put zonal on side
  res@gsnZonalMeanXMinF = 240.			; set minimum X-axis value for zonal mean plot  
  res@gsnZonalMeanXMaxF = 315.			; set maximum X-axis value for zonal mean plot  
  res@gsnZonalMeanYRefLine = 273.15		; set reference line X-axis value
  
  plot=gsn_csm_contour_map(wks,ts(0,:,:),res)   ; plot temp contours 
;***********************************************
end