anchor.diff:    ! If not, we're doing something, so start the stop watch
     call time_startClock(time_clockSCAT)
   
  -  ! Store the seeds in the randum number generator, and set ours
  +  ! Store the seeds in the random number generator, and set ours
     call recuut(tmpSeed1,tmpSeed2)
     call recuin(scatter_seed1,scatter_seed2)
   
  
positive.code:   end if
  
    if(updateE) then
      call part_updatePartEnergy(2)
    end if
  
    ! Restore seeds in random generator
  
negative.code:       llostp(j) = .true.
        phi       = zero
      else
        ! Calculate new energy, and scattering angle
        if(abs(dPP) >= pieni) then ! Update the energy
          ejv(j)  = sqrt(((one+dPP)*ejfv(j))**2 + nucm(j)**2)
          updateE = .true.
        end if
  
        ! Update particle trajectory
  