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

[Radek Smolec]

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/e18685a6/attachment.htm>