The documentation should include examples of how to use various features of the package.
For example:

• How to deal with artifacts when plotting contour lines on stereographic plots (e.g. sea-ice concentration)
• How to do seasonal decomposition of time series data
• How to use the different accessors (xarray, polar plots, etc.)

https://github.com/markusritschel/my-code-base/blob/cc61abb06f2587e45cf50e60d3a25a76b8d5ea7a/my_code_base/esd/utils.py#L130-L140

Use pyproject.toml for package configuration and install:
https://packaging.python.org/en/latest/guides/writing-pyproject-toml/
https://setuptools.pypa.io/en/latest/userguide/pyproject_config.html

Read more about packaging:
https://packaging.python.org/en/latest/flow/

Provide a decorator for plt.savefig that feeds values to the metadata argument. For example:

• the current git commit hash
• the script producing the plot, plus the line of the function or the plt.savefig statement
• version of the code/project

Why do we integrate over autocorr[i] * autocorr[i+1]?

            τ += dt * (autocorr[i] + autocorr[i + 1]) / 2   # TODO: Why is the integral not just over autocorr[i]?

https://github.com/markusritschel/my-code-base/blob/4e779beca5640f4f3bbf720b8170e4258a214858/src/my_code_base/stats/timeseries.py#L268