anchor.diff:  
     implicit none
   
  -  integer j, fUnit
  +  integer j,fUnit,ieff,ieffdpop
     real(kind=fPrec) gammax,gammay,xdisp,dnormx,dnormy,driftx,drifty,xnorm,xpnorm,xangle,ynorm,ypnorm,&
  -    yangle
  +    yangle,c_dpop,dpopmin,dpopmax
   
   #ifdef HDF5
     ! For tracks2
  
positive.code:         n_tot_absorbed = n_tot_absorbed + 1
          iturn_last_hit = part_hit_before_turn(j)
          iturn_absorbed = part_hit_turn(j)
          if(iturn_last_hit.eq.0) then
            iturn_last_hit = iturn_absorbed
            iturn_survive  = iturn_absorbed - iturn_last_hit
          end if
        end if
  
  
negative.code:   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,zpj,xmax_pencil,ymax_pencil,xmax_nom,ymax_nom,nom_aperture,     &
      scale_bx,scale_by,c_tilt(2),c_offset,c_aperture,c_rotation
  
  #ifdef G4COLLIMATION
    integer :: g4_lostc
  