Omega glitch after remeshing
- Dominant language
- Fortran
- Stars
- 249
- Forks
- 82
- Avg merge
- 1d 3h
- Merged PRs (30d)
- 3
Description
I've been doing some tests with red giant models including rotation (right now using the default settings, so only accounting for conservation of angular momentum) and I'm seeing some glitches that I think are introduced by remeshing.
To illustrate, here are some Omega profiles in the core region of a somewhat evolved red giant (log g ~3):
As you can see, this region gets remeshed after model number 1545 and then Omega exhibits a glitch for every model after this point. If I run these same models with `okay_to_remesh = .false.` then this glitch never develops.
This, (as well as the fact that in the remeshed models the inner zone seems to overestimate Omega whereas the outer zone seems to underestimate it) suggests to me that there's some issue with how the angular momentum is getting distributed when the model is remeshed.
Here's a work dir with a starting model right before this glitch develops for ease of replication:
[rotation_mesh_test.zip](https://github.com/user-attachments/files/31415651/rotation_mesh_test.zip)
Version information:
- MESA version: 26.04.1
- MESA SDK VERSION: mesasdk-x86_64-linux-26.3.2
- OS: Ubuntu 24.04.4 LTS
Any insight would be much appreciated!
Contributor guide
Assessment
This issue has not been assessed yet.