anchor.diff:  
     ! Intialise random generator with offset_seed
     c_offsettilt_seed = abs(c_offsettilt_seed)
  -  rnd_lux = 3
  -  rnd_k1  = 0
  -  rnd_k2  = 0
  -  call rluxgo(rnd_lux, c_offsettilt_seed, rnd_k1, rnd_k2)
  +  call rluxgo(3, c_offsettilt_seed, 0, 0)
   
     ! Generate random tilts (Gaussian distribution plus systematic)
     if(c_rmstilt_prim > zero .or. c_rmstilt_sec > zero .or. c_systilt_prim /= zero .or. c_systilt_sec /= zero) then
  
positive.code:   ! Alternatively, we can use ranecu instead, which is capable of continuing a chain of random numbers from
    ! a given set of seeds.
    ! It is probably unnecessary to use different random seeds here in the first place.
    rnd_lux = 3
    rnd_k1  = 0
    rnd_k2  = 0
    call rluxgo(rnd_lux, rnd_seed, rnd_k2, rnd_k2)
    dummy = rndm5(1) ! Reset rndm5 too
  
  end subroutine collimate_start
  
negative.code:   ! Initialize random number generator
    if(rnd_seed == 0) rnd_seed = time_getSysClock()
    if(rnd_seed <  0) rnd_seed = abs(rnd_seed)
    call rluxgo(3, rnd_seed, 0, 0)
    write(outlun,*) 'INFO>  rnd_seed: ', rnd_seed
  
    ! Call distribution routines only if collimation block is in fort.3
  