anchor.diff:    call alloc(elens_radial_jguess  ,           nelens_new,              -1, 'elens_radial_jguess'  )
   end subroutine elens_expand_arrays_lenses
   
  +subroutine elens_expand_arrays_rad_profiles(nelens_profiles_new)
  +  implicit none
  +  integer, intent(in) :: nelens_profiles_new
  +  call alloc(elens_radial_filename,  mFileName,     nelens_profiles_new,  " ", 'elens_radial_filename' )
  +  call alloc(elens_radial_profile_R,         0,     nelens_profiles_new, zero, 'elens_radial_profile_R', &
  +                                             0,     nelens_radial_profiles )
  +  call alloc(elens_radial_profile_J,         0,     nelens_profiles_new, zero, 'elens_radial_profile_J', &
  +                                             0,     nelens_radial_profiles )
  +  call alloc(elens_radial_profile_nPoints,          nelens_profiles_new,    0, 'elens_radial_profile_nPoints' )
  +end subroutine elens_expand_arrays_rad_profiles
  +
   ! ================================================================================================ !
   !  Parse Elens Line
   !  Rewritten from code from DATEN by VKBO
  
positive.code:   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
  
    ! variables to save elens parameters for tracking etc.
    integer, allocatable, save          :: elens_type(:)            ! integer for elens type (nelens)
  
negative.code:                                                                   !    i.e. linear interpolation
    integer, allocatable, save          :: elens_radial_jguess(:)   ! bin for guessed search (nelens)
    ! - file handling and data storage:
    character(len=:), allocatable, save :: elens_radial_filename(:)        ! names (nelens_radial_profiles)
    real(kind=fPrec), allocatable, save :: elens_radial_profile_R(:,:)     ! [mm] (0:elens_radial_dim,nelens_radial_profiles)
    real(kind=fPrec), allocatable, save :: elens_radial_profile_J(:,:)     ! [A]  (0:elens_radial_dim,nelens_radial_profiles)
    integer, allocatable, save          :: elens_radial_profile_nPoints(:) ! number of points in current radial profile (nelens_radial_profiles)
  
  contains
  
  