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:       end if
        if ( elens_lFull(ielens(iElem)) ) then
          if ( elens_lZeroThick(ielens(iElem)) ) then
            write(lout,"(a)") "ELENS> Elens '"//trim(bez(iElem))//"' is actually a wire."
          else
            write(lout,"(a)") "ELENS> Elens '"//trim(bez(iElem))//"' is full (i.e. R1=0.0, R2>0.0)."
          end if
  