def test_scenfile():
datadir = os.path.join(os.path.dirname(__file__), 'rcp45/')
# Purpose of this test is to determine whether the SCEN file for RCP4.5
# which does not include CFCs, years before 2000, or emissions from every
# year from 2000 to 2500, equals the emissions file from RCP4.5
# after reconstruction.
# The .SCEN and .XLS files at http://www.pik-potsdam.de/~mmalte/rcps
# sometimes differ in the 4th decimal place. Thus we allow a tolerance of
# 0.0002 in addition to machine error in this instance.
E1 = magicc.scen_open(datadir + 'RCP45.SCEN')
assert np.allclose(E1, rcp45.Emissions.emissions, rtol=1e-8, atol=2e-4)
E2 = magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs=False)
assert np.allclose(E2,
rcp45.Emissions.emissions[:, :24],
rtol=1e-8,
atol=2e-4)
E3 = magicc.scen_open(datadir + 'RCP45.SCEN', startyear=1950)
assert np.allclose(E3,
rcp45.Emissions.emissions[185:, :],
rtol=1e-8,
atol=2e-4)
# Test `_import_emis_file`
assert magicc._import_emis_file('rcp26') == rcp3pd.Emissions
assert magicc._import_emis_file('rcp6') == rcp6.Emissions
assert magicc._import_emis_file('rcp85') == rcp85.Emissions
with pytest.raises(ValueError):
magicc._import_emis_file('rcp19')
test_files = os.path.join(os.path.dirname(__file__), "scenfiles")
scenfile_2000 = os.path.join(test_files, "WORLD_ONLY.SCEN")
scenfile_2010 = os.path.join(test_files, "WORLD_ONLY_2010.SCEN")
# Test CFCs inclusion.
with pytest.raises(ValueError):
magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs=np.zeros(0))
E4 = magicc.scen_open(scenfile_2000,
startyear=2000,
include_cfcs=np.ones((51, 16)))
assert E4[0, -1] == 1
with pytest.raises(ValueError):
magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs="foo")
# Test filling of history and harmonisation.
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2010)
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2000, harmonise=1950)
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2000, harmonise=2060)
E5 = magicc.scen_open(scenfile_2000, harmonise=2010)
assert E5[0, 1] == rcp45.Emissions.co2_fossil[0]
def test_scenfile():
datadir = os.path.join(os.path.dirname(__file__), 'rcp45/')
# Purpose of this test is to determine whether the SCEN file for RCP4.5
# which does not include CFCs, years before 2000, or emissions from every
# year from 2000 to 2500, equals the emissions file from RCP4.5
# after reconstruction.
# The .SCEN and .XLS files at http://www.pik-potsdam.de/~mmalte/rcps
# sometimes differ in the 4th decimal place. Thus we allow a tolerance of
# 0.0002 in addition to machine error in this instance.
E1 = magicc.scen_open(datadir + 'RCP45.SCEN')
assert np.allclose(E1, rcp45.Emissions.emissions, rtol=1e-8, atol=2e-4)
E2 = magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs=False)
assert np.allclose(E2, rcp45.Emissions.emissions[:,:24], rtol=1e-8,
atol=2e-4)
E3 = magicc.scen_open(datadir + 'RCP45.SCEN', startyear=1950)
assert np.allclose(E3, rcp45.Emissions.emissions[185:,:], rtol=1e-8,
atol=2e-4)
# Test `_import_emis_file`
assert magicc._import_emis_file('rcp26') == rcp3pd.Emissions
assert magicc._import_emis_file('rcp6') == rcp6.Emissions
assert magicc._import_emis_file('rcp85') == rcp85.Emissions
with pytest.raises(ValueError):
magicc._import_emis_file('rcp19')
test_files = os.path.join(os.path.dirname(__file__), "scenfiles")
scenfile_2000 = os.path.join(test_files, "WORLD_ONLY.SCEN")
scenfile_2010 = os.path.join(test_files, "WORLD_ONLY_2010.SCEN")
# Test CFCs inclusion.
with pytest.raises(ValueError):
magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs=np.zeros(0))
E4 = magicc.scen_open(
scenfile_2000, startyear=2000, include_cfcs=np.ones((51, 16)))
assert E4[0, -1] == 1
with pytest.raises(ValueError):
magicc.scen_open(datadir + 'RCP45.SCEN', include_cfcs="foo")
# Test filling of history and harmonisation.
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2010)
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2000, harmonise=1950)
with pytest.raises(ValueError):
magicc.scen_open(scenfile_2000, harmonise=2060)
E5 = magicc.scen_open(scenfile_2000, harmonise=2010)
assert E5[0, 1] == rcp45.Emissions.co2_fossil[0]