anchor.diff:    end if
   
     ! Parse files with radial profiles
  -   do j=1,melens_radial_profiles
  +   do j=1,nelens_radial_profiles
       inquire(file=elens_radial_filename(j), exist=exist)
       if(.not. exist) then
         write(lerr,"(a)") "ELENS> ERROR Problems with file with radial profile: "//trim(elens_radial_filename(j))
  
positive.code: 
    ! set current density at r(0) as at r(1), to properly compute the cumulative curve
    elens_radial_profile_J(0,ifile)=elens_radial_profile_J(1,ifile)
    
    write(lout,"(a,i0,a)") "ELENS> ...acquired ",elens_radial_profile_nPoints(ifile)," points."
    if(st_quiet < 2) then
      ! Echo parsed data (unless told to be quiet!)
  
negative.code:   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
  