anchor.diff:    use crcoall
     use mod_alloc
     use numerical_constants, only : zero
  -
  -  implicit none
  -
  -  integer, allocatable, save  :: ielens(:) !(nele)    ! index of elens
  -  integer, save               :: nelens=0             ! number of elenses lenses actually in memory
  -  integer, parameter          :: elens_kz=29          ! kz of electron lenses
  -  integer, parameter          :: elens_ktrack=63      ! ktrack of electron lenses
  +  implicit none
  +
  +  integer, allocatable, save  :: ielens(:) !(nele)        ! index of elens
  +  integer, save               :: nelens=0                 ! number of elenses lenses actually in memory
  +  integer, save               :: nelens_radial_profiles=0 ! number of radial profiles available in memory
  +  integer, parameter          :: elens_kz=29              ! kz of electron lenses
  +  integer, parameter          :: elens_ktrack=63          ! ktrack of electron lenses
   
     ! variables to save elens parameters for tracking etc.
     integer, allocatable, save          :: elens_type(:)            ! integer for elens type (nelens)
  
positive.code:                                                                   !    i.e. linear interpolation
    integer, allocatable, save          :: elens_radial_jguess(:)   ! bin for guessed search (nelens)
    ! - file handling and data storage:
    integer, parameter     :: nelens_radial_profiles=20 ! max number of radial profiles
    integer, save          :: melens_radial_profiles    ! radial profiles available in memory
    integer, parameter     :: elens_radial_dim=500      ! max number of points in radial profiles
    character(len=mFileName), save:: elens_radial_filename(nelens_radial_profiles) ! names
    real(kind=fPrec), save :: elens_radial_profile_R(0:elens_radial_dim,nelens_radial_profiles) ! [mm]
    real(kind=fPrec), save :: elens_radial_profile_J(0:elens_radial_dim,nelens_radial_profiles) ! [A]
    integer, save          :: elens_radial_profile_nPoints(nelens_radial_profiles)=0
  
  contains
  
  
negative.code:   character(len=:), allocatable   :: lnSplit(:)
    character(len=mStrLen) tmpch
    real(kind=fPrec) tmpflt
    integer nSplit, iElem, j, tmpi1, tmpi2, tmpi3
    logical spErr, tmpl, lfound
  
    call chr_split(inLine, lnSplit, nSplit, spErr)
    if(spErr) then
  