anchor.diff:    ! (x,y)<r1 or ( (xx>r2) || (yy>r2) ): no kick from lens
     if (.not. ( (ax.lt.r1lteps.and.ay.lt.r1lteps) .or. (ax.gt.r2gteps.or.ay.gt.r2gteps) ) ) then
        
  -    if (st_debug) then
  -      write(lout,'(2(a,1pe22.15))')'CHEBY> CHEBY_KICK_FOX computing at XA=',XA,' - YA=',YA
  -      write(lout,'(2(a,1pe22.15))')'CHEBY>                when R1=',cheby_r1(icheby(ix)),&
  -                                                          'and R2=',cheby_r2(icheby(ix))
  -    end if
  -      
  +    if (st_debug) write(lout,'(2(a,1pe23.16))')'CHEBY> CHEBY_KICK_FOX computing kick when R1[mm]= ',&
  +         cheby_r1(icheby(ix)),' and R2[mm]= ',cheby_r2(icheby(ix))
  +    
       ! compute kick
       
       ! normalised variables
  
positive.code: 
  ! Do not remove or modify the comment below.
  !     DADAL AUTOMATIC INCLUSION
    
  end subroutine cheby_kick_fox
  
  
  
negative.code:   use numerical_constants, only : zero, one, c180e0, pi, c1m15, two, pieni
    use physical_constants, only: clight
    use mod_lie_dab, only : lnv, idao, rscrri, iscrda
    use mod_common_da
    
    implicit none
  