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:   call chr_cast(lnSplit(1),tmpR,spErr)
    call chr_cast(lnSplit(2),tmpJ,spErr)
    ii=ii+1
    if(ii>elens_radial_dim) then
      iErr = 2
      write(lerr,"(a,i0,a,i0)") "ELENS> ERROR too many points in radial profile: ",ii, &
           ". Max is ",elens_radial_dim
      goto 30
    end if
    elens_radial_profile_nPoints(ifile) = ii
    elens_radial_profile_R(ii,ifile) = tmpR
  
negative.code:       tmpch = trim(lnSplit(8))
      
        ! Check if profile has already been requested:
        lfound = .false.
        do tmpi1=1,nelens_radial_profiles
          if(tmpch == elens_radial_filename(tmpi1)) then
            elens_iRadial(ielens(iElem)) = tmpi1
            lfound = .true.
            exit
          end if
        end do
        if(.not.lfound) then
          ! Unsuccessful search
          nelens_radial_profiles = nelens_radial_profiles+1
          call elens_expand_arrays_rad_profiles(nelens_radial_profiles)
          elens_iRadial(ielens(iElem)) = nelens_radial_profiles
          elens_radial_filename(nelens_radial_profiles) = tmpch
        end if
      end if
  
  