Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Steady state emissions from concentrations tool #34

Merged
merged 5 commits into from
Aug 1, 2018

Conversation

chrisroadmap
Copy link
Collaborator

This goes some way to helping with #28 as it provides a simple way to calculate natural emissions of CH4 and N2O for varying lifetime runs.

Example for a pulse emission of methane:

import numpy as np
from fair.forward import fair_scm
from fair.constants import lifetime
from fair.tools import steady

# get steady state natural emissions for a 4-year methane lifetime and 
# preindustrial concentration of 722ppb
steady_ch4_4yr = steady.emissions(species="CH4", lifetime=4, C=722.)
steady_n2o     = steady.emissions(species="N2O", C=275.)

# load up default lifetimes but override methane to 4 years
lt = lifetime.aslist
lt[1] = 4.0

# 10 Gt methane pulse
emissions = np.zeros((100, 40))
emissions[0,3] = 10000.

# turn off forcing for everything except methane
scale = np.zeros(13)
scale[1] = 1.0
C_pi = np.zeros((31))
C_pi[0:3] = [278., 722., 275.]

# run FAIR
C,F,T = fair_scm(emissions, lifetimes=lt, natural = [steady_ch4_4yr, steady_n2o], F_volcanic=0., F_solar=0., scale=scale, fixPre1850RCP=False)

@chrisroadmap chrisroadmap changed the title Steady state emissions from concentrations to.l Steady state emissions from concentrations tool Jul 10, 2018
@chrisroadmap chrisroadmap merged commit e67b6bc into master Aug 1, 2018
@chrisroadmap chrisroadmap deleted the equilibriate_natural branch August 14, 2018 12:43
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant