At the end of the Moist numerics a series of Saturation related code are executed: - GEOS_QsatICE - GEOS_DQSAT - GEOS_QsatLQU Original code is in file [@GMAO_Shared/GEOS_Shared/GEOS_Utilities.F90](https://github.com/GEOS-ESM/GMAO_Shared/blob/0ab88e82398072426d573bc5cb484939cc118837/GEOS_Shared/GEOS_Utilities.F90#L31) This should be saved aside of the normal Moist code since it's a shared code. Document heavily. Parent: https://github.com/GEOS-ESM/SMT-Nebulae/issues/57 --- - [ ] Port `QSATICEX` to `pyMoist` - [ ] Port `QSATLQU` to `pyMoist` - [ ] Port `DQSAT` to `pyMoist`