da_set_randomcv.inc

1 subroutine da_set_randomcv(cv_size, rcv)
2 
3    !---------------------------------------------------------------------------
4    ! Purpose: Allocate components of control variable.
5    !
6    ! Method:  Allocate component in turn.
7    !---------------------------------------------------------------------------
8 
9    implicit none
10 
11    integer, intent(in) :: cv_size
12    real, intent(inout) :: rcv(1:cv_size)     ! Control variable structure.
13 
14    integer             :: i
15    real                :: mean_cv, rms_cv, std_dev_cv
16    real                :: sum_cv, sum_cv2
17 
18    ! [1] Initialize random number generator and scalars:
19 
20    call da_random_seed
21 
22    sum_cv = 0.0
23    sum_cv2 = 0.0
24 
25    ! [2] Calculate random numbers with Gaussian distribution:
26 
27    do i = 1, cv_size
28       call da_gauss_noise(rcv(i))
29 
30       sum_cv = sum_cv + rcv(i)
31       sum_cv2 = sum_cv2 + rcv(i) * rcv(i)
32    end do
33 
34    mean_cv = sum_cv / real(cv_size)
35    rms_cv = sqrt(sum_cv2 / real(cv_size))
36    std_dev_cv = sqrt(rms_cv * rms_cv - mean_cv * mean_cv)
37 
38    write(unit=stdout)
39    write(unit=stdout,fmt='(a)')' Gaussian (Normal) noise statistics:'
40    write(unit=stdout,fmt='(a,f15.5)')' Mean = ',mean_cv
41    write(unit=stdout,fmt='(a,f15.5)')' RMS = ', rms_cv
42    write(unit=stdout,fmt='(a,f15.5)')' STD DEV = ', std_dev_cv
43 
44 end subroutine da_set_randomcv
45 
46