anchor.diff:    use geant4
   #endif
   
  -  implicit none
  -
  -  integer :: j
  +  real(kind=fPrec), intent(in) :: stracki
  +
  +  integer j
  +  real(kind=fPrec) average,sigma
   
   #ifdef HDF5
     ! For tracks2
  
positive.code:   if (n_impact.gt.0) then
      average = sum/n_impact
  
      if (sqsum/n_impact.ge.average**2) then
        sigma = sqrt(sqsum/n_impact - average**2)
      else
        sigma = zero
  
negative.code: 
    integer j, iSlice, nSlices
    logical onesided, linside(napx)
    real(kind=fPrec) jawLength,jawAperture,jawOffset,jawTilt(2),x_Dump,xpDump,y_Dump,ypDump,s_Dump,   &
      xmax,ymax,calc_aperture,ldrift,c_nex2,c_ney2,Nap1pos,Nap2pos,Nap1neg,Nap2neg,tiltOffsPos1,      &
      tiltOffsPos2,tiltOffsNeg1,tiltOffsNeg2,beamsize1,beamsize2,betax1,betax2,betay1,betay2,alphax1, &
      alphax2,alphay1,alphay2,minAmpl
  
  #ifdef G4COLLIMATION
    integer :: g4_lostc
  