anchor.diff:    epsilon=c1m15
     gteps=one+epsilon
     lteps=one-epsilon
  +  cl2ori=c1m7
   
     iLens=ielens(ix)
     XFFSET=elens_offset_x(iLens)
  
positive.code: end subroutine normaliseRadialProfile
  
  ! ================================================================================================ !
  !  Last modified: 2019-04-11
  !  compute kick from electron lens
  ! ================================================================================================ !
  subroutine elens_kick(i,ix,n)
  
negative.code:   integer, save               :: elens_iSet_def=0           ! default iSet
    logical, save               :: elens_lFox_def=.true.      ! default lFox
    integer, save               :: elens_radial_mpoints_def=2 ! default number of points for interpolating radial profiles from ASCI files
    integer, save               :: elens_radial_mpoints_ori=3 ! default minimum number of points for interpolating radial profiles from ASCI files in case we are very close to the origin
  
    ! beam of the lens
    real(kind=fPrec), save      :: elens_beam_mass_def=pmae   ! default mass of lens beam [MeV/c2]
  