pyMesa
pyMesa copied to clipboard
Constants return incorrect value
Closely related to https://github.com/rjfarmer/pyMesa/issues/4 Running
import pyMesa as pym
const_lib, const_def = pym.loadMod("const")
print(const_def.Msun)
print(const_def.kerg)
prints
1.9884098706980504e+33
9.948624515087117
The first output is correct, but the second isn't. This may be related to gfort2py: the same issue shows up if I directly use gfort2py to read the MESA libraries.
I am using gfort2py 2.6.1, pyMesa 2.0.1, mesasdk 24.7.1 (in particular gfortran 15.2.1), and MESA 24.08.1 (although the bug persists with previous versions).
Thanks this is just the bug in #4 again, the code seems to have gotten lost at some point.
Can you try 9b7825f436a529f2168724c5c138763873f38c56 which bumps the gofrt2py version with the fix?