anchor.diff:                                                                     ! 1 : uniform profile
                                                                      ! 2 : Gaussian profile
                                                                      ! 3 : radial profile from file
  +                                                                   ! 4 : degenerate case: wire
     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)
  
positive.code:         iErr = .true.
          return
        end if
      case default
        write(lerr,"(a)") "ELENS> ERROR Elens type '"//trim(lnSplit(2))//"' not recognized. Remember to use all UPPER CASE."
        iErr = .true.
  
negative.code:     !    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.and..not.elens_lZeroThick(ielens(ix)) ) then ! R1<=rr<R2: type-dependent kick
            
          select case (elens_type(ielens(ix)))
              
  