anchor.diff:  !FOX  D V RE INT ELENOR ;
   !FOX  D V RE INT ELESIG ;
   !FOX  D V RE INT ELEBET ;
  -!FOX  D V RE INT LINLIN ;
  +!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 ;
  
positive.code:   call comt_daStart
  !-----------------------------------------------------------------------
  
    XCLO=elens_offset_x(ielens(ix))
    YCLO=elens_offset_y(ielens(ix))
    ELE_R1=elens_r1(ielens(ix))
    ELE_R2=elens_r2(ielens(ix))
    ELENOR=elens_geo_norm(ielens(ix))
    ELESIG=elens_sig(ielens(ix))
    ELEFR1=elens_radial_fr1(ielens(ix))
    ELEBET=elens_beta_e(ielens(ix))
    
    ! 1) apply offset of e-lens
    !    xx = x(proton) - elens_offset_x
  
negative.code:   !    frr  if R1 < r < R2
    !    1    if r > R2
    call dapek(RR,hh,RRA)
    if (RRA.gt.elens_r1(iLens)) then ! rr <= r1 -> no kick from elens
      if (RRA.lt.elens_r2(iLens)) then ! r1 < rr < r2
        select case (elens_type(iLens))
        case (1)
          ! UNIFORM: eLens with uniform profile
          ! formula: (r^2-r1^2)/(r2^2-r1^2)
  