corrected output filename for d_l_gw

itrharrison · Feb 28, 2022 · 468d794 · 468d794
1 parent 548c948
commit 468d794
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/hi_class_interface.py b/hi_class_interface.py
@@ -244,7 +244,8 @@ def get_class_outputs(block, c_classy, config):
         # !!! why necessary to divide?
         alpha_mz = ([c_classy.alpha_m_at_z(zi)/c_classy.Omega_smg() for zi in z])
         block[distances, 'alpha_mz'] = alpha_mz
-        block[distances, 'd_l_gw'] = np.exp(0.5*cumtrapz(alpha_mz/(1. + z), z, initial=0)) # actually the ratio
+        block[distances, 'd_l_gw_on_d_l_em'] = np.exp(0.5*cumtrapz(alpha_mz/(1. + z), z, initial=0))
+        block[distances, 'd_l_gw'] = d_l * np.exp(0.5*cumtrapz(alpha_mz/(1. + z), z, initial=0))
 
     #Save some auxiliary related parameters
     block[distances, 'age'] = c_classy.age()