Python package to perform calculations with the FaIR simple climate model
Project description
FaIR
Finite Amplitude Impulse-Response simple climate-carbon-cycle model
Installation
Make sure you have Python 3.6+ and pip installed
From terminal/command prompt pip install fair
Usage
FaIR takes emissions of greenhouse gases, aerosol and ozone precursors, and converts these into greenhouse gas concentrations, radiative forcing and temperature change.
There are two ways to run FaIR:
Carbon dioxide emissions only with all other radiative forcings specified externally (specify useMultigas=False in the call to fair_scm);
All species included in the RCP emissions datasets, with, optionally, solar and volcanic forcing still specified externally. For convenience, the RCP datasets are provided in the RCP subdirectory and can be imported:
from fair.forward import fair_scm from fair.RCPs import rcp85 emissions = rcp85.Emissions.emissions C,F,T = fair_scm(emissions=emissions)
The main engine of the model is the fair_scm function in forward.py. This function can be imported into a Python script or iPython session. The most important keyword to fair_scm is the emissions. This should be either a (nt, 40) numpy array (in multigas mode) or (nt,) numpy array (in CO2 only mode), where nt is the number of model timesteps. The outputs are a tuple of (C, F, T) arrays which are GHG concentrations ((nt, 31) in multigas mode, (nt,) in CO2-only mode), forcing ((nt, 13) or (nt,)) and temperature change (nt,). The index numbers corresponding to each species will be given in tables 1 to 3 of the revised version of the Smith et al. paper reference below (we hope to make this object-oriented in the future). For now, note that the input emissions follow the ordering of the RCP datasets, which are included under fair/RCPs, and the GHG concentrations output are in the same order, except that we don’t output the year, only use one column for total CO2, and the short-lived species (input indices 5 to 11 inclusive) are not included, reducing the number of columns from 40 to 31. In multigas mode the forcing output indices are:
CO2
CH4
N2O
Minor GHGs (CFCs, HFCs etc)
Tropospheric ozone
Stratospheric ozone
Stratospheric water vapour from methane oxidation
Contrails
Aerosols
Black carbon on snow
Land use
Volcanic
Solar
For further information, see the example ipython notebook contained in the GitHub repo at https://github.com/OMS-NetZero/FAIR.
References:
Smith, C. J., Forster, P. M., Allen, M., Leach, N., Millar, R. J., Passerello, G. A., and Regayre, L. A.: FAIR v1.3: A simple emissions-based impulse response and carbon cycle model, Geosci. Model Dev., https://doi.org/10.5194/gmd-11-2273-2018, 2018.
Millar, R. J., Nicholls, Z. R., Friedlingstein, P., and Allen, M. R.: A modified impulse-response representation of the global near-surface air temperature and atmospheric concentration response to carbon dioxide emissions, Atmos. Chem. Phys., 17, 7213-7228, https://doi.org/10.5194/acp-17-7213-2017, 2017.
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
