anchor.diff:    call f_close(fUnit)
     write(lout,"(a,i0,a)") "ELENS> ...acquired ",elens_radial_profile_nPoints(ifile)," points."
   
  +  ! check array of x-values is sensible
  +  if(.not.checkArray(elens_radial_profile_R(1:elens_radial_profile_nPoints(ifile),ifile),elens_radial_profile_nPoints(ifile))) then
  +    iErr=3
  +    write(lerr,"(a)") "ELENS> ERROR radial profile has problem with values of radius."
  +    goto 30
  +  end if
  +  
     if(st_quiet < 2) then
       ! Echo parsed data (unless told to be quiet!)
       write(lout,"(a,i0)") "ELENS> Radial profile as from file "//&
  
positive.code:   end if
    call chr_cast(lnSplit(1),tmpR,spErr)
    call chr_cast(lnSplit(2),tmpJ,spErr)
    if(tmpR<=elens_radial_profile_R(ii,ifile)) then
      iErr = 1
      write(lerr,"(a,i0)") "ELENS> ERROR radius not in increasing order at ii=",ii
      goto 30
    end if
    if(tmpJ>=0.0) then
      ii=ii+1
      if(ii>elens_radial_dim) then
  
negative.code:   use mod_settings
    use string_tools
    use mod_units
    use utils, only: checkArray
  
    implicit none
  
  