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	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