anchor.diff:      sigma   = zero
     end if
   
  -  if(cn_impact(icoll).gt.0) then
  +  if(cn_impact(icoll) > 0) then
       caverage(icoll) = csum(icoll)/cn_impact(icoll)
   
  -    if((caverage(icoll)**2).gt.(csqsum(icoll)/cn_impact(icoll))) then
  +    if(caverage(icoll)**2 > csqsum(icoll)/cn_impact(icoll)) then
         csigma(icoll) = 0
       else
         csigma(icoll) = sqrt(csqsum(icoll)/cn_impact(icoll) - caverage(icoll)**2)
  
positive.code:   real(kind=fPrec), allocatable, private, save :: sqsum_ay(:)
  
    ! Arrays allocated to npart
    integer,          allocatable, private, save :: part_hit_before_pos(:)
    integer,          allocatable, private, save :: part_hit_before_turn(:)
    integer,          allocatable, private, save :: part_hit_pos(:)
    integer,          allocatable, private, save :: part_hit_turn(:)
    integer,          allocatable, private, save :: part_abs_pos(:)
    integer,          allocatable, public,  save :: part_abs_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_indiv(:)
    real(kind=fPrec), allocatable, private, save :: part_linteract(:)
    real(kind=fPrec), allocatable, private, save :: part_impact(:)
    real(kind=fPrec), allocatable, private, save :: rcx0(:)
    real(kind=fPrec), allocatable, private, save :: rcxp0(:)
    real(kind=fPrec), allocatable, private, save :: rcy0(:)
  
negative.code:   integer j
    real(kind=fPrec) average,sigma,zpj,sum,sqsum
  
    n_impact   = 0
    n_absorbed = 0
    sum        = zero
    sqsum      = zero
  
  #ifdef G4COLLIMATION
    do j=1,napx
  