[Mesa-users] Load saved model vs create PMS model

[Pablo Marchant]

Hi Roque,

one thing I notice is that your inlist has

    new_rotation_flag = .true.
    change_rotation_flag = .true.
    set_initial_omega = .true.
    new_omega =  9.09025d-6

If you start with a pre-main sequence star you will set that value of omega
to a very bloated star. That might be the source of your problems.

Cheers

On Tue, Jun 9, 2020 at 8:41 PM Roque Caballero via Mesa-users <
mesa-users at lists.mesastar.org> wrote:

> Dear all,
>
> I'm trying to reproduce the MESA 2019 exercise (link
> <http://doi.org/10.5281/zenodo.3374959>) proposed by Professor Marc
> Pinsonneault and I'm getting some problems for which I'm not able to
> understand the reason.
>
> The inlist of that exercise used a a pre-computed PMS model. If a run the
> simulation with the unmodified inlist file, it works. On the other hand, if
> I configure the inlist for not using the pre-computed model but for
> generating it,
>
> *    create_pre_main_sequence_model = .true.*
> *    load_saved_model = .false.*
>
> I get the following error after a few seconds when I run the simulation:
>
>
> *  -2.0889505372223345E-012  -2.0395779779221372E-012
> -1.9901408914903830E-012  -1.9405185782435689E-012*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -5.3010299956639804D+00*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -5.6020599913279616D+00*
> *backup       1*
> *                           1st model retry log10(dt/yr)
> -7.6020599913279616D+00*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -7.9030899869919429D+00*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -8.2041199826559232D+00*
> *backup       1*
> *                           1st model retry log10(dt/yr)
> -1.0806179973983886D+01*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -1.1107209969647867D+01*
> * retry       1*
> *                           1st model retry log10(dt/yr)
> -1.1408239965311848D+01*
> *                                                     dt
> 7.7046264609375016D-08*
> *                                     min_timestep_limit
> 9.9999999999999995D-07*
>
> * stopping because of problems dt < min_timestep_limit*
>
>
> *terminated evolution: cannot find acceptable model*
>
> *termination code: min_timestep_limit*
> *DATE: 2020-06-09*
> *TIME: 20:31:14*
>
> Could any one help me to understand why is this happening? BTW, I'm using
> the MESA version 11701, the same used for this particular exercise.
>
> Many thanks in advance for your support.
>
> Best regards,
>
> Roque
>
> PS: The attached inlist file is the one producing the error
>
>
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
>

-- 
Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200609/1b5447b7/attachment.htm>