[Mesa-users] How to get Prad and Pgas from RSP models with message: P <= Prad

[Radek Smolec]

Hello,

You have pushed RSP beyond the limits it was tested and designed for. It
fails for models with very high temperatures ~10 000K, so we are well
beyond the blue edge of the instability strip (you may see that for a bit
cooler models the growth rates are all negative).
RR Lyrae pulsations are not expected there and RSP model builder fails in
this regime.
Radek

wt., 29 gru 2020 o 20:18 Richard Battistini <richardbattistini2 at gmail.com>
napisał(a):

> Hello there, for the original inlist, there are not any significant
> modifications besides the number of modes.
>
> As for the physical parameters, I have attached a file (input.dat) of the
> list of the parameters that are read in, ordered label, Z, X, mass,
> luminosity, effective temperature.
> I have also included the terminal output after running the grid.
>
>
> On Tue, Dec 29, 2020 at 2:01 PM Radek Smolec <smolec at camk.edu.pl> wrote:
>
>> Hello Richard,
>>
>> With 11701 the original RR Lyrae test_suite case runs smoothly on my
>> computer.
>> What are your modifications to the inlists? It seems that for some
>> physical parameters you are using RSP is not able to construct the envelope
>> and hence there is no output.
>> Cheers,
>> Radek
>>
>> wt., 29 gru 2020 o 19:20 Richard Battistini via Mesa-users <
>> mesa-users at lists.mesastar.org> napisał(a):
>>
>>> Hello all,
>>>
>>> I am running mesa-r11701. Whenever I run some models in the
>>> "rsp_RR_Lyrae" test case, I have the message:
>>>
>>> ```
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> *gfortran -fopenmp -o ../star run_star_extras.o run_star.o  run.o
>>>  -L../../../../lib -lstar -lgyre -lionization -latm -lcolors -lnet -leos
>>> -lkap -lrates -lneu -lchem -linterp_2d -linterp_1d -lnum -lmtx -lconst
>>> -lutils -lf2crlibm -lcrlibm  `mesasdk_lapack_link` `mesasdk_blas_link`
>>> `mesasdk_hdf5_link`  `mesasdk_pgplot_link` -lz run
>>> inlist_rsp_RR_Lyrae_headerDATE: 2020-12-22TIME: 18:32:22
>>>                      version_number       11701 read
>>> ../../inlist_rsp_common read inlist_rsp_RR_Lyrae create initial RSP
>>> modelSTOP build_rsp_model P <= Prad          -1   14.338359330232963
>>>  18.034728922739919        9195.9049210975136      failed in
>>> do_rsp_buildDATE: 2020-12-22TIME: 18:32:23finished
>>> inlist_rsp_RR_Lyrae_headerDATE: 2020-12-22TIME: 18:32:23finished inlists*
>>> ```
>>>
>>> Whenever this message is printed, there are no LOG files or any LINA
>>> calculations outputted for the model. Is there a way to still grab
>>> information about the radiation pressure, gas pressure, and zone for these
>>> "failed" models so that I may plot the radiation pressure against zone and
>>> the gas pressure against zone?
>>>
>>> Thanks,
>>> *Richard Battistini*
>>> *Undergraduate in Physics *
>>> _______________________________________________
>>> 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/20201229/7c5e8957/attachment.htm>