[Mesa-users] Questions about RSP

[Jeremy Sakstein]

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?

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/d0661f73/attachment-0001.html>