[Mesa-users] Questions about RSP

[Bill Paxton]

> 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.

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