anchor.diff:  !FOX    FRR=FRR*((RV-ELEBET*BETREL)/(ONE-ELEBET*BETREL)) ;
         end if
       endif
  +    write(lout,'(a)')'ELENS> ELENS_KICK_FOX radial kick 03'
   !FOX  YY(1)=YY(1)-(FRR*XI)/RR ;
  +    write(lout,'(a)')'ELENS> ELENS_KICK_FOX radial kick 04'
   !FOX  YY(2)=YY(2)-(FRR*YI)/RR ;
     endif
   
     call comt_daEnd
  +
  +  call dapek(FRR,hh,FRRA)
  +  write(lout,'(2(a,1pe22.15))')'ELENS> ELENS_KICK_FOX computed at RRA=',RRA,' - FRRA=',FRRA
  +
   end subroutine elens_kick_fox
   
   #ifdef CR
  
positive.code:   logical, save          :: elens_lAllowUpdate(nelens) = .true.  ! Flag for disabling updating of kick,
                                                                   !  i.e. after DYNK has touched thetaR2,
                                                                   !  the energy update is disabled.
    logical, save          :: elens_lFox(nelens)         = .false. ! the kick from the elens should be taken into account
                                                                   !  when searching for the closed orbit with DA algebra
  #ifdef CR
    logical, save          :: elens_lAllowUpdate_CR(nelens)
  
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))
  