anchor.diff:  !FOX  D V RE INT ELENOR ;
   !FOX  D V RE INT ELESIG ;
   !FOX  D V RE INT ELEBET ;
  -!FOX  D V RE INT LINLIN ;
  +!FOX  D V DA INT LINLIN NORD NVAR ;
  +!FOX  D V DA INT TMPTMP NORD NVAR ;
  +!FOX  D V RE INT TMPCOF ;
  +!FOX  D V RE INT RKK ;
   !FOX  D V RE INT ONE ; D V RE INT TWO ; D V RE INT HALF ;
   !FOX  D V RE INT RRA ;
   !FOX  E D ;
  
positive.code:       case (3)
          ! RADIAL PROFILE: eLens with radial profile as from file
          ! formula: (cumul_J(r)-cumul_J(r1))/(cumul_J(r2)-cumul_J(r1))
          linlin=lininterp( RRA, &
                elens_radial_profile_R(0:elens_radial_profile_nPoints(elens_iRadial(ielens(ix))),elens_iRadial(ielens(ix))), &
                elens_radial_profile_J(0:elens_radial_profile_nPoints(elens_iRadial(ielens(ix))),elens_iRadial(ielens(ix))), &
                elens_radial_profile_nPoints(elens_iRadial(ielens(ix)))+1)
  !FOX    FRR=(LINLIN-ELEFR1)/ELENOR ;
  
        case default
          write(lout,"(a,i0,a)") "ELENS> ERROR in kickelens: elens_type=",elens_type(ielens(ix))," not recognized. "
          write(lout,"(a)")      "ELENS>       Possible values for type are: 1, 2 and 3"
          call prror
        end select
  
negative.code:   !    frr  if R1 < r < R2
    !    1    if r > R2
    call dapek(RR,hh,RRA)
    if (RRA.gt.elens_r1(iLens)) then ! rr <= r1 -> no kick from elens
      if (RRA.lt.elens_r2(iLens)) then ! r1 < rr < r2
        select case (elens_type(iLens))
        case (1)
          ! UNIFORM: eLens with uniform profile
          ! formula: (r^2-r1^2)/(r2^2-r1^2)
  