anchor.diff:    real(kind=fPrec), save :: remitx_dist,remity_dist,remitx_collgap,remity_collgap
   
     logical, save :: firstcoll
  -  integer rnd_lux,rnd_k1,rnd_k2
   
     integer, save :: myix
   
  
positive.code:   ! Initialize random number generator
    if(rnd_seed == 0) rnd_seed = time_getSysClock()
    if(rnd_seed <  0) rnd_seed = abs(rnd_seed)
    rnd_lux = 3
    rnd_k1  = 0
    rnd_k2  = 0
    call rluxgo(rnd_lux, rnd_seed, rnd_k1, rnd_k2)
    write(outlun,*) 'INFO>  rnd_seed: ', rnd_seed
  
    ! Call distribution routines only if collimation block is in fort.3
  
negative.code: 
    ! Intialise random generator with offset_seed
    c_offsettilt_seed = abs(c_offsettilt_seed)
    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
  