def METHOD_NAME( profile_2d, dx1, dx2, deg1, deg2 ):
    """
    Compute error bound for spline approximation of 2D analytical profile.
    Parameters
    ----------
    profile_2d : 2D analytical profile
        Must provide 'max_norm( n1,n2 )' method to compute max norm of its
        mixed derivative of degree (n1,n2) over domain.
    dx1 : float
        Grid spacing along 1st dimension.
    dx2 : float
        Grid spacing along 2nd dimension.
    deg1 : int
        Spline degree along 1st dimension.
    deg2 : int
        Spline degree along 2nd dimension.
    Result
    ------
    max_error : float
        Error bound: max-norm of error over domain should be smaller than this.
    """
    # Max norm of highest partial derivatives in x1 and x2 of analytical profile
    max_norm1 = profile_2d.max_norm( deg1+1, 0      )
    max_norm2 = profile_2d.max_norm( 0     , deg2+1 )
    # Error bound on function value
    max_error = f_tihomirov_error_bound( dx1, deg1, max_norm1 ) \
              + f_tihomirov_error_bound( dx2, deg2, max_norm2 )
    # Empirical correction: for linear interpolation increase estimate by 5%
    if (deg1 == 1 or deg2 == 1):
        max_error = 1.05 * max_error
    return max_error