[Mesa-users] "good match for period" stopping RSP_cepheid nonlinear calculations

[Richard Battistini]

Hello again,

I followed the suggestions Radek, but even after doing so, my nonlinear
calculations stop after 10 periods with the same terminal output.

I thought it was because the controls were still left on in the file
inlist_rsp_common, but this file does not even exist in this version of
MESA.

Is there something I am missing here?

On Wed, May 19, 2021 at 3:24 PM Radek Smolec <smolec at camk.edu.pl> wrote:

> Hi Richard,
>
> As you noticed, the example in test_suite integrates only for 10 periods,
> and checks the period with value passed via inlist:
>    x_integer_ctrl(1) = 10 ! which period to check
>    x_ctrl(1) = 3.92305d0 ! expected period (in days)
> using run_star_extras. If you want to compute many cycles and possibly for
> different model parameters, I suggest starting with a copy of standard mesa
> working directory. Then, copy paste the inlists from test_suite, remove the
> above two lines as well as max_model_number. To integrate for a given
> number of cycles specify RSP_max_num_periods in &controls section of the
> inlist.
> Best wishes,
> Radek
>
> śr., 19 maj 2021 o 19:36 Richard Battistini via Mesa-users <
> mesa-users at lists.mesastar.org> napisał(a):
>
>> Hi Rob,
>>
>> After doing the search, I come across these lines in
>> src/run_star_extras.f90
>>
>>          if (s% total_energy /= 0d0 .and. abs(rel_run_E_err) > 1d-5) then
>>             write(*,*) '*** BAD rel_run_E_error ***', &
>>             s% cumulative_energy_error/s% total_energy
>>          else if (abs(s% rsp_period/(24*3600) - target_period) > 1d-2)
>> then
>>             write(*,*) '*** BAD ***', s% rsp_period/(24*3600) -
>> target_period, &
>>                s% rsp_period/(24*3600), target_period
>>          else
>>             write(*,*) 'good match for period', &
>>                s% rsp_period/(24*3600), target_period
>>          end if
>>
>> Which, if I understand correctly, means that if the rsp_period is close
>> enough to the target period, the string 'good match for period' is printed
>> and the calculations terminate.
>>
>>
>> On Wed, May 19, 2021 at 1:17 PM Rob Farmer <robert.j.farmer37 at gmail.com>
>> wrote:
>>
>>> Hi,
>>>
>>> Have you tried searching all files in the rsp_cephied folder for the
>>> string 'good match for period' ? What do you learn if you do?
>>>
>>> Rob
>>>
>>> On Wed, 19 May 2021 at 17:53, Richard Battistini via Mesa-users <
>>> mesa-users at lists.mesastar.org> wrote:
>>>
>>>> Hello,
>>>>
>>>> I am running mesa-r15140 on Ubuntu 64 bit, Virtual Machine. I am
>>>> trying to run a maximum of 4000 periods of nonlinear calculations for the
>>>> RSP_Cepheid test case. Yet, even when I believe all the relevant stopping
>>>> conditions, (max_model_number, x_integer_ctrl(1), x_ctrl(1)) is commented
>>>> out, the model does not run to 4000 periods. It stops around 10 periods
>>>> with the message,
>>>>
>>>>  rel_run_E_err   1.8311209038468803E-011
>>>>  good match for period   3.9239601544671396        3.9230499999999999
>>>>
>>>>
>>>> I have also attached the complete terminal output and the inlist for
>>>> your convenience. I was wondering if someone could help understand what I
>>>> am missing here.
>>>>
>>>> Thanks in advance,
>>>> Richard
>>>> _______________________________________________
>>>> mesa-users at lists.mesastar.org
>>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>>
>>>> _______________________________________________
>> 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/20210520/575e21a2/attachment.htm>