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

[Richard Battistini]

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
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210519/98af92ae/attachment.htm>