Evan Bauer
Evan Bauer
`erf` (e.g. https://gcc.gnu.org/onlinedocs/gfortran/ERF.html) does not produce bit-for-bit identical results across compilers, so is not currently appropriate for general usage in MESA in the occasional circumstance where one might want to...
Anne suggests that the documentation for each test case should give users an indication at the start of the documentation of whether it falls into one or both of the...
Right now, our default is to do linear interpolation in composition (X,Z) for opacity tables. Even though cubic interpolation returns more realistic opacities, it returns much worse derivatives. Turning on...
Add in modifications due to magnetic fields as calculated in https://ui.adsabs.harvard.edu/abs/2019ApJ...870L...5H/abstract Adrian Fraser has a python implementation for this that we should be able to port into fortran for MESA....
https://github.com/MESAHub/mesa/blob/1baa4348b5bb1c12a487e9caab6264daef919af4/turb/private/thermohaline.f90#L116-L119
When the composition diffusivity enters into thermohaline calculations (e.g. when `thermohaline_option = Brown_Garaud_Stellmach_13`), currently we assume a H/He-dominated mixture. See https://github.com/MESAHub/mesa/blob/d215b38166687af0351c621f47ddf622415625d7/turb/private/thermohaline.f90#L128-L153 We already have all of the physics in the...
On the distillation branch, Rob pointed out to me that there's a risk of stack overflow when using explicit shape automatic arrays, for example to declare an extra array for...
See discussion on c58b3fc7cae532008ef437231033dd86718764b2, which was reverted. The ideal ion free energy in Skye gives negative entropy at sufficiently high density. This is a problem not just for Skye, but...
I happened to have these settings while working today: ``` profile_interval = 1 history_interval = 10 ``` Instead of writing a profile every step, MESA only wrote one every 10...