[mesa-users] Zero Brunt-Vaisala freq. in outer regions of GYRE output files

[Bill Paxton]

On Jun 10, 2014, at 1:07 PM, Richard Townsend wrote:

> 
> On Jun 10, 2014, at 3:03 PM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
> 
>> 
>> On Jun 10, 2014, at 12:40 PM, Richard Townsend wrote:
>> 
>>> Is there a good reason why the Brunt can’t be calculated in the atmosphere? Since the acoustic cutoff frequency depends on the Brunt, this issue is quite important.
>> 
>> Confusion alert....    try setting this in your &controls inlist: 
>> 
>>       add_atmosphere_to_pulse_info = .true.
>> 
>> see star/private/create_atm.f and pulsation_info.f for details.
>> 
>> -b
>> 
>> 
>> 
> 
> Still confused — in store1_gyre_point (inside pulsation_info.f) the BV is explicitly set to zero. That doesn’t seem right.


ahhhh...... work to be done.  it seems that N2 is calculated in the atmosphere for FGONG, but not for GYRE.   See "store_atm_info" in pulsation_info.

please fix the routine that does gyre output and send me the updates.  

thanks,
b


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140610/f183ae21/attachment.html>