var/da/da_tools/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    if (trace_use) call da_trace_entry("da_set_randomcv")
  19
  20    ! [1] Initialize random number generator and scalars:
  21
  22    !moved to da_solve.inc
  23    !call da_random_seed
  24
  25    sum_cv = 0.0
  26    sum_cv2 = 0.0
  27
  28    ! [2] Calculate random numbers with Gaussian distribution:
  29
  30    do i = 1, cv_size
  31       call da_gauss_noise(rcv(i))
  32
  33       sum_cv = sum_cv + rcv(i)
  34       sum_cv2 = sum_cv2 + rcv(i) * rcv(i)
  35    end do
  36
  37    mean_cv = sum_cv / real(cv_size)
  38    rms_cv = sqrt(sum_cv2 / real(cv_size))
  39    std_dev_cv = sqrt(rms_cv * rms_cv - mean_cv * mean_cv)
  40
  41    write(unit=stdout,fmt=*)
  42    write(unit=stdout,fmt='(a)')' Gaussian (Normal) noise statistics:'
  43    write(unit=stdout,fmt='(a,f15.5)')' Mean    = ', mean_cv
  44    write(unit=stdout,fmt='(a,f15.5)')' RMS     = ', rms_cv
  45    write(unit=stdout,fmt='(a,f15.5)')' STD DEV = ', std_dev_cv
  46
  47    if (trace_use) call da_trace_exit("da_set_randomcv")
  48
  49 end subroutine da_set_randomcv
  50
  51