fxt
i second that notion. it'll cost in terms of resources though ...
it is available in r12778 as a run_star_extras % ls /Users/fxt/mesa_work/r12778/star/test_suite/other_physics_hooks/src % more /Users/fxt/mesa_work/r12778/star/test_suite/other_physics_hooks/src/convective_bdy/convective_bdy.inc perhaps a similar run_star_extras can be implemented in more recent releases.
and i remembered this exchange: ``` On May 26, 2020, at 10:28 AM, Bill Paxton wrote: please take a look at star/test_suite/other_physics_hooks where this is an implementation of convective_bdy_weight done...
whoa! freaking awesome!
running 246ac13b now. happy to dig into some of the test suite cases that may be suspect.
without defaults for all parameters, mesa will not work out of the box. users, especially less experiences users, will need to make decisions on things they likely do not know....
i agree requests for this kind of output keeps coming up. as you replied, a run_star_extras for the 2019 summer school was the last time i worked on it. maybe...
for a GIGO model, no negative or NAN input quantities - notionally T, Rho, composition. on the output side of GIGO, maybe add no negative entropies and specific heats.
another option on the input side is to only consider values within a valid range - no temperatures less than 10 K or larger than 1e13 K, no densities larger...
usually between steps i do things like x(i) = max(x(i), floor) where here floor is 1d-30 and optionally renormalize the positive definite mass fractions.