subroutine ranmean2(xens,ave,n,nens)

	implicit none

	real*8, intent(in)  :: xens(n,nens)
	real*8, intent(out) :: ave(n)
	real*8 rnens
	integer i, nens, n

	rnens=nens
	ave(:)=xens(:,1)
	do i=2,nens
	 ave(:)=ave(:)+xens(:,i)
	enddo
	ave=(1./(rnens))*ave
	end subroutine ranmean2