! GeInitRndm.f90 ! ! Initialize the standard f90 random generator. ! All elements of the seed array are initialized ! to a common seed different from 0. If the user ! variable rndmSeed equals 0 no seed gets loaded. ! ! Detlef Irmscher, Thomas Ullrich, Uni Heidelberg, September 1993 ! Last update: 17.3.95 tu ! subroutine GeInitRndm use genesis_m, self=>GeInitRndm implicit none integer :: isize integer, allocatable :: seed(:) if (rndmSeed /= 0) then call random_seed(size=isize) ! get size of seed array allocate(seed(isize)) seed = rndmSeed ! initialize seed array call random_seed(put=seed) deallocate(seed) ! free seed array endif end subroutine GeInitRndm