anchor.diff:      !    0    if r <= R1
       !    frr  if R1 < r < R2
       !    1    if r >= R2
  -    if (rr.gt.elens_r1(ielens(ix)).or.(elens_lZeroThick(ielens(ix)).and.rr.eq.elens_r1(ielens(ix)))) then ! rr<=R1 -> no kick from elens
  -      if (rr.lt.elens_r2(ielens(ix))) then ! R1<rr<R2
  +    if (rr>elens_r1(ielens(ix)).or.(elens_lZeroThick(ielens(ix)).and.rr.eq.elens_r1(ielens(ix)))) then ! rr<=R1 -> no kick from elens
  +      if (rr<elens_r2(ielens(ix))) then ! R1<rr<R2
           select case (elens_type(ielens(ix)))
           case (1)
             ! UNIFORM: eLens with uniform profile
  
positive.code:   character(len=mStrLen) tmpch
    real(kind=fPrec) tmpflt
    integer nSplit, iElem, j, tmpi1, tmpi2, tmpi3
    logical spErr, tmpl, lfound, tmplFox
    real(kind=fPrec), save :: tmpEmin, tmpSigDpp
    integer, save :: tmpSet, tmpMpoints
  
    call chr_split(inLine, lnSplit, nSplit, spErr)
    if(spErr) then
  
negative.code:       end if
      end do
  
      if (elens_r1(j)<zero.or.elens_r2(j)<zero.or.elens_sig(j)<zero) then
        if (elens_iSet(j).eq.3 ) then
          elens_optVal(j)=elens_optVal(j)/real(elens_nUpdates(j),fPrec)
        elseif (elens_iSet(j).eq.4 ) then
  