Fix for species_configs file read from CSV

This version fixes an error that occurred if either CH4 or Aerosol-cloud interactions were missing from the list of species defined and the list of species configs in f.fill_species_configs() was read in from a CSV file.

OHC restart fix and py3.7 dropped

An inconsistency with ocean heat content before and after a restart has now been fixed. It was not affecting results since OHC is a post-processed quantity, only reporting.

One module dependency was failing tests with python 3.7, so support for python 3.7 has been dropped.

Zenodo DOI patch

The move to the new Zenodo API made pooch fail in the fill_from_rcmip() convenience function. The path to the emissions files from RCMIP have been updated to directly link to URLs rather than DOIs.

xarray fix and update to documentation

xarray 2023.9.0 introduced a change that broke FaIR with python 3.9+. Following user feedback, documentation has been updated, including, most importantly, an example using the latest available calibration. The positive semi-definite RuntimeWarnings should have also disappeared when running stochastically, or at least can be disabled.

Support for python 3.6 is removed due to dependency issues, and 3.11 added.

v2.1 + documentation

v2.1.0 is an overhaul of FaIR, and includes the following:

• object-oriented design
• interactive methane lifetime
• optional internal variability
• mass parallelisation

A new interface, and new features

v2.1.0 is an overhaul of FaIR, and includes the following:

• object-oriented design
• interactive methane lifetime
• optional internal variability

Add module interface for the SSPs

For a number of years, the RCP scenarios have been easily importable into FaIR using something like

from fair.RCPs import rcp45
from fair.forward import fair_scm

conc, forc, temp = fair_scm(rcp45.Emissions.emissions)

You can now do the same with the SSP scenarios from CMIP6:

from fair.SSPs import ssp119
from fair.forward import fair_scm

conc, forc, temp = fair_scm(ssp119.Emissions.emissions)

The following are available: ssp119, ssp126, ssp245, ssp370, ssp370_lowNTCF, ssp434, ssp460, ssp534, ssp585

More information is available from O'Neill et al. (2016) and Riahi et al. (2017), and the IPCC's Sixth Assessment Working Group 1 report. The emissions scenarios themselves were compiled by Zeb Nicholls and Jared Lewis and available at https://doi.org/10.5281/zenodo.4589756. Please cite Nicholls et al. (2020) if you use them,

The SSPs are supplied as a convenience as the actual SSP scenarios run for AR6 were done a little differently, but will need a whole configuration change to run as it was there, so it should be noted that results will probably not agree exactly with AR6 or observed history. This is just for people who want to easily get at the scenarios without digging too much into the AR6 repository.

Restarts for Geoffroy temperature function and GIR carbon cycle

An old feature given a new lease of life.

Updates to ozone forcing

This micro version introduces an updated treatment to the ozone forcing, which allows for a dependence on N2O concentrations, and a temperature-dependent ozone feedback (this was previously included with the stevenson treatment, but has now been simplified to be a linear function of temperature). The new treatment follows CMIP6 and AerChemMIP experiments:

Thornhill et al., 2021: https://acp.copernicus.org/articles/21/853/2021/acp-21-853-2021.pdf
Thornhill et al., accepted: https://acp.copernicus.org/preprints/acp-2019-1207/
Skeie et al., 2020: https://www.nature.com/articles/s41612-020-00131-0

Default coefficients that link the emissions to the CMIP6 model forcing in Skeie et al. are provided. In fair_scm() use tropO3_forcing=thornhill-skeie to use the new treatment.

Ozone forcing is not split into tropospheric and stratospheric components. For ease of communication with existing versions of FaIR, we maintain the forcing column indices, but all ozone forcing is put into the tropospheric bin (index 4 in 13-species FaIR, or index 31 in 45-species FaIR).

45-species FaIR

This version expands the multispecies version of FaIR to report 45 (previously 41) forcing species. Some fixes for running concentration driven runs with the 45-species FaIR have been implemented as well as including specified forcing and Geoffroy temperature options for inverse FaIR. scmdata interface has been removed.