function hrms,f,xvertical=xvertical,yvertical=yvertical
;
;  Horizontal rms of 3-d scalar
;
	s=size(f)
	if n_elements(xvertical) ne 0 then begin
	  h=fltarr(s(1))
	  for n=0,s(1)-1 do h(n)=rms(f(n,*,*))
	endif else if n_elements(yvertical) ne 0 then begin
	  h=fltarr(s(2))
	  for n=0,s(2)-1 do h(n)=rms(f(*,n,*))
	end else begin
	  h=fltarr(s(3))
	  for n=0,s(3)-1 do h(n)=rms(f(*,*,n))
	end
	return,h
end