anchor.diff:    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
  +  logical spErr, tmpl, lfound
   
     call chr_split(inLine, lnSplit, nSplit, spErr)
     if(spErr) then
  
positive.code: 
    ! check that, in case R1/R2 are declared by the user in terms of n-sigma, the provided info is consistent
    do j=1,nelens
      if (elens_r1(j).lt.zero .or. elens_r2(j).lt.zero .or. elens_sig(j).lt.zero) then
        ! printout:
        ! - find name of elens
        do jj=1,nele
  
negative.code:       write(lout,"(a,1pe22.15)") "ELENS> - new value [mm]=",elens_r1(j)
      end if
      ! compute R2 out of normalised settings
      if (elens_r2(j)<zero) then
        oldVal=elens_r2(j)
        if (elens_emin(j)>zero) then ! betatron cut
          elens_r2(j)=abs(elens_r2(j))*sqrt(elens_optVal(j)*elens_emin(j)/(e0f/nucm0))
        else ! momentum cut
          elens_r2(j)=abs(elens_r2(j))*(elens_optVal(j)*elens_sigdpp(j))
  