anchor.diff:          elens_optVal(j)=elens_optVal(j)/sqrt(real(elens_nUpdates(j),fPrec))
         end if
       end if
  -   
  +
  +    lPrint=.false.
       ! compute R1 out of normalised settings
       if (elens_r1(j)<zero) then
         oldVal=elens_r1(j)
  
positive.code:       write(lout,"(a,i0,a)") "ELENS> Recomputing R2 of e-lens #",j," named "//trim(bez(jj))//": "
        write(lout,"(a,1pe22.15)") "ELENS> - original value [sig]=",oldVal
        write(lout,"(a,1pe22.15)") "ELENS> - new value [mm]=",elens_r2(j)
      end if
      ! compute electron sigma out of normalised settings
      if (elens_sig(j)<zero) then
  
negative.code:       write(lout,"(a,i0,a)") "ELENS> Recomputing R1 of e-lens #",j," named "//trim(bez(jj))//": "
        write(lout,"(a,1pe22.15)") "ELENS> - original value [sig]=",oldVal
        write(lout,"(a,1pe22.15)") "ELENS> - new value [mm]=",elens_r1(j)
        lPrint=.true.
      end if
      ! compute R2 out of normalised settings
      if (elens_r2(j)<zero) then
  