anchor.diff:          rate = rate - nkill
   
           ! Additional stochasitic death?
  -        call random_number(r)
  +        r = genrand_real2()
           if (abs(rate) > r) then
               if (rate > 0.0_dp) then
                   nkill = nkill + 1
  
positive.code:         spawned_walker_population = 0
  
          ! Initialise random numbers.
          call random_seed()
  
          ! Now we need to set the reference determinant.
          ! We choose the determinant with the lowest Hamiltonian matrix element.
  
negative.code:             rate = abs(Tau*hamil(iwalker,j))
              nspawn = int(rate)
              rate = rate - nspawn
              r = genrand_real2()
              if (rate > r) nspawn = nspawn + 1
  
              ! Create particles.
  