Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(in) | :: | mean | |||
real(kind=rp), | intent(in) | :: | sigma |
FUNCTION random_norm(mean,sigma)
REAL(rp), INTENT(IN) :: mean
REAL(rp), INTENT(IN) :: sigma
REAL(rp) :: random_norm
REAL(rp) :: rand1, rand2
call RANDOM_NUMBER(rand1)
call RANDOM_NUMBER(rand2)
random_norm = SQRT(-2.0_rp*LOG(1.0_rp-rand1))*COS(2.0_rp*C_PI*rand2);
END FUNCTION random_norm