osse/EnKF/ranmean2.F

        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

1	subroutine ranmean2(xens,ave,n,nens)
2
3	implicit none
4
5	real*8, intent(in) :: xens(n,nens)
6	real*8, intent(out) :: ave(n)
7	real*8 rnens
8	integer i, nens, n
9
10	rnens=nens
11	ave(:)=xens(:,1)
12	do i=2,nens
13	ave(:)=ave(:)+xens(:,i)
14	enddo
15	ave=(1./(rnens))*ave
16	end subroutine ranmean2
17