Automatic order selection

[GreatEmerald]

We could select the harmonic order automatically with an information criterion: fit models with different orders, then compare the model fit and choose one that is simplest. This would slow down runtime but might help against overfitting.