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	afe	1.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