From 89200e43cbc9db7848797016dca8836e199c878d Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Mon, 20 Nov 2023 12:25:29 +0100 Subject: [PATCH 01/11] Add new stellar mass - BH mass plot --- colibre/auto_plotter/black_holes.yml | 33 ++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) diff --git a/colibre/auto_plotter/black_holes.yml b/colibre/auto_plotter/black_holes.yml index 008c8bbd..e4a5c3b0 100644 --- a/colibre/auto_plotter/black_holes.yml +++ b/colibre/auto_plotter/black_holes.yml @@ -67,6 +67,39 @@ stellar_mass_black_hole_mass_50: - filename: GalaxyStellarMassBlackHoleMass/Graham2023_S.hdf5 - filename: GalaxyStellarMassBlackHoleMass/McConnell2013_Fit.hdf5 +stellar_mass_black_hole_mass_binned_50: + type: "scatter" + legend_loc: "upper left" + x: + quantity: "apertures.mass_star_50_kpc" + units: Solar_Mass + start: 1e6 + end: 1e12 + y: + quantity: "black_hole_masses.max" + units: Solar_Mass + start: 1e3 + end: 1e11 + median: + plot: true + log: true + adaptive: true + number_of_bins: 30 + start: + value: 1e6 + units: Solar_Mass + end: + value: 1e12 + units: Solar_Mass + metadata: + title: Stellar Mass-Black Hole Mass relation (50 kpc Stellar Mass) + caption: SMBHM relation. The binned observational data are obtained by combining BH mass measurements for galaxies of different morphologies (Graham 2023) and passive vs. active galaxies (Terrazas et al. 2017) with dependencies on stellar mass of morphological types (Moffett et al. 2016, from the GAMA survey) and quenched fractions (Gilbank et al. 2010, SDSS). + section: Black Holes + observational_data: + - filename: GalaxyStellarMassBlackHoleMassBinned/Terrazas2017_binned.hdf5 + - filename: GalaxyStellarMassBlackHoleMassBinned/Graham2023_binned.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/McConnell2013_Fit.hdf5 + stellar_mass_black_hole_mass_passive_30: type: "scatter" select_structure_type: 10 From c45245328f875ca7ab92129adae07a20e2a7a39d Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Thu, 21 Mar 2024 10:47:38 +0000 Subject: [PATCH 02/11] Add new data to BH accretion evolution plot --- colibre/scripts/bh_accretion_evolution.py | 202 +++++++++++++++++++++- 1 file changed, 194 insertions(+), 8 deletions(-) diff --git a/colibre/scripts/bh_accretion_evolution.py b/colibre/scripts/bh_accretion_evolution.py index f1218787..d24db32f 100644 --- a/colibre/scripts/bh_accretion_evolution.py +++ b/colibre/scripts/bh_accretion_evolution.py @@ -16,6 +16,7 @@ from swiftpipeline.argumentparser import ScriptArgumentParser from velociraptor.observations import load_observation +from velociraptor.observations import load_observations arguments = ScriptArgumentParser( description="Creates a BH accretion rate density evolution plot, with added observational data." @@ -75,15 +76,197 @@ ) simulation_labels.append(name) -# Observational data plotting +path_to_obs_data = f"{arguments.config.config_directory}/{arguments.config.observational_data_directory}" -observational_data = glob.glob( - f"{arguments.config.config_directory}/{arguments.config.observational_data_directory}/data/BlackHoleAccretionHistory/*.hdf5" +observation_lines = [] +observation_labels = [] + +Aird2015 = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Aird2015.hdf5"] +)[0] +observation_lines.append( + ax.plot( + Aird2015.x.value, + Aird2015.y.value, + color='brown', + zorder=-10000, + alpha=0.6, + linewidth=1, + )[0] +) +observation_labels.append(Aird2015.citation) + +Annana2020_low = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Annana2020_low.hdf5"] +)[0] +observation_lines.append( + ax.plot( + Annana2020_low.x.value, + Annana2020_low.y.value, + label=Annana2020_low.citation, + zorder=-10000, + color='orange', + alpha=0.6, + linewidth=1, + )[0] +) +observation_labels.append(Annana2020_low.citation) + +Annana2020_high = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Annana2020_high.hdf5"] +)[0] +observation_lines.append( + ax.plot( + Annana2020_high.x.value, + Annana2020_high.y.value, + label=Annana2020_high.citation, + color='olive', + zorder=-10000, + alpha=0.6, + linewidth=1, + )[0] +) +observation_labels.append(Annana2020_high.citation) + +Pouliasis2024_low = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Pouliasis2024_low.hdf5"] +)[0] +observation_lines.append( + ax.fill_between( + Pouliasis2024_low.x.value, + Pouliasis2024_low.y.value - Pouliasis2024_low.y_scatter[0].value, + Pouliasis2024_low.y.value + Pouliasis2024_low.y_scatter[1].value, + color='magenta', + label=Pouliasis2024_low.citation, + zorder=-10000, + alpha=0.2, + ) ) +observation_labels.append(Pouliasis2024_low.citation) -for index, observation in enumerate(observational_data): - obs = load_observation(observation) - obs.plot_on_axes(ax) +Pouliasis2024_high = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Pouliasis2024_high.hdf5"] +)[0] +observation_lines.append( + ax.fill_between( + Pouliasis2024_high.x.value, + Pouliasis2024_high.y.value - Pouliasis2024_high.y_scatter[0].value, + Pouliasis2024_high.y.value + Pouliasis2024_high.y_scatter[1].value, + color='teal', + label=Pouliasis2024_high.citation, + zorder=-10000, + alpha=0.2, + ) +) +observation_labels.append(Pouliasis2024_high.citation) + +ax.set_xlabel("Redshift $z$") +ax.set_ylabel(r"BH Accretion Rate Density [M$_\odot$ yr$^{-1}$ Mpc$^{-3}$]") + +Delvacchio2014 = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Delvacchio2014.hdf5"] +)[0] +observation_lines.append( + ax.errorbar( + Delvacchio2014.x.value, + Delvacchio2014.y.value, + xerr=None if Delvacchio2014.x_scatter is None else Delvacchio2014.x_scatter.value, + yerr=None if Delvacchio2014.y_scatter is None else Delvacchio2014.y_scatter.value, + label=Delvacchio2014.citation, + linestyle="none", + marker="o", + color="maroon", + elinewidth=0.75, + markeredgecolor="none", + markersize=3, + zorder=-10, + )[0] +) +observation_labels.append(Delvacchio2014.citation) + +DSilva2023 = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/DSilva2023.hdf5"] +)[0] +observation_lines.append( + ax.errorbar( + DSilva2023.x.value, + DSilva2023.y.value, + xerr=None if DSilva2023.x_scatter is None else DSilva2023.x_scatter.value, + yerr=None if DSilva2023.y_scatter is None else DSilva2023.y_scatter.value, + label=DSilva2023.citation, + linestyle="none", + marker="o", + color="darkblue", + elinewidth=0.75, + markeredgecolor="none", + markersize=3, + zorder=-10, + )[0] +) +observation_labels.append(DSilva2023.citation) + +Yang2023 = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Yang2023.hdf5"] +)[0] +observation_lines.append( + ax.errorbar( + Yang2023.x.value, + Yang2023.y.value, + xerr=None if Yang2023.x_scatter is None else Yang2023.x_scatter.value, + yerr=None if Yang2023.y_scatter is None else Yang2023.y_scatter.value, + label=Yang2023.citation, + linestyle="none", + marker="o", + color="goldenrod", + elinewidth=0.75, + markeredgecolor="none", + markersize=3, + zorder=-10, + ) +) +observation_labels.append(Yang2023.citation) + +Shen2020_low = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Shen2020_low.hdf5"] +)[0] +observation_lines.append( + ax.errorbar( + Shen2020_low.x.value, + Shen2020_low.y.value, + xerr=None if Shen2020_low.x_scatter is None else Shen2020_low.x_scatter.value, + yerr=None if Shen2020_low.y_scatter is None else Shen2020_low.y_scatter.value, + label=Shen2020_low.citation, + linestyle="none", + marker="o", + color="lime", + elinewidth=0.5, + markeredgecolor="none", + markersize=2, + zorder=-10, + ) +) +observation_labels.append(Shen2020_low.citation) + +Shen2020_high = load_observations( + [f"{path_to_obs_data}/data/BlackHoleAccretionHistory/Shen2020_high.hdf5"] +)[0] +observation_lines.append( + ax.errorbar( + Shen2020_high.x.value, + Shen2020_high.y.value, + xerr=None if Shen2020_high.x_scatter is None else Shen2020_high.x_scatter.value, + yerr=None if Shen2020_high.y_scatter is None else Shen2020_high.y_scatter.value, + label=Shen2020_high.citation, + linestyle="none", + marker="o", + color="coral", + elinewidth=0.5, + markeredgecolor="none", + markersize=2, + zorder=-10, + ) +) +observation_labels.append(Shen2020_high.citation) ax.set_xlabel("Redshift $z$") ax.set_ylabel(r"BH Accretion Rate Density [M$_\odot$ yr$^{-1}$ Mpc$^{-3}$]") @@ -109,9 +292,12 @@ ax.tick_params(axis="x", which="minor", bottom=False) ax.set_xlim(1.02, 0.07) -ax.set_ylim(6e-7, 6e-4) +ax.set_ylim(1e-7, 6e-4) -observation_legend = ax.legend(markerfirst=True, loc="lower left") +#observation_legend = ax.legend(markerfirst=True, loc="lower left") +observation_legend = ax.legend( + observation_lines, observation_labels, markerfirst=True, loc="lower center", fontsize=4, ncol=1 +) ax.add_artist(observation_legend) From c665dca2aeaa8b2f143bc0b7d46d69c0b442e79d Mon Sep 17 00:00:00 2001 From: Filip Husko Date: Thu, 21 Mar 2024 10:50:05 +0000 Subject: [PATCH 03/11] Run formatting --- colibre/scripts/bh_accretion_evolution.py | 59 +++++++++++++---------- 1 file changed, 34 insertions(+), 25 deletions(-) diff --git a/colibre/scripts/bh_accretion_evolution.py b/colibre/scripts/bh_accretion_evolution.py index d24db32f..f6a85b35 100644 --- a/colibre/scripts/bh_accretion_evolution.py +++ b/colibre/scripts/bh_accretion_evolution.py @@ -88,7 +88,7 @@ ax.plot( Aird2015.x.value, Aird2015.y.value, - color='brown', + color="brown", zorder=-10000, alpha=0.6, linewidth=1, @@ -105,7 +105,7 @@ Annana2020_low.y.value, label=Annana2020_low.citation, zorder=-10000, - color='orange', + color="orange", alpha=0.6, linewidth=1, )[0] @@ -120,7 +120,7 @@ Annana2020_high.x.value, Annana2020_high.y.value, label=Annana2020_high.citation, - color='olive', + color="olive", zorder=-10000, alpha=0.6, linewidth=1, @@ -136,7 +136,7 @@ Pouliasis2024_low.x.value, Pouliasis2024_low.y.value - Pouliasis2024_low.y_scatter[0].value, Pouliasis2024_low.y.value + Pouliasis2024_low.y_scatter[1].value, - color='magenta', + color="magenta", label=Pouliasis2024_low.citation, zorder=-10000, alpha=0.2, @@ -152,7 +152,7 @@ Pouliasis2024_high.x.value, Pouliasis2024_high.y.value - Pouliasis2024_high.y_scatter[0].value, Pouliasis2024_high.y.value + Pouliasis2024_high.y_scatter[1].value, - color='teal', + color="teal", label=Pouliasis2024_high.citation, zorder=-10000, alpha=0.2, @@ -170,8 +170,12 @@ ax.errorbar( Delvacchio2014.x.value, Delvacchio2014.y.value, - xerr=None if Delvacchio2014.x_scatter is None else Delvacchio2014.x_scatter.value, - yerr=None if Delvacchio2014.y_scatter is None else Delvacchio2014.y_scatter.value, + xerr=None + if Delvacchio2014.x_scatter is None + else Delvacchio2014.x_scatter.value, + yerr=None + if Delvacchio2014.y_scatter is None + else Delvacchio2014.y_scatter.value, label=Delvacchio2014.citation, linestyle="none", marker="o", @@ -235,14 +239,14 @@ Shen2020_low.y.value, xerr=None if Shen2020_low.x_scatter is None else Shen2020_low.x_scatter.value, yerr=None if Shen2020_low.y_scatter is None else Shen2020_low.y_scatter.value, - label=Shen2020_low.citation, - linestyle="none", - marker="o", - color="lime", - elinewidth=0.5, - markeredgecolor="none", - markersize=2, - zorder=-10, + label=Shen2020_low.citation, + linestyle="none", + marker="o", + color="lime", + elinewidth=0.5, + markeredgecolor="none", + markersize=2, + zorder=-10, ) ) observation_labels.append(Shen2020_low.citation) @@ -256,14 +260,14 @@ Shen2020_high.y.value, xerr=None if Shen2020_high.x_scatter is None else Shen2020_high.x_scatter.value, yerr=None if Shen2020_high.y_scatter is None else Shen2020_high.y_scatter.value, - label=Shen2020_high.citation, - linestyle="none", - marker="o", - color="coral", - elinewidth=0.5, - markeredgecolor="none", - markersize=2, - zorder=-10, + label=Shen2020_high.citation, + linestyle="none", + marker="o", + color="coral", + elinewidth=0.5, + markeredgecolor="none", + markersize=2, + zorder=-10, ) ) observation_labels.append(Shen2020_high.citation) @@ -294,9 +298,14 @@ ax.set_xlim(1.02, 0.07) ax.set_ylim(1e-7, 6e-4) -#observation_legend = ax.legend(markerfirst=True, loc="lower left") +# observation_legend = ax.legend(markerfirst=True, loc="lower left") observation_legend = ax.legend( - observation_lines, observation_labels, markerfirst=True, loc="lower center", fontsize=4, ncol=1 + observation_lines, + observation_labels, + markerfirst=True, + loc="lower center", + fontsize=4, + ncol=1, ) ax.add_artist(observation_legend) From fbbebb589758e265f3cc54ef2ce41a506c15c2bd Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Thu, 21 Mar 2024 16:57:56 +0000 Subject: [PATCH 04/11] Remove duplicate line --- colibre/scripts/bh_accretion_evolution.py | 1 - 1 file changed, 1 deletion(-) diff --git a/colibre/scripts/bh_accretion_evolution.py b/colibre/scripts/bh_accretion_evolution.py index f6a85b35..df6f9660 100644 --- a/colibre/scripts/bh_accretion_evolution.py +++ b/colibre/scripts/bh_accretion_evolution.py @@ -15,7 +15,6 @@ from swiftpipeline.argumentparser import ScriptArgumentParser -from velociraptor.observations import load_observation from velociraptor.observations import load_observations arguments = ScriptArgumentParser( From 89e7f9b22e5f4a410382a196c3f94f486dfc63b5 Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Thu, 21 Mar 2024 16:58:55 +0000 Subject: [PATCH 05/11] Correct directory name for binned BH mass data Co-authored-by: Evgenii Chaikin --- colibre/auto_plotter/black_holes.yml | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/colibre/auto_plotter/black_holes.yml b/colibre/auto_plotter/black_holes.yml index e4a5c3b0..3a07045a 100644 --- a/colibre/auto_plotter/black_holes.yml +++ b/colibre/auto_plotter/black_holes.yml @@ -96,8 +96,8 @@ stellar_mass_black_hole_mass_binned_50: caption: SMBHM relation. The binned observational data are obtained by combining BH mass measurements for galaxies of different morphologies (Graham 2023) and passive vs. active galaxies (Terrazas et al. 2017) with dependencies on stellar mass of morphological types (Moffett et al. 2016, from the GAMA survey) and quenched fractions (Gilbank et al. 2010, SDSS). section: Black Holes observational_data: - - filename: GalaxyStellarMassBlackHoleMassBinned/Terrazas2017_binned.hdf5 - - filename: GalaxyStellarMassBlackHoleMassBinned/Graham2023_binned.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Terrazas2017_binned.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Graham2023_binned.hdf5 - filename: GalaxyStellarMassBlackHoleMass/McConnell2013_Fit.hdf5 stellar_mass_black_hole_mass_passive_30: From 1430587b476233b1de1675f02cea274d4af12de3 Mon Sep 17 00:00:00 2001 From: Filip Husko Date: Fri, 22 Mar 2024 12:36:41 +0000 Subject: [PATCH 06/11] update submodule --- observational_data | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/observational_data b/observational_data index 58ae6627..a4f9ff9a 160000 --- a/observational_data +++ b/observational_data @@ -1 +1 @@ -Subproject commit 58ae66274d7a2dab07e700d3d55536d8e08915c5 +Subproject commit a4f9ff9ac2d5287d224067a5c459a17bb5ebe6b8 From 425b46d6104caca83b722c0f5bf7669b8fdf66f7 Mon Sep 17 00:00:00 2001 From: Filip Husko Date: Fri, 22 Mar 2024 12:46:13 +0000 Subject: [PATCH 07/11] update submodule #2 --- observational_data | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/observational_data b/observational_data index a4f9ff9a..28a17cda 160000 --- a/observational_data +++ b/observational_data @@ -1 +1 @@ -Subproject commit a4f9ff9ac2d5287d224067a5c459a17bb5ebe6b8 +Subproject commit 28a17cdab54a417090710a17962c33ac00d7fefb From ab6cce05b19781250a54fdcdd1f79dfb4f31ef2b Mon Sep 17 00:00:00 2001 From: Filip Husko Date: Fri, 22 Mar 2024 13:03:36 +0000 Subject: [PATCH 08/11] update submodule #3 --- observational_data | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/observational_data b/observational_data index 28a17cda..0c425e8d 160000 --- a/observational_data +++ b/observational_data @@ -1 +1 @@ -Subproject commit 28a17cdab54a417090710a17962c33ac00d7fefb +Subproject commit 0c425e8dd38cbe60b135419697ea1b1e4b469f1f From ff60e6a38bdcfb3b4026ad9ff92ee4de1c420d45 Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Wed, 1 May 2024 15:34:43 +0100 Subject: [PATCH 09/11] Add high-redshift data for BH-stellar mass relation --- colibre/auto_plotter/black_holes.yml | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/colibre/auto_plotter/black_holes.yml b/colibre/auto_plotter/black_holes.yml index 3a07045a..5bd6c2c9 100644 --- a/colibre/auto_plotter/black_holes.yml +++ b/colibre/auto_plotter/black_holes.yml @@ -32,6 +32,9 @@ stellar_mass_black_hole_mass_30: - filename: GalaxyStellarMassBlackHoleMass/Graham2023_E.hdf5 - filename: GalaxyStellarMassBlackHoleMass/Graham2023_S.hdf5 - filename: GalaxyStellarMassBlackHoleMass/McConnell2013_Fit.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Schramm2013.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Ding2020.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Tanaka2024.hdf5 stellar_mass_black_hole_mass_50: type: "scatter" @@ -66,6 +69,9 @@ stellar_mass_black_hole_mass_50: - filename: GalaxyStellarMassBlackHoleMass/Graham2023_E.hdf5 - filename: GalaxyStellarMassBlackHoleMass/Graham2023_S.hdf5 - filename: GalaxyStellarMassBlackHoleMass/McConnell2013_Fit.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Schramm2013.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Ding2020.hdf5 + - filename: GalaxyStellarMassBlackHoleMass/Tanaka2024.hdf5 stellar_mass_black_hole_mass_binned_50: type: "scatter" From 6581e752faf00aab3c2067d8b5c152de5b9851db Mon Sep 17 00:00:00 2001 From: FHusko <53567522+FHusko@users.noreply.github.com> Date: Wed, 1 May 2024 15:40:59 +0100 Subject: [PATCH 10/11] Add new BH mass function data --- colibre/auto_plotter/black_holes.yml | 16 ++++++++++------ 1 file changed, 10 insertions(+), 6 deletions(-) diff --git a/colibre/auto_plotter/black_holes.yml b/colibre/auto_plotter/black_holes.yml index 5bd6c2c9..f4827b50 100644 --- a/colibre/auto_plotter/black_holes.yml +++ b/colibre/auto_plotter/black_holes.yml @@ -470,9 +470,11 @@ black_hole_mass_function: - filename: BlackHoleMassFunction/Marconi2004_Data.hdf5 - filename: BlackHoleMassFunction/Shankar2004_Data.hdf5 - filename: BlackHoleMassFunction/MerloniHeinz2008_Data.hdf5 - - filename: BlackHoleMassFunction/Shankar2009_Data.hdf5 - - filename: BlackHoleMassFunction/Vika2009_Data.hdf5 - - filename: BlackHoleMassFunction/Caramete2010_Data.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_Mstar_observed.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_Mstar_unbiased.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_sigma_observed.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_sigma_unbiased.hdf5 + - filename: BlackHoleMassFunction/Gallo2019_Data.hdf5 black_hole_mass_function_adaptive: type: "adaptivemassfunction" @@ -496,6 +498,8 @@ black_hole_mass_function_adaptive: - filename: BlackHoleMassFunction/Marconi2004_Data.hdf5 - filename: BlackHoleMassFunction/Shankar2004_Data.hdf5 - filename: BlackHoleMassFunction/MerloniHeinz2008_Data.hdf5 - - filename: BlackHoleMassFunction/Shankar2009_Data.hdf5 - - filename: BlackHoleMassFunction/Vika2009_Data.hdf5 - - filename: BlackHoleMassFunction/Caramete2010_Data.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_Mstar_observed.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_Mstar_unbiased.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_sigma_observed.hdf5 + - filename: BlackHoleMassFunction/Shankar2016_sigma_unbiased.hdf5 + - filename: BlackHoleMassFunction/Gallo2019_Data.hdf5 From 4c34cee05ab8a015d9fca503533ce500fcf8337f Mon Sep 17 00:00:00 2001 From: Filip Husko Date: Thu, 23 May 2024 10:43:34 +0100 Subject: [PATCH 11/11] Update submodule --- observational_data | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/observational_data b/observational_data index e532543b..09f431d2 160000 --- a/observational_data +++ b/observational_data @@ -1 +1 @@ -Subproject commit e532543b1b24052b86eddb5c13f7c66634dc80fd +Subproject commit 09f431d2d88fadd0ee56aec3b5ce86b071631257