anchor.diff:  
   !  private
   
  -!+cd collpara
     integer, parameter :: max_ncoll  = 100
  -  !integer, parameter :: maxn       = 20000
  + !integer, parameter :: maxn       = 20000
     integer, parameter :: numeff     = 32
     integer, parameter :: numeffdpop = 29
     integer, parameter :: nc         = 32
   
  -!+cd collMatNum
  -! EQ 2016 added variables for collimator material numbers
     integer, parameter :: nmat  = 14
     integer, parameter :: nrmat = 12
   
  -!+cd database
  -!GRD THIS BLOC IS COMMON TO MAINCR, DATEN, TRAUTHIN AND THIN6D
  -  logical, save :: do_coll = .false.
  -  logical, save :: do_select
  -  logical, save :: do_nominal
  -  logical, save :: dowrite_dist
  -  logical, save :: do_oneside
  -  logical, save :: dowrite_impact
  -  logical, save :: dowrite_secondary
  -  logical, save :: dowrite_amplitude
  -  logical, save :: radial
  -  logical, save :: systilt_antisymm
  -  logical, save :: dowritetracks
  -  logical, save :: cern
  -  logical, save :: do_mingap
  -
  -!SEPT2005 for slicing process
  -  integer, save :: nloop    = 1
  -  integer, save :: rnd_seed = 0
  -  integer, save :: c_offsettilt_seed
  -  integer, save :: ibeam
  -  integer, save :: jobnumber
  -  integer, save :: do_thisdis = 0
  -  integer, save :: pencil_distr
  +  ! Logical Flags
  +  logical, public,  save :: do_coll           = .false.
  +  logical, private, save :: do_select         = .false.
  +  logical, private, save :: do_nominal        = .false.
  +  logical, private, save :: dowrite_dist      = .false.
  +  logical, private, save :: do_oneside        = .false.
  +  logical, private, save :: dowrite_impact    = .false.
  +  logical, private, save :: dowrite_secondary = .false.
  +  logical, private, save :: dowrite_amplitude = .false.
  +  logical, private, save :: radial            = .false.
  +  logical, private, save :: systilt_antisymm  = .false.
  +  logical, private, save :: dowritetracks     = .false.
  +  logical, private, save :: cern              = .false.
  +  logical, private, save :: do_mingap         = .false.
  +
  +  integer, private, save :: nloop      = 1
  +  integer, private, save :: rnd_seed   = 0
  +  integer, private, save :: ibeam      = 1
  +  integer, private, save :: jobnumber  = 0
  +  integer, private, save :: do_thisdis = 0
   
     real(kind=fPrec), private, save :: myenom = zero
     real(kind=fPrec), private, save :: mynex  = zero
  
positive.code: 
    select case(lnSplit(1))
  
    case("DOCOLL")
      do_coll = .true.
  
    case("ENERGY")
      if(nSplit /= 2) then
  
negative.code:       iErr = .true.
        return
      end if
      filename_dis = trim(lnSplit(2))
  
    case("NSIG_FAM")
      if(nSplit /= 3) then
  