anchor.diff:      cheby_maxOrder(ifile)=max(ii,jj,cheby_maxOrder(ifile))
       if ( cheby_maxOrder(ifile) > cheby_maxNterms-1 ) then
         write(lerr,"(2(a,i0))") "CHEBY> ERROR order of map too high: ",cheby_maxOrder(ifile), &
  -                                            " >",cheby_maxNterms-1
  +                                            " > ",cheby_maxNterms-1
         write(lerr,"(a)")       "CHEBY>       please increase cheby_maxNterms parameter"
         goto 30
       end if
  
positive.code: 
  
  ! ================================================================================================ !
  !  Last modified: 2019-02-25
  !  Rewritten by VKBO, June 2018
  !  Read file with coefficients for chebyshev polynomials
  !  ifile is index of file in table of chebyshev files
  
negative.code:     ! Check if profile has already been requested:
      chIdx = -1
      do tmpi1=1,ncheby_tables
        if(trim(tmpch) == trim(cheby_filename(tmpi1))) then
          cheby_itable(icheby(iElem)) = tmpi1
          chIdx = tmpi1
          exit
  