anchor.diff:    integer, intent(in) :: ix
     integer, intent(in) :: n
   
  -  real(kind=fPrec) xx, yy, ax, ay, rr, dxp, dyp
  +  real(kind=fPrec) xx, yy, ax, ay, rr, dxp, dyp, uu, vv
     real(kind=fPrec) theta, angle_rad, epsilon, gteps, lteps, r1lteps, r2gteps
  -  integer          jj
  +  integer          ii, jj, nn
     logical          lrotate
   
  +  real(kind=fPrec) :: Tx (0:cheby_maxOrder(cheby_itable(icheby(ix)))), &
  +                      Ty (0:cheby_maxOrder(cheby_itable(icheby(ix)))), &
  +                      Tpx(0:cheby_maxOrder(cheby_itable(icheby(ix)))), &
  +                      Tpy(0:cheby_maxOrder(cheby_itable(icheby(ix))))
  +  
     epsilon=c1m15
     gteps=one+epsilon
     lteps=one-epsilon
  
positive.code:     if ( (ax < r1lteps .and. ay < r1lteps) .or. (ax > r2gteps .or. ay > r2gteps) ) cycle
  
      ! compute kick from cheby map
      call cheby_getKick( xx, yy, dxp, dyp, cheby_itable(icheby(ix)) )
      if (abs(dxp)<epsilon) dxp=zero
      if (abs(dyp)<epsilon) dyp=zero
      if ( dxp == zero .and. dyp == zero ) cycle
  
negative.code: end subroutine cheby_potentialMap
  
  
  subroutine cheby_getPotential( xx, yy, zz, iTable )
  
    ! A. Mereghetti (CERN, BE-ABP-HSS)
  