anchor.diff:  
     if (st_debug) then
       call dapek(FRR,hh,FRRA)
  -    write(lout,'(2(a,1pe23.16))')'ELENS> ELENS_KICK_FOX computed at RRA=',RRA,' - FRRA=',FRRA
  +    write(lout,'(2(a,1pe23.16))')'ELENS> ELENS_KICK_FOX computed at RRA= ',RRA,' - FRRA= ',FRRA
       write(lout,'(a)')'ELENS> ELENS_KICK_FOX closed orbit AFTER elens:'
       call dapri(XX(1),6)
       call dapri(XX(2),6)
  
positive.code:       write(lout,"(a,i0,a)") "ELENS> Recomputing sigma of e-beam 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_sig(j)
      end if
      if(st_debug) then
        write(lout,"(a)"         ) "ELENS> ...using the following parameters:"
        write(lout,"(a,1pe22.15)") "ELENS> - normalised emittance [m rad]=",elens_emin(j)
        write(lout,"(a,1pe22.15)") "ELENS> - geometrical emittance [m rad]=",elens_emin(j)/(e0f/nucm0)
  
negative.code:         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)
  