init_random.F File Reference

#include "implicit_f.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	init_random ()

Function/Subroutine Documentation

init_random()

subroutine init_random

Definition at line 30 of file init_random.F.

!$COMMENT
!       init_random description
!       initialization of RAND_STRUCT structure     
!       
!       INIT_RANDOM organization :
!       - check -rxalea or -rseed cmd line is used
!       - check the value or -rxalea / -rseed and print
!         a warning if the value is not good
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE random_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        IF(rand_struct%CMD) THEN
 
            WRITE(iout,1000)
 
            IF(rand_struct%ALEA.AND.(.NOT.rand_struct%SEED) ) THEN
                !   apply default value to seed
                rand_struct%SEED=.true.
                rand_struct%SEED_NBR=0.d+00
                WRITE(iout,'(A)') ' '
                WRITE(iout,'(A)') '** WARNING : -RSEED OPTIONAL OPTION IS NOT USED'
                WRITE(iout,'(A)') '              DEFAULT VALUE WILL BE APPLIED'
                WRITE(iout,'(A)') ' '
            ENDIF
            IF(rand_struct%SEED) THEN
                IF(rand_struct%SEED_NBR<0.d+00.OR.rand_struct%SEED_NBR>1.d+00) THEN
                    !   error, 0 <= seed < 1
                    WRITE(iout,'(A)') ' '
                    WRITE(iout,'(A)') '** WARNING : -RSEED VALUE MUST BE DEFINED BETWEEN'
                    WRITE(iout,'(a)') '              0. and 1.'     
                    WRITE(IOUT,'(a)') ' '
                ENDIF
            ENDIF
            IF(RAND_STRUCT%ALEA) THEN
.OR.                IF(RAND_STRUCT%ALEA_NBR<0.d+00RAND_STRUCT%ALEA_NBR>10.d+30) THEN
                    !   error, 0 <= alea < +infinity
                    WRITE(IOUT,'(a)') ' '
                    WRITE(IOUT,'(a)') '** WARNING : -RALEA VALUE MUST BE DEFINED BETWEEN'
                    WRITE(iout,'(A)') '              0 AND +INFINITY'       
                    WRITE(iout,'(A)') ' '
                ENDIF
            ELSE
                rand_struct%CMD=.false.
                WRITE(iout,'(A)') ' '
                WRITE(iout,'(A)') '** WARNING : -RSEED OPTIONAL OPTION MUST BE USED WITH -RXALEA OPTION'
                WRITE(iout,'(A)') '             -rseed option will be ignored'
                WRITE(IOUT,'(a)') ' '
            ENDIF
        ENDIF
 
 1000 FORMAT(//
     .'     random noise       '/
     .'     ------------       ')
 
        RETURN