[Mesa-users] Equivalents to use_zero_dLdm_outer_BC in newer MESA releases

[Benjamin Herfray]

Hello MESA users,

I've been trying to run simulations of stellar winds during type-I X-ray bursts using version r24.03.1. In order to try to counter problems with density inversions developing near the outer layers of the atmosphere I would like to impose boundary conditions on the luminosity—specifically, I would like to impose dL/dm = 0 at the outer boundary.​ This is similar to the methods described in Yu & Weinberg, 2018 (https://doi.org/10.3847/1538-4357/aad045) and Guichandut & Cumming, 2021 (https://doi.org/10.3847/1538-4357/ace43c) who used older versions of MESA.

Older versions of MESA (I believe up to r15140) have the control use_zero_dLdm_outer_BC to impose L(1) = L(2), however this was removed in subsequent versions with no explanation of the change in the changelog.

Does anyone know why this was removed, and is there a combination of controls present in the current version of MESA that I can use to have the equivalent effect as use_zero_dLdm_outer_BC? Or would I have to implement this as a custom condition through run_star_extras?

Thank you,
Benjamin Herfray
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240709/b45814c5/attachment.htm>