anchor.diff:  
     ! IN "+CD DBTRTHIN" and "+CD DBDATEN"
   !  real(kind=fPrec) remitx_dist,remity_dist,
  -! &     remitx_collgap,remity_collgap
  -!  common  /remit/ remitx_dist, remity_dist,
  -! &     remitx_collgap,remity_collgap
   
     integer, private :: k
   !
  
positive.code:   integer, allocatable, save :: cn_absorbed(:) !(max_ncoll)
    real(kind=fPrec), allocatable, save :: caverage(:) !(max_ncoll)
    real(kind=fPrec), allocatable, save :: csigma(:) !(max_ncoll)
  !  common /collsummary/ caverage,csigma,cn_impact,cn_absorbed
  
  ! Change the following block to npart
  ! This is the array that the generated distribution is placed into
  
negative.code:   real(kind=fPrec) :: mingap, gap_h1, gap_h2, gap_h3, gap_h4
    integer :: coll_mingap_id
  
  ! IN "+CD DBTRTHIN", "+CD DBDATEN" and "+CD DBTHIN6D"
  ! logical cut_input
  
  ! IN "+CD DBTRTHIN" and "+CD DBDATEN"
    real(kind=fPrec), save :: remitx_dist,remity_dist,remitx_collgap,remity_collgap
  
    logical, save :: firstcoll,found,onesided
    integer rnd_lux,rnd_k1,rnd_k2
  
  