Conversation
…ofile columns flux_limit_lambda and flux_limit_R (ported from r9793)
Porting the flux-limiter prescription from r9793 to the latest MESA version
|
Let's make sure there are no floating point errors associated with this. I don't know if we need to save any star information in the photos, but let's see what happens to : all use dPrad_dm form of the T-gradient eqn. https://testhub.mesastar.org/feature%2Freintroduce_flux_limiter/commits/9f8a91c. I'll add a fpe check as well. |
star/defaults/controls.defaults
Outdated
| ! :: | ||
|
|
||
| use_dPrad_dm_form_of_T_gradient_eqn = .false. | ||
| use_flux_limiting_with_dPrad_dm_form = .true. |
There was a problem hiding this comment.
reminder: switch back to false before merge
|
This pr needs testing and a changelog before it can be merged. |
| !lum_rad_div_L_Edd_sub_fourPrad_div_PchiT ! density increases outward if this is > 0 | ||
| ! see Joss, Salpeter, and Ostriker, "Critical Luminosity", ApJ 181:429-438, 1973. | ||
| Frad_div_cUrad | ||
|
|
There was a problem hiding this comment.
need to add :
flux_limit_lambda
flux_limit_R
|
I put The |
4 participants
see #771. where @JRGarza has reimplemented a flux limiter for use with the dPrad_dm form of the T-gradient equation.
Resolves #769
Related to #625