anchor.diff:      else ! r1 < r2 <= rr
   !FOX    FRR=ONE ;
       endif
  +    write(lout,'(a)')'ELENS> ELENS_KICK_FOX radial kick 01'
       ! 'radial kick'
   !FOX    FRR=(((ELE_R2*ELE_R2)/RR)*FRR)*(MTCDA/(ONE+DPDA)) ;
  +    write(lout,'(a)')'ELENS> ELENS_KICK_FOX radial kick 02'
       if(elens_lThetaR2(ielens(ix))) then
         if(elens_I(ielens(ix)) < zero) then
   !FOX    FRR=FRR*((RV+ELEBET*BETREL)/(ONE+ELEBET*BETREL)) ;
  
positive.code: !FOX    LINLIN=TMPCOF ;
  !FOX    TMPTMP=RR;
          do kk=2,mPoints
             TMPCOF=COF(kk)
             RKK=real(kk-1,fPrec)
  !FOX       LINLIN=LINLIN+(TMPTMP*TMPCOF)/RKK ;
  !FOX       TMPTMP=TMPTMP*RR ;
          end do
  !FOX    FRR=(LINLIN-ELEFR1)/ELENOR ;
          call dealloc(cof,'cof')
  
negative.code:         ! UNIFORM: eLens with uniform profile
          ! formula: (r^2-r1^2)/(r2^2-r1^2)
  !FOX    FRR=((RR+ELE_R1)*(RR-ELE_R1))/ELENOR ;
  
        case (2)
          ! GAUSSIAN: eLens with Gaussian profile
          ! formula: (exp(-r1^2/2sig^2)-exp(-r^2/2sig^2))/(exp(-r1^2/2sig^2)-exp(-r2^2/2sig^2))
  