osse/EnKF/random2.F

c***  returns a vector of random values N(variance=1,mean=0)
        subroutine random2(work1,n)

        implicit none

        integer, intent(in) :: n
        real,   intent(out) :: work1(n)
        real,   allocatable :: work2(:)
        real,   pi

        allocate (work2(n))

        pi=3.14159
        call random_number(work1)
        call random_number(work2)
        work1=sqrt(-2.*log(work1))*cos(2.*pi*work2)

        deallocate(work2)
        end subroutine random2
1	c*** returns a vector of random values N(variance=1,mean=0)
2	subroutine random2(work1,n)
3
4	implicit none
5
6	integer, intent(in) :: n
7	real, intent(out) :: work1(n)
8	real, allocatable :: work2(:)
9	real, pi
10
11	allocate (work2(n))
12
13	pi=3.14159
14	call random_number(work1)
15	call random_number(work2)
16	work1=sqrt(-2.log(work1))cos(2.piwork2)
17
18	deallocate(work2)
19	end subroutine random2