anchor.diff:  
     select case(lnSplit(1))
   
  -  case("DOCOLL")
  -    do_coll = .true.
  +  case("DO_COLL")
  +    if(nSplit /= 2) then
  +      write(lout,"(a,i0)") "COLL> ERROR DOCOLL expects 1 value, got ",nSplit-1
  +      write(lout,"(a)")    "COLL>       DOCOLL true|false"
  +      iErr = .true.
  +      return
  +    end if
  +    call chr_cast(lnSplit(2), do_coll, iErr)
   
     case("ENERGY")
       if(nSplit /= 2) then
  
positive.code:   real(kind=fPrec), private, save :: jaw_fit(6,6) = zero
    real(kind=fPrec), private, save :: jaw_ssf(2)   = zero
  
  
  
    !SEPT2005,OCT2006 added offset
    real(kind=fPrec), private, save :: &
    &xbeat,xbeatphase,ybeat,ybeatphase,                                &
    &c_rmstilt_prim,c_rmstilt_sec,c_systilt_prim,c_systilt_sec,        &
    &c_rmsoffset_prim,c_rmsoffset_sec,c_sysoffset_prim,                &
    &c_sysoffset_sec,c_rmserror_gap,ndr,                            &
    &driftsx,driftsy,pencil_offset,pencil_rmsx,pencil_rmsy,            &
    &sigsecut3,sigsecut2,enerror,bunchlength
  
    ! From collimation_comnul
    real(kind=fPrec), public,  save :: emitnx0_dist    = zero
  
negative.code: 
  !  private
  
    integer, parameter :: max_ncoll  = 100
   !integer, parameter :: maxn       = 20000
    integer, parameter :: numeff     = 32
    integer, parameter :: numeffdpop = 29
    integer, parameter :: nc         = 32
  
    integer, parameter :: nmat  = 14
    integer, parameter :: nrmat = 12
  
    ! 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
  