!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  
subroutine RandomNumGen()
        use simulation
!implicit none
!       integer :: i,j, N, use_last
        integer, dimension (8) :: TimeForRand
        integer, dimension (2) :: rand_seed
        real ( kind = 8 )  :: x
        call DATE_AND_TIME(values=TimeForRand)     ! Get the current time
        rand_seed(1) = TimeForRand(4) * (360000*TimeForRand(5) + 6000*TimeForRand(6) + 100*TimeForRand(7) + TimeForRand(8))
        call RANDOM_SEED(PUT=rand_seed)          
        call random_number(harvest = x)
        RandomNum = x
    return
end