Added references to WD constants

Author: SimonStevenson

src/constants.h

@@ -430,19 +430,19 @@ constexpr double Q_HYDROGEN_BURNING                     = 6.4E18 * MSOL_TO_G / (
 constexpr double L0_RITTER_HIGH_Z                       = 1995262.3;                                                // Luminosity constant which depends on metallicity in Ritter 1999, eq 10
 constexpr double L0_RITTER_LOW_Z                        = 31622.8;
 constexpr double WD_HE_SHELL_MCRIT_DETONATION           = 0.05;                                                     // Minimum shell mass of He for detonation. Should be composed of helium (so, exclude burnt material), but not implemented yet. Ruiter+ 2014.
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_RG_SS_0      = -6.84;
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_RG_SS_1      = 1.349;
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SS_MF_0      = -8.115;
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SS_MF_1      = 2.29;
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SF_Dt_0      = -8.313;
-constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SF_Dt_1      = 1.018;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_0      = -8.3302;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_1      = 2.8825;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_2      = -0.9802;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_0        = -9.2176;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_1        = 3.5732;
-constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_2        = -1.2138;
-constexpr double WD_MP                                  = 5.7E-4;
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_RG_SS_0      = -6.84;                                                    // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_RG_SS_1      = 1.349;                                                    // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SS_MF_0      = -8.115;                                                   // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SS_MF_1      = 2.29;                                                     // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SF_Dt_0      = -8.313;                                                   // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_PIERSANTI_SF_Dt_1      = 1.018;                                                    // Constant from Piersanti et al. 2014 see Table A1
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_0      = -8.3302;                                                  // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_1      = 2.8825;                                                   // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_REDGIANT_2      = -0.9802;                                                  // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_0        = -9.2176;                                                  // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_1        = 3.5732;                                                   // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_LOG_MT_LIMIT_NOMOTO_STABLE_2        = -1.2138;                                                  // Quadratic fit to results from Nomoto et al. 2007 Table 5
+constexpr double WD_MP                                  = 5.7E-4;                                                   // White dwarf mass parameter determining where mass-radius relation changes for low-mass white dwarfs
 constexpr double WD_BELCZYNSKI_SN_CONSTANT              = 1.34;                                                     // Constant from Eq 62, Belczynski+ 2008
 constexpr double WD_BELCZYNSKI_SN_LINEAR                = 4.0E8;                                                    // Linear term factor from Eq 62, Belczynski+ 2008
 constexpr double WD_BELCZYNSKI_IMMEDIATE_FLASH          = 1.64E-6;                                                  // Accretion limit from eq 61, Belczynski+ 2008.