[Mesa-users] Relaxing Entropy from Profile Pairs?

[Charles Gibson]

Hello,

I am an undergraduate student working on taking a star evolved in MESA from
an SPH code (StarSmasher) and inputting it back into MESA to observe its
evolution. I have been working in the directory

star/test_suite/relax_composition_j_entropy

to ensure that I can accurately recreate a model from a previous mesa run.
So far, I have succeeded in replicating the star's evolution through this
method.

However, the StarSmasher simulations that I am using do not supply direct
values for entropy. I read in Appendix B of this paper
<https://iopscience.iop.org/article/10.3847/1538-4365/aaa5a8#apjsaaa5a8app2>
that
I can supply different pairs of profiles including (rho, T), (P_gas, T),
and (rho, e) in the entropy's place. These are all quantities that I can
easily obtain.

When I run the relax model using (rho, T) in place of entropy, the run
fails when processing the entropy file. I have it formatted the same as the
entropy file, but with three values on each line instead of two. Is there a
variable in the inlist that I have to alter to indicate that I am changing
the type of supplied entropy data? If it helps, I am using v10398.

Additionally, as a side question, when generating the number of shells and
each of their q values (how much mass is located outside the shell), is
there a specific formula that determines how I should determine those
values? I haven't gotten to the point where I am extracting values from
StarSmasher yet, but once I do, I want to make sure that I determine the
values so that MESA can effectively read them.

Thank you,
Charlie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230402/78fb264e/attachment.htm>