open-atmos packages 
Before going into specific descriptions of our proposed solutions, let us introduce open-source Pythonic packages that will serve as examples and show-cases. All of them are developed and maintained by open-atmos contributors .
Packages started over 5 years ago, hence, serve as good examples for discussion about maintainability.
PySDM 
PySDM is a package for simulating the dynamics of a population of particles undergoing diffusional and collisional growth (and breakage). The package features a Pythonic high-performance (multithreaded CPU & CUDA GPU) implementation of the Super-Droplet Method (SDM) Monte-Carlo algorithm for representing collisional growth Shima et al. [2009], hence the name.
The animation from the landing page PySDM Documentation
PyMPDATA 
PyMPDATA is a Numba-accelerated multithreaded Pythonic implementation of the MPDATA algorithm of Smolarkiewicz [1984] used in geophysical fluid dynamics for solving convection-diffusion PDEs. PyMPDATA supports integration in 1D, 2D, and 3D structured meshes with optional coordinate transformations.
More information and examples can be found in PyMPDATA Documentation
open-atmos-jupyter-utils 
open-atmos-jupyter-utils is a Python package providing Jupyter notebook utility routines for:
presenting Matplotlib plots as either SVG vector graphics or animated GIFs, embedding them within the notebooks, and rendering correctly in GitHub’s Rich Jupyter Notebook diffs
save-as buttons below each figure (triggering Google-Drive downloads on Colab)
execution of unmodified notebook code for automated testing (e.g., within pytest fixtures)
pip-installation of external packages on Colab safeguarded against alterations of Google-shipped packages
Can be installed via pip
pip install open-atmos-jupyter-utils