[Mesa-users] Questions about RSP

[Jeremy Sakstein]

Thanks Bill, hope the school is going well.

On Wed, Aug 14, 2019 at 11:03 AM Bill Paxton <paxton at kitp.ucsb.edu> wrote:

> On Aug 14, 2019, at 7:06 AM, Jeremy Sakstein <sakstein at sas.upenn.edu>
> wrote:
>
> Hi Bill,
>
> I just noticed that if I update MESA the rsp files have been changed. Is
> this enough for me to use cgrav(k) now or should I wait for an official
> release?
>
>
> Hi Jeremy,
>
> Please wait a bit longer.  There should be a new public release by the end
> of the month.   As you may know, we’ve delayed making the release until
> after the MESA Summer School and that’s going on right now.
>
> -Bill
>
>
>
> Cheers,
>
> Jeremy
>
> On Tue, Jul 16, 2019 at 5:00 PM Jeremy Sakstein <sakstein at sas.upenn.edu>
> wrote:
>
>> Thanks for the information Bill. I’ll look forward to the next release.
>>
>> On Tue, Jul 16, 2019 at 4:10 PM Bill Paxton <paxton at kitp.ucsb.edu> wrote:
>>
>>>
>>> On Jul 16, 2019, at 12:07 PM, Jeremy Sakstein <sakstein at sas.upenn.edu>
>>> wrote:
>>>
>>> Hi everyone,
>>>
>>> I'm playing around with the RSP test_suite cases and have two questions
>>> if that is okay:
>>>
>>> 1) In the pulsations paper it says that other hooks are still used. Does
>>> that include other_cgrav? In particular, is it used in e.g. equation (1)? I
>>> am trying to figure this out from the rsp files in star/private but it is
>>> not clear to me.
>>>
>>>
>>> rsp currently uses standard_cgrav, but in the next release it will start
>>> using cgrav(k) at point k, and that can be set by other_cgrav.
>>>
>>> 2) With regards to the new rsp quantities in the history file
>>> e.g. rsp_DeltaR and rsp_DeltaMag appear to be positive definite. At least
>>> when I plot them I only see positive values. Is this correct?
>>>
>>>
>>> you’ll find the following in rsp.f90.  LMAX and LMIN are the max and min
>>> luminosities during the most recent cycle.
>>>
>>>             s% rsp_DeltaMAG = 2.5d0*log10_cr(LMAX/LMIN)
>>> similarly for rsp_DeltaR where RMAX and RMIN are max and min during the
>>> most recent cycle.
>>>
>>>             s% rsp_DeltaR = RMAX-RMIN
>>>
>>>
>>> If so, I have been able to reproduce some figures in the paper showing
>>> the magnitude varying over one cycle  (e.g. fig 6) by simply picking a
>>> cycle and subtracting [image: -2.5\log(L_i)] where [image: L_i] is the
>>> input luminosity from my inlist. Is this the only way to do it or is there
>>> an internal variable I can plot.
>>>
>>> 3) Which quantity in the history_columns corresponds to the radial
>>> velocity plotted in the paper?
>>>
>>>
>>> rsp uses the variable v for radial velocities just as is done in the
>>> standard mesa/star with v_flag true.
>>>
>>> bill
>>>
>>>
>>>
>>>
>>>
>>> Cheers,
>>>
>>> Jeremy
>>> _______________________________________________
>>> 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/20190814/687e697a/attachment.html>