diff options
-rw-r--r-- | src/EOS_Hybrid_Startup.F | 28 |
1 files changed, 14 insertions, 14 deletions
diff --git a/src/EOS_Hybrid_Startup.F b/src/EOS_Hybrid_Startup.F index 24111d4..d298ff0 100644 --- a/src/EOS_Hybrid_Startup.F +++ b/src/EOS_Hybrid_Startup.F @@ -61,20 +61,20 @@ end if -c$$$ all these are now set by the polytropic EOS - -c$$$ These constants match the values in RNSID: consts.h - -c$$$ m_solar_cgs = 1.987d33 -c$$$ c_cgs = 2.9979d10 -c$$$ G_cgs = 6.6732d-8 -c$$$ -c$$$ m_solar_geom = G_cgs / c_cgs**2 * m_solar_cgs -c$$$ p_geom_factor = G_cgs / c_cgs**4 * m_solar_geom**2 -c$$$ rho_geom_factor = p_geom_factor * c_cgs**2 -c$$$ rho_geom_factor_inv = 1.d0 / rho_geom_factor -c$$$ -c$$$ eos_k_cgs = eos_k * rho_geom_factor**gamma_ini / p_geom_factor +cccc all these are now set by the polytropic EOS + +cccc These constants match the values in RNSID: consts.h + +cccc m_solar_cgs = 1.987d33 +cccc c_cgs = 2.9979d10 +cccc G_cgs = 6.6732d-8 +cccc +cccc m_solar_geom = G_cgs / c_cgs**2 * m_solar_cgs +cccc p_geom_factor = G_cgs / c_cgs**4 * m_solar_geom**2 +cccc rho_geom_factor = p_geom_factor * c_cgs**2 +cccc rho_geom_factor_inv = 1.d0 / rho_geom_factor +cccc +cccc eos_k_cgs = eos_k * rho_geom_factor**gamma_ini / p_geom_factor eos_k_supernuclear_cgs = eos_k_cgs * (rho_nuc * rho_geom_factor_inv)** . (eos_gamma - eos_gamma_supernuclear) |