[Mesa-users] Questions about RSP

[Bill Paxton]

> 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 <mailto: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 <mailto:paxton at kitp.ucsb.edu>> wrote:
> 
>> On Jul 16, 2019, at 12:07 PM, Jeremy Sakstein <sakstein at sas.upenn.edu <mailto: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 <mailto:mesa-users at lists.mesastar.org>
>> https://lists.mesastar.org/mailman/listinfo/mesa-users <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/3c31492f/attachment.html>