pro avinput,data,index=index,origin=origin,xxx,out,rrrr,full,$
	xy=xy,new=new,fit_param=fit_param,yfit=yfit,$
	noplot=noplot,multiple=multiple,unc_data=unc_data,$
	unc_sig=unc_sig,unc_fullsig=unc_fullsig,skip=skip
;+
; NAME
;  avinput  
;
; PURPOSE
;  To calculate the average intensity in a radial direction, for 
;  a three dimensional data cube.
;
; CALLING SEQUENCE
;  avinput,data,index=index,origin=origin,xxx,out,rrrr,full,xy=xy
; INPUT
;       
; OUTPUT
;
; OPTIONAL OUTPUT
;	
; OPTIONAL KEYWORD INPUT
;   	
; RESTRICTION
;   none
; HISTORY
;   Written CRF 12 Nov 95 : caf@mssl.ucl.ac.uk
;-

aalpha=n_elements(data(1,1,*))

azimuth_avg,data(*,*,0),xxx,index=index,out1,rrrr,full1,xy=xy,/noplot,origin=origin
pp=n_elements(out1)
pp1=n_elements(full1)
full=dblarr(aalpha,pp1)
out=dblarr(aalpha,pp)
out(0,*)=out1
full(0,*)=full1

for asa=1,aalpha-1. do begin
azimuth_avg,data(*,*,asa),xxx,out1,rrrr,index=index,full1,xy=xy,/noplot,origin=origin
out(asa,*)=out1
help,full,full1
full(asa,*)=full1
endfor

end