[Mesa-users] Questions about RSP
paxton at kitp.ucsb.edu
Tue Jul 16 16:09:23 EDT 2019
> 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 where 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.
> mesa-users at lists.mesastar.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users