anchor.diff:              return
            end if
   
  +            ! Check that the second dimensions of the SplineCoefs table and the Coefs table match.  If they don't, it may indicate
  +            ! that the number of columns in each table in the AF input file are different.  This will result in array bounds issues
  +            ! in the CubicSplineInitM routine.
  +         if ( size(p%Table(iTable)%SplineCoefs, DIM=2) /= size(p%Table(iTable)%Coefs, DIM=2) ) then
  +            call SetErrStat ( ErrID_Fatal, 'Number of columns in the SplineCoefs table and Coefs tables do not match in size.'// &
  +                  ' Check that all tables in airfoil input file "'//TRIM( InitInput%FileName )//'" have the same number of columns.', &
  +                  ErrStat, ErrMsg, RoutineName ) 
  +            call Cleanup()
  +            return
  +         end if
               
               ! Compute the one set of coefficients of the piecewise polynomials for the irregularly-spaced data.
               ! Unlike the 2-D interpolation in which we use diffent knots for each airfoil coefficient, we can do
  
positive.code:                CALL ParseVarWDefault ( FileInfo, CurLine, 'filtCutOff', p%Table(iTable)%UA_BL%filtCutOff, 0.5_ReKi, ErrStat2, ErrMsg2, UnEc )
                    CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
                    
                 p%ColUAf = NumCoefs + 1 ! column for f_st
                 NumCoefs = p%ColUAf + 1 ! precompute f_st and cl_fs for the HGM model
                 
                 IF (ErrStat >= AbortErrLev) THEN
                    CALL Cleanup()
  
negative.code:       ENDIF
  
        DO iTable=1,p%NumTabs
           NumCoefsTab = NumCoefsIn  ! Reset this counter for each table
  
           CALL ParseVar ( FileInfo, CurLine, 'Re', p%Table(iTable)%Re, ErrStat2, ErrMsg2, UnEc )
              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  