anchor.diff:  
     integer j,fUnit,ieff,ieffdpop
     real(kind=fPrec) gammax,gammay,xdisp,dnormx,dnormy,driftx,drifty,xnorm,xpnorm,xangle,ynorm,ypnorm,&
  -    yangle,c_dpop,dpopmin,dpopmax
  +    yangle,c_dpop,dpopmin,dpopmax,nspx,nspy
   
   #ifdef HDF5
     ! For tracks2
  
positive.code:   real(kind=fPrec), private, save :: scale_by0 = zero
    real(kind=fPrec), private, save :: bx_dist   = zero
    real(kind=fPrec), private, save :: by_dist   = zero
    real(kind=fPrec), private, save :: nspx      = zero
    real(kind=fPrec), private, save :: nspy      = zero
  
  #ifdef HDF5
    ! Variables to save hdf5 dataset indices
  
negative.code: #endif
  
    integer j
    real(kind=fPrec) gammax,gammay,xdisp,nspxd,xj,xpj,yj,ypj,pj,nspx,nspy
  
  #ifdef HDF5
    integer hdfturn,hdfpid,hdftyp
  