anchor.diff:    integer, intent(in) :: i
     integer, intent(in) :: ix
     
  -  integer          :: jj, idaa
  +  integer          :: jj, idaa, iLens, iRadial, nBin, nPoints, mPoints, kMin, kMax, kk
     integer          :: hh(lnv)=0
  -  real(kind=fPrec) linlin, rra, xclo, yclo, ele_r1, ele_r2, elenor, elesig, elefr1, elebet
  +  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
  
positive.code:   !    frr  if R1 < r < R2
    !    1    if r > R2
    call dapek(RR,hh,RRA)
    if (RRA.gt.elens_r1(ielens(ix))) then ! rr <= r1 -> no kick from elens
      if (RRA.lt.elens_r2(ielens(ix))) then ! r1 < rr < r2
        select case (elens_type(ielens(ix)))
        case (1)
          ! UNIFORM: eLens with uniform profile
          ! formula: (r^2-r1^2)/(r2^2-r1^2)
  
negative.code: !FOX  D V RE INT ELENOR ;
  !FOX  D V RE INT ELESIG ;
  !FOX  D V RE INT ELEBET ;
  !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 ;
  