diff --git a/src/MINT/MINT_data_columns_MS_list.def b/src/MINT/MINT_data_columns_MS_list.def index c2b161bc3ba91b3fa9be8be7af920c1117d33f49..d5147ba0beadf8350a75f7958e37003c37cae373 100644 --- a/src/MINT/MINT_data_columns_MS_list.def +++ b/src/MINT/MINT_data_columns_MS_list.def @@ -75,6 +75,7 @@ X(K2, NOTHING) \ X(TIDAL_E2, NOTHING) \ X(TIDAL_E_FOR_LAMBDA, NOTHING) \ + X(MOMENT_OF_INERTIA_FACTOR, NOTHING) \ X(TIMESCALE_KELVIN_HELMHOLTZ, NOTHING) \ X(TIMESCALE_DYNAMICAL, NOTHING) \ X(TIMESCALE_NUCLEAR, NOTHING) \ diff --git a/src/MINT/MINT_data_tables.def b/src/MINT/MINT_data_tables.def index db1a7f1e0b49d4e72e6778a51eea712f9da626c5..a8cbff7bf636c689bf1befc5f4ba0c7842f1c74f 100644 --- a/src/MINT/MINT_data_tables.def +++ b/src/MINT/MINT_data_tables.def @@ -25,7 +25,7 @@ X(MS,"MS/grid_ZZ_chebyshev_withE.dat",&MINT_load_MS_grid,TRUE) \ X(HG,NULL,NULL,FALSE) \ X(GB,NULL,NULL,FALSE) \ - X(CHeB,"CHeB/grid_ZZ_chebyshev_withE.dat",&MINT_load_CHeB_grid,FALSE) \ + X(CHeB,NULL, NULL, FALSE) \ X(EAGB,NULL,NULL,FALSE) \ X(TPAGB,NULL,NULL,FALSE) \ X(HeMS,NULL,NULL,FALSE) \ diff --git a/src/MINT/MINT_generic_data_types.def b/src/MINT/MINT_generic_data_types.def index ac1cd6ace4ee3b80598d4c434f3b2846d8479334..7cc547b48fc6c6924fd5d9b750e25e02cc576ab2 100644 --- a/src/MINT/MINT_generic_data_types.def +++ b/src/MINT/MINT_generic_data_types.def @@ -27,6 +27,7 @@ X(K2) \ X(TIDAL_E2) \ X(TIDAL_E_FOR_LAMBDA) \ + X(MOMENT_OF_INERTIA_FACTOR) \ X(TIMESCALE_KELVIN_HELMHOLTZ) \ X(TIMESCALE_DYNAMICAL) \ X(TIMESCALE_NUCLEAR) \ diff --git a/src/MINT/MINT_stellar_structure_CHeB.c b/src/MINT/MINT_stellar_structure_CHeB.c index f28b30a22959f38b7dea46d021d5aac667353981..859b7596a9fafb60b7ccfe8d3353e9ecc01430d6 100644 --- a/src/MINT/MINT_stellar_structure_CHeB.c +++ b/src/MINT/MINT_stellar_structure_CHeB.c @@ -282,8 +282,13 @@ Stellar_type MINT_stellar_structure_CHeB(struct stardata_t * RESTRICT const star newstar->age, newstar->epoch); - newstar->k2 = get_data(MINT_CHeB_K2); + /* + * The line above should be replaced with + * newstar->k2 = get_data(MINT_CHeB_MOMENT_OF_INERTIA_FACTOR); + * once that quantity is available in CHeB grids + * The K2 above is MESA's and not I/(MR^2) + */ newstar->E2 = available[MINT_CHeB_TIDAL_E2] ? get_data(MINT_CHeB_TIDAL_E2) : diff --git a/src/MINT/MINT_stellar_structure_MS.c b/src/MINT/MINT_stellar_structure_MS.c index 3378a094b49d15935bd7d51b5b089e24d4cb1a24..363a033959980c330ce3c66ff8bd1fd81badd59e 100644 --- a/src/MINT/MINT_stellar_structure_MS.c +++ b/src/MINT/MINT_stellar_structure_MS.c @@ -293,7 +293,11 @@ Stellar_type MINT_stellar_structure_MS(struct stardata_t * RESTRICT const starda newstar->starnum, stardata->model.dt); - newstar->k2 = get_data(MINT_MS_K2); + newstar->k2 = get_data(MINT_MS_MOMENT_OF_INERTIA_FACTOR); + if(vb)printf("Moment of inertia factor k2=I/(MR^2) : star = %d, k22 = %g\n", + newstar->starnum, + stardata->k2); + newstar->E2 = available[MINT_MS_TIDAL_E2] ? get_data(MINT_MS_TIDAL_E2) :