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

[Richard Battistini]

Hi Radek,

I just figured out the issue. Commenting out max_model_number was not the
solution, but setting it to -1 and rerunning was the solution.

Thanks,

On Thu, May 20, 2021 at 11:23 AM Radek Smolec <smolec at camk.edu.pl> wrote:

> Hi Richard,
>
> Are you using 15140? Have you copied your working directory from
> $MESA_DIR/star/work/ and did ./mk?
> Also, please send the inlists from your working directory.
> Radek
>
> czw., 20 maj 2021 o 17:11 Richard Battistini <richardbattistini2 at gmail.com>
> napisał(a):
>
>> 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/041c3d70/attachment.htm>