anchor.diff:    real(kind=fPrec), allocatable, private, save :: sqsum_ay(:)
   
     ! Arrays allocated to npart
  -  integer,          allocatable, private, save :: part_hit_pos(:)         ! Hit flag for last hit
  -  integer,          allocatable, private, save :: part_hit_turn(:)        ! Hit flag for last hit
  -  integer,          allocatable, private, save :: part_abs_pos(:)         ! Absorbed in element
  -  integer,          allocatable, public,  save :: part_abs_turn(:)        ! Absorbed in turn
  -  integer,          allocatable, private, save :: part_hit_before_pos(:)
  -  integer,          allocatable, private, save :: part_hit_before_turn(:)
  +  integer,          allocatable, private, save :: part_hit_pos(:)   ! Hit flag for last hit
  +  integer,          allocatable, private, save :: part_hit_turn(:)  ! Hit flag for last hit
  +  integer,          allocatable, private, save :: part_abs_pos(:)   ! Absorbed in element
  +  integer,          allocatable, public,  save :: part_abs_turn(:)  ! Absorbed in turn
     integer,          allocatable, private, save :: part_select(:)
     integer,          allocatable, private, save :: nabs_type(:)
     integer,          allocatable, private, save :: nhit_type(:)
     real(kind=fPrec), allocatable, private, save :: part_linteract(:)
  -  real(kind=fPrec), allocatable, private, save :: part_indiv(:)            ! Divergence of impacting particles
  -  real(kind=fPrec), allocatable, private, save :: part_impact(:)           ! Impact parameter (0 for inner face)
  +  real(kind=fPrec), allocatable, private, save :: part_indiv(:)     ! Divergence of impacting particles
  +  real(kind=fPrec), allocatable, private, save :: part_impact(:)    ! Impact parameter (0 for inner face)
     real(kind=fPrec), allocatable, private, save :: rcx0(:)
     real(kind=fPrec), allocatable, private, save :: rcxp0(:)
     real(kind=fPrec), allocatable, private, save :: rcy0(:)
  
positive.code: 
    real(kind=fPrec), intent(in) :: c_aperture
    real(kind=fPrec), intent(in) :: c_rotation
  
    integer j
  
  
negative.code: 
    integer j, iSlice, nSlices
    logical onesided, linside(napx), isAbs, isHit
    real(kind=fPrec) nsig,c_length,jawLength,jawAperture,jawOffset,jawTilt(2),x_Dump,xpDump,y_Dump,   &
      ypDump,s_Dump,xmax,ymax,calc_aperture,zpj,xmax_pencil,ymax_pencil,xmax_nom,ymax_nom,            &
      nom_aperture,scale_bx,scale_by,c_tilt(2),c_offset,c_aperture,c_rotation,cry_bendangle,cry_tilt, &
      cry_tilt0
  
    call time_startClock(time_clockCOLL)
  
  