MESAHub / MESAHub/mesa

Bumps when loading custom rates

Open
#421 11 comments 0 reactions 0 assignees View on GitHub
bug rates-net
Dominant language
Fortran
Stars
249
Forks
82
Avg merge
1d 3h
Merged PRs (30d)
3

Description

Interpolating over the rate table seems to be adding bumps in the rate (this is the issue reported by Francisco De Geronimo)

For the attached triple alpha rate see the plot for what happens at low T. Sure the rate has only 23 points so is under-sampled but I would have expected a smoother curve from MESA.

![cache_rate](https://user-images.githubusercontent.com/10981326/172832199-b44cd2e2-6ee0-4259-901d-a1e7bd35aaa5.png)

Input rate
[r_he4_he4_he4_to_c12.txt](https://github.com/MESAHub/mesa/files/8869676/r_he4_he4_he4_to_c12.txt)

MESA's cached rate

[r_he4_he4_he4_to_c12_0.txt](https://github.com/MESAHub/mesa/files/8869678/r_he4_he4_he4_to_c12_0.txt)

Maybe its an issue with having only a few points near the boundary? As at higher T's the curve is smooth.

Contributor guide

Open the contributing guide

Research direction

Start by reproducing the low-temperature bumps using the attached input rate and MESA's cached rate, then inspect the Fortran code responsible for loading and interpolating custom rates. Done means identifying and correcting the interpolation behavior so the cached curve is smooth at low temperatures while preserving the higher-temperature behavior.

Written by the indexing model from the issue text.

Assessment

Tech stack
fortran
Domain
data
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.