anchor.diff:    integer, save               :: elens_iSet_def=0           ! default iSet
     logical, save               :: elens_lFox_def=.true.      ! default lFox
     integer, save               :: elens_radial_mpoints_def=2 ! default number of points for interpolating radial profiles from ASCI files
  +
  +  ! beam of the lens
  +  real(kind=fPrec), save      :: elens_beam_mass_def=pmae   ! default mass of lens beam [MeV/c2]
  +  real(kind=fPrec), save      :: elens_beam_chrg_def=-one   ! default charge of lens beam [e]
     
     ! variables to save elens parameters for tracking etc.
     integer, allocatable, save          :: elens_type(:)            ! integer for elens type (nelens)
  
positive.code: !FOX    FRR=ZERO ;
        call dapok(FRR,hh,FRRA)
      end if
  !FOX  YY(1)=YY(1)-(FRR*XI)/RR ;
  !FOX  YY(2)=YY(2)-(FRR*YI)/RR ;
    end if
  
    if (st_debug) then
  
negative.code:   call alloc(elens_iSet           , nelens_new,           elens_iSet_def, 'elens_iSet'           )
    call alloc(elens_optVal         , nelens_new,                     zero, 'elens_optVal'         )
    call alloc(elens_nUpdates       , nelens_new,                        0, 'elens_nUpdates'       )
    call alloc(elens_beam_mass      , nelens_new,      elens_beam_mass_def, 'elens_beam_mass'      )
    call alloc(elens_beam_chrg      , nelens_new,      elens_beam_chrg_def, 'elens_beam_chrg'      )
    call alloc(elens_iRadial        , nelens_new,                        0, 'elens_iRadial'        )
    call alloc(elens_radial_fr1     , nelens_new,                     zero, 'elens_radial_fr1'     )
    call alloc(elens_radial_fr2     , nelens_new,                     zero, 'elens_radial_fr2'     )
  