[mesa-users] trying to run RLO calculations for hot Jupiter systems

[Brian Jackson]

Rob,
Thanks for catching that mistake. I've fixed the planetary model
(attached). Unfortunately, that didn't seem to fix the problem I have,
though. Same error.

Brian

On Tue, Feb 28, 2017 at 4:47 PM, Robert Farmer <rjfarmer at asu.edu> wrote:
> Hi
> The nuclear networks (and hence the isotopes each model has defined) in the
> two model files are different, can you try in the planet inlist
>
> new_net_name = 'o18_and_ne22.net'
> change_initial_net = .true.
>
> so the planet model will have the same network as the star
>
> Rob
>
>
> On Tue, Feb 28, 2017 at 4:30 PM, Brian Jackson <bjackson at boisestate.edu>
> wrote:
>>
>> Dear MESA Users,
>> In my group, we're trying to continue our analysis of planetary
>> Roche-lobe overflow (RLO).
>>
>> For that, we generate starting planetary and stellar models
>> (attached). The star model is generated using the 1M_pre_ms_to_WD test
>> suite with only "max_age = 20e6" added to inlist_1.0 file. The
>> planetary model is generated using the make_planets test suite with
>> only "max_age = 20d6" changed in inlist_evolve_1.0_MJ_10.0_ME_2.0_RJ
>> file.
>>
>> We make some modifications to inlists for the binary_both_stars test
>> suite to load in the saved model files but otherwise use the standard
>> test suite.
>>
>> The following shows the last few lines from MESA before it stops
>> running our calculation:
>>
>> get_xa_for_accretion: accretion species mass fractions do not add to 1.0
>>                                      sum(xa(1:species))
>> 0.0000000000000000D+00
>>                                                      h1           1
>> 0.0000000000000000D+00
>>                                                     he3           2
>> 0.0000000000000000D+00
>>                                                     he4           3
>> 0.0000000000000000D+00
>>                                                     c12           4
>> 0.0000000000000000D+00
>>                                                     n14           5
>> 0.0000000000000000D+00
>>                                                     o16           6
>> 0.0000000000000000D+00
>>                                                    ne20           7
>> 0.0000000000000000D+00
>>                                                    mg24           8
>> 0.0000000000000000D+00
>>                    1                   1     2.000000000E+07
>> 9.244518984E-04    -9.900000000E+01    -3.245769587E-01
>>                    2                   1     2.000000000E+07
>> 9.999853565E-01    -9.900000000E+01    -6.799606683E-01
>>  have reached retry limit so now backup
>>  have reached retry limit so now backup
>>
>>                                                      dt
>> 6.3566365930376634D-07
>>                                      min_timestep_limit
>> 9.9999999999999995D-07
>>
>>  stopping because of convergence problems dt < min_timestep_limit
>>
>> terminated evolution: convergence problems
>>
>> termination code: min_timestep_limit
>>
>> ---
>>
>> The species mass fractions look like nonsense, even though they look
>> fine in the initial model files. And strangely, they change depending
>> on the initial orbital period we choose.
>>
>> I am running Mac OS X 10.9.5, MESA version 8845, and the MESA SDK
>> version from April 4 2016. I have attached the binary inlists and
>> initial models we use to run the calculation.
>>
>> Please let me know if you need additional information and if you have
>> any advice.
>>
>> Many thanks,
>> Brian Jackson
>>
>>
>> ------------------------------------------------------------------------------
>> Check out the vibrant tech community on one of the world's most
>> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: planet_init.mod
Type: application/octet-stream
Size: 213840 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/9baef859/attachment.obj>