osse/EnKF/ranvar.F

        subroutine ranvar(xens,ave,n,nens,var)

        implicit none

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

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