[Mesa-users] Running super-solar metallicity EGPs

[Aaron Dotter]

Hi Christopher,

I've tried running your inlist and you're running into a hard-coded limit
in the create_initial_model routine that stops when T < 150 K at some point
in the model.  In this case it's fairly easy to find where the problem
occurs since it never gets out of the create_intial_model stage of the
process.  However, from there it may be difficult to figure out exactly
what's going wrong.  In general, when tracking down an error, grepping for
the error message in the source code is a decent way to get started.

You might start with the original inlists and initial_Z=0.02 value in
make_planets and gradually increase from (say 0.03, 0.04, etc.) just to see
where you start running into problems.  Although the Freedman low-T
opacities are more than up to the task of going to Z=0.1, you may find that
other parts of the code are not.  Other than that, I suggest you experiment
with the radius of the initial model for a given mass to see if that helps.

Aaron



On Mon, Jan 7, 2019 at 12:28 PM C.R. Seay <seaychristopher.r at gmail.com>
wrote:

> Whoops! Sure thing. Sorry for the vague initial email!
>
> The relevant inlist is inlist_create, MESA version is r10398.
>
> Thank you,
> --Christopher
>
>
> On Mon, Jan 7, 2019 at 4:38 PM Aaron Dotter <aaron.dotter at gmail.com>
> wrote:
>
>> Hi,
>>
>> We'll need more information to help you.  Please send your inlists and
>> tell us what MESA version you're using.
>>
>> Thanks,
>> Aaron
>>
>> On Mon, Jan 7, 2019 at 10:36 AM C.R. Seay via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>>
>>> Hello all,
>>>
>>> I'm trying to run a hot Jupiter with 3x, 5x metallicity but my models
>>> are failing to converge. I'm also not sure how to debug here. (Actually, as
>>> an aside, if anyone can give some tips with debugging, in general, that
>>> would be appreciated!)
>>>
>>> For example, my 3x solar model doesn't converge because the "temp is too
>>> low in integration." I've read the 2013 paper, but I'm not quite sure what
>>> that's really referring to, or what parameters I need to mess around with
>>> to change it. I'm assuming it's referring to the adiabatic temperature
>>> profile that MESA uses to build the model.
>>>
>>> Thanks for any help!
>>> --Christopher
>>>
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>
>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190107/87a7f86b/attachment.html>