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

[Radek Smolec]

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/20210519/79f6c389/attachment.htm>