function haver,f,yvertical=yvertical,array=array,gz=gz
;
;  Horizontal average of 3-d scalar
;
	COMMON cvert, cxvert
	s=size(f)
	if s(0) eq 2 then begin
	  h=total(f,2)/s(2)
	  if keyword_set(array) then $
	    h=rebin(reform(h,s(1),1),s(1),s(2))
	end else if n_elements(xvertical) ne 0 or keyword_set(cxvert) then begin
	  h=total(total(f,1),1)/(s(1)*s(2))
	  if keyword_set(array) then $
	    h=rebin(reform(h,1,1,s(3)),s(1),s(2),s(3))
	end else begin
	  h=total(total(f,3),1)/(s(1)*s(3))
	  if keyword_set(array) then $
	    h=rebin(reform(h,1,s(2),1),s(1),s(2),s(3))
	end
	if keyword_set(gz) then h=h(4:n_elements(h)-5)
	return,h
end