anchor.diff:      !    frr  if R1 < r < R2
       !    1    if r >= R2
       if ( rr > elens_r1(ielens(ix))*lteps .or. elens_lFull(ielens(ix)) ) then ! rr<R1: no kick from elens
  -      if ( rr < elens_r2(ielens(ix))*lteps ) then ! R1<=rr<R2: type-dependent kick
  +      if ( rr < elens_r2(ielens(ix))*lteps.and..not.elens_lZeroThick(ielens(ix)) ) then ! R1<=rr<R2: type-dependent kick
             
           select case (elens_type(ielens(ix)))
               
  
positive.code:                                                                    ! 1 : uniform profile
                                                                     ! 2 : Gaussian profile
                                                                     ! 3 : radial profile from file
    real(kind=fPrec), allocatable, public, save :: elens_theta_r2(:) ! kick strength at R2 [mrad] (nelens)
    real(kind=fPrec), allocatable, save :: elens_r2(:)               ! outer radius R2 [mm] (nelens)
    real(kind=fPrec), allocatable, save :: elens_r1(:)               ! inner radius R1 [mm] (nelens)
  
negative.code:         iErr = .true.
          return
        end if
      case("WIRE")
        elens_type(ielens(iElem)) = 4
      case default
        write(lerr,"(a)") "ELENS> ERROR Elens type '"//trim(lnSplit(2))//"' not recognized. Remember to use all UPPER CASE."
        iErr = .true.
  