anchor.diff:        arcdx = 2.5_fPrec
         arcbetax = c180e0
   
  -      if(xj.le.0.) then
  +      if(xj <= 0.) then
           xdisp = xj + (pj-myenom)/myenom * arcdx * sqrt(tbetax(ie)/arcbetax)
         else
           xdisp = xj - (pj-myenom)/myenom * arcdx * sqrt(tbetax(ie)/arcbetax)
  
positive.code:       call h5_finaliseWrite(coll_hdf5_survival)
      else
  #endif
        write(survival_unit,'(2i7)') iturn, nsurvive
  #ifdef HDF5
      end if
  #endif
  
  #ifdef ROOT
      if(root_flag .and. root_Collimation.eq.1) then
        call SurvivalRootWrite(iturn, nsurvive)
      end if
  #endif
  
      if (iturn.eq.numl) then
        nsurvive_end = nsurvive_end + nsurvive
      end if
    end if
  
negative.code:   integer, private, save :: impact_unit, tracks2_unit, pencilbeam_distr_unit, coll_ellipse_unit, all_impacts_unit
    integer, private, save :: FLUKA_impacts_unit, FLUKA_impacts_all_unit, coll_scatter_unit, FirstImpacts_unit, RHIClosses_unit
    integer, private, save :: twisslike_unit, sigmasettings_unit, distsec_unit, efficiency_unit, efficiency_dpop_unit
    integer, private, save :: coll_summary_unit, amplitude_unit, amplitude2_unit, betafunctions_unit, orbitchecking_unit
    integer, private, save :: CollPositions_unit, all_absorptions_unit, efficiency_2d_unit
    integer, private, save :: collsettings_unit, outlun
  
  