anchor.diff:    use mod_common_main
     use mathlib_bouncer
     use numerical_constants, only : zero, one
  -  use utils, only : lininterp
  +  use utils, only : huntBin, polcof
     use mod_lie_dab, only : lnv, idao, rscrri, iscrda
     use mod_common_track, only : xxtr, yytr, issss, comt_daStart, comt_daEnd
     use mod_common_da
  
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:   integer, intent(in) :: i
    integer, intent(in) :: ix
    
    integer          :: jj, idaa, iLens, iRadial, nBin, nPoints, mPoints, kMin, kMax, kk
    integer          :: hh(lnv)=0
    real(kind=fPrec) :: rra, xclo, yclo, ele_r1, ele_r2, elenor, elesig, elefr1, elebet, rkk, tmpcof
    real(kind=fPrec), allocatable :: cof(:)
  
    common/daele/alda,asda,aldaq,asdaq,smida,xx,yy,dpda,dpda1,sigmda,ej1,ejf1,rv
    save
  