anchor.diff:          case (2)
             ! GAUSSIAN: eLens with Gaussian profile
             if ( elens_lFull(ielens(ix)) ) then
  -            if ( rr <= c1m7 ) then
  +            if ( rr <= cl2ori ) then
                 ! formula: r*Rref/2sig^2/(exp(-r1^2/2sig^2)-exp(-r2^2/2sig^2))
                 frr=(((rr/elens_sig(ielens(ix)))*(elens_rref(ielens(ix))/elens_sig(ielens(ix))))/two)/elens_geo_norm(ielens(ix))
               else
  
positive.code:       case (2)
          ! GAUSSIAN: eLens with Gaussian profile
          if ( elens_lFull(iLens) ) then
            if ( RRA <= c1m7 ) then
              ! formula: r*Rref/2sig^2/(exp(-r1^2/2sig^2)-exp(-r2^2/2sig^2))
  !FOX        FRR=(((RR/ELESIG)*(ELE_RR/ELESIG))/TWO)/ELENOR ;
              if (st_debug) write(lout,'(a)') "ELENS> ELENS_KICK_FOX: elens type GAUSSIAN, linearised formula"
  
negative.code:   integer, save               :: elens_iSet_def=0           ! default iSet
    logical, save               :: elens_lFox_def=.true.      ! default lFox
    integer, save               :: elens_radial_mpoints_def=2 ! default number of points for interpolating radial profiles from ASCI files
    integer, save               :: elens_radial_mpoints_ori=3 ! default minimum number of points for interpolating radial profiles from ASCI files in case we are very close to the origin
  
    ! beam of the lens
    real(kind=fPrec), save      :: elens_beam_mass_def=pmae   ! default mass of lens beam [MeV/c2]
  