[Mesa-users] Strange behavior for other_rsp_build_model

[Jeremy Sakstein]

Thanks Radec.

For my purposes, it turns out I just needed convergence to a single mode.

Cheers,

Jeremy

On Thu, Mar 19, 2020 at 4:42 PM Radek Smolec <smolec at camk.edu.pl> wrote:

> Hi Jeremy,
>
> Since the same model (the same physical parameters) seems to arrive at two
> different pulsation states, F and 1O, the initialization of the model is
> likely different for the two cases you describe - so called either-or
> scenario, see Section 2.2.4 in MESA V.
> To check whether it is indeed the case, you should enable more output, in
> particular look at the initial stage of model integration and compare the
> run of eg. radius vs. time. If it is different then that is it, and it
> remains to investigate why the inclusion of this particular line changes
> the initial kick.
> Cheers,
> Radek
>
> śr., 18 mar 2020 o 19:11 Jeremy Sakstein <sakstein at hawaii.edu> napisał(a):
>
>> Hi everyone,
>>
>> I have found some rather strange behavior when using my
>> own other_rsp_build_model.
>>
>> The routine below is a hack that Bill suggested a while ago. As far as I
>> am aware, it calls the usual routines that run when you call RSP except
>> that I was able to output what I need earlier and reduce the run-time.
>>
>> Now here is the strange thing: if I fix all other variables in the inlist
>> and don't include the line to use the other_rsp_build_model routine the
>> star pulsates in the fundamental mode. If I turn this on it appears to
>> pulsate in the first overtone.
>>
>> Am I just doing something silly here?
>>
>> Cheers,
>>
>> Jeremy
>>
>> subroutine js_other_rsp_build_model(id, ierr)
>>          integer, intent(in) :: id ! star id if available; 0 otherwise
>>
>>          integer, intent(out) :: ierr ! 0 means AOK.
>>          type (star_info), pointer :: s
>>          ierr = 0
>>          call star_ptr(id, s, ierr)
>>          if (ierr /= 0) return
>>          call star_do1_rsp_build(s,ierr)
>>
>>          write(*,*) 'linear periods:', s% rsp_lina_periods(1) / 86400, s%
>> rsp_l\
>> ina_periods(2)/86400, s% rsp_lina_periods(3)/86400
>>
>>          write(*,*) 'log(g):', s% photosphere_logg
>>          write(*,*) 'Radius:', 10.0**s% log_surface_radius
>>          !stop
>>          !write(*,*) 'no implementation for other_rsp_build_model'
>>
>>       end subroutine js_other_rsp_build_model
>> _______________________________________________
>> 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/20200325/29f450df/attachment.htm>