Updated seasonal model code

- Introduced a new function `momlevel.trend.deseason`
  that is based on the existing `seasonal_model` code
- New code supports both dask arrays and in-memory
  numpy arrays
- Options to return the fitted model, residuals from
  the model, and the coefficients of the model

Towards a faster seasonal model routine

Updated seasonal model code that works with coeffs

Updated seasonal model