anchor.diff:          call sixin_echoVal("L", elens_len(ielens(iElem)),"ELENS",iLine)
           call sixin_echoVal("I", elens_I(ielens(iElem)),  "ELENS",iLine)
           call sixin_echoVal("Ek",elens_Ek(ielens(iElem)), "ELENS",iLine)
  +      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 (i.e. R1=R2=0.0)."
  +        else
  +          write(lout,"(a)") "ELENS> Elens '"//trim(bez(iElem))//"' is full (i.e. R1=0.0, R2>0.0)."
  +        end if
  +      else 
  +        if ( elens_lZeroThick(ielens(iElem)) ) then
  +          write(lout,"(a)") "ELENS> Elens '"//trim(bez(iElem))//"' is a zero-thickness hollow electron lens (i.e. R1=R2>0.0)."
  +        else
  +          write(lout,"(a)") "ELENS> Elens '"//trim(bez(iElem))//"' is a regular hollow electron lens (i.e. R1>0.0, R2>0.0, R1<R2)."
  +        end if
         end if
       end if
    
  
positive.code: ! ================================================================================================ !
  subroutine elens_kick(i,ix,n)
  
    use mod_common, only : beta0, napx
    use mod_common_main, only : xv1, xv2, yv1, yv2, moidpsv, rvv
    use mathlib_bouncer
    use numerical_constants, only : zero, one, two, c1m15, c1m7
  
negative.code:          
        case (2)
          ! GAUSSIAN: eLens with Gaussian profile
          if ( elens_lFull(iLens) ) then
            if ( RRA <= c1m7 ) then
              ! formula: r*r2/2sig^2/(exp(-r1^2/2sig^2)-exp(-r2^2/2sig^2))
  !FOX        FRR=(((RR/ELESIG)*(ELE_R2/ELESIG))/TWO)/ELENOR ;
  