anchor.diff:    ! R1/R2 expressed in normalised units (either betatron or dispersion)
     real(kind=fPrec), save      :: elens_emin_def=-one      ! default normalised emittance [m rad]
     real(kind=fPrec), save      :: elens_sigdpp_def=-one    ! default rms of distribution in delta_p / p []
  +  integer, save               :: elens_iSet_def=0         ! default iSet
  +  
   
     ! variables to save elens parameters for tracking etc.
     integer, allocatable, save          :: elens_type(:)            ! integer for elens type (nelens)
  
positive.code:         write(lerr,"(a)") "ELENS>       please choose one among min|max|ave|qve"
          call prror
        end if
      end if
    end do
    
  
negative.code:       else ! momentum cut
          elens_sig(j)=abs(elens_sig(j))*(elens_optVal(j)*elens_sigdpp(j))
        end if
        elens_sig(j)=elens_sig(j)*c1e3
        ! ...and printout:
        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)
        write(lout,"(a,1pe22.15)") "ELENS> - momentum of reference particle [MeV/c]=",e0f
        write(lout,"(a,1pe22.15)") "ELENS> - mass of reference particle [MeV/c2]=",nucm0
        write(lout,"(a,1pe22.15)") "ELENS> - rms of delta distributione []=",elens_sigdpp(j)
        write(lout,"(a,1pe22.15)") "ELENS> - beta/disp function [m]=",elens_optVal(j)
      end if
  
      ! compute geometrical normalisation factor
  