[mesa-users] negative age and rotation values in "OSC" file
Juan Carlos Suárez
jcsuarez at iaa.es
Fri Feb 10 13:51:14 EST 2017
Hi!
The last version
JC
> On 10 Feb 2017, at 19:29, RICHARD H D TOWNSEND <townsend at astro.wisc.edu> wrote:
>
> Hi Juan-Carlos —
>
> Very interesting! Can you remind me which version of MESA you’re running with?
>
> cheers,
>
> Rich
>
>> On Feb 10, 2017, at 12:26 PM, Juan Carlos Suárez <jcsuarez at iaa.es> wrote:
>>
>> Dear Richard,
>>
>> We have modified a little bit the previous “inlist” in order to include some atmosphere (the default simple atmosphere)
>> since we are interested in working with the photosphere radius. So we decided to include the following options in the
>> inlist (attached)
>> ! if true, write atmosphere to pulse files
>>
>> add_atmosphere_to_pulse_info = .TRUE.
>> ! if true, add point for r=0 to pulse files
>>
>> add_center_point_to_pulse_info = .FALSE.
>> ! if true, add k=1 cell to pulse files
>>
>> keep_surface_point_for_pulse_info = .false.
>> ! add double points at discontinuities
>>
>> add_double_points_to_pulse_info = .TRUE.
>>
>> Curiously, it seems that now the problem of the negative ages and rotation values is solved, but we
>> do not know why.
>>
>> On the other hand unfortunately the computation halted with an error
>>
>> failed to open LOGS/profile686.data
>> star_finish_step ierr 28
>> after_step_loop ierr 28
>>
>> That is, it halted at model 686, while without the above options MESA reached the 905 models with no errors.
>>
>> In addition to this, with the above options it seems that MESA has doubled the number
>> of mesh points (from 4K to 8K). Is that what we should expect? The 4K additional
>> points are located at the atmosphere?
>>
>> Thanks!
>> JC
>>
>>
>>> On 9 Feb 2017, at 20:38, Juan Carlos Suárez <jcsuarez at iaa.es> wrote:
>>>
>>> Yes here it is!
>>> JC
>>>> On 9 Feb 2017, at 20:05, RICHARD H D TOWNSEND <townsend at astro.wisc.edu> wrote:
>>>>
>>>> Hi Juan —
>>>>
>>>> Can you post the inlist file(s) that you’re using?
>>>>
>>>> Many thanks,
>>>>
>>>> Rich
>>>>
>>>>> On Feb 9, 2017, at 12:43 PM, Juan Carlos Suárez <jcsuarez at iaa.es> wrote:
>>>>>
>>>>> Dear colleagues
>>>>>
>>>>> I am trying to compute a 1.5 Msun evolutionary model (with solar metallicity) from PMS to TAMS approx.
>>>>> We have set the options to obtain the .OSC files to be read by our oscillation codes.
>>>>>
>>>>> We found that for all the models the age provided has negative values.
>>>>>
>>>>> Moreover, although no option for rotation was included in the inlist (i.e. rotation should be zero)
>>>>> we obtain some models with rotation 0 but others with negative non-zero values
>>>>>
>>>>> Any idea?
>>>>>
>>>>> Thanks in advance!
>>>>>
>>
>> --------------------------------------------------------------------------------------------------------------------
>> Dr. Juan Carlos Suárez
>> Ramon y Cajal fellow at Física Teórica y del Cosmos Dept. University of Granada
>> Facultad de Ciencias. 18071. Granada. Spain
>> jcsuarez at ugr.es | +34 958249056 | www.ugr.es/~fqm292/
>> Associate researcher at Stellar Physics Dept. IAA-CSIC
>> Glorieta de la Astronomia s/n. 18008. Granada. Spain
>> jcsuarez at iaa.es | +34 958230619 | www.iaa.es
>> --------------------------------------------------------------------------------------------------------------------
>>
>>
>>
>> <inlist_project>
>
--------------------------------------------------------------------------------------------------------------------
Dr. Juan Carlos Suárez
Ramon y Cajal fellow at Física Teórica y del Cosmos Dept. University of Granada
Facultad de Ciencias. 18071. Granada. Spain
jcsuarez at ugr.es | +34 958249056 | www.ugr.es/~fqm292/
Associate researcher at Stellar Physics Dept. IAA-CSIC
Glorieta de la Astronomia s/n. 18008. Granada. Spain
jcsuarez at iaa.es | +34 958230619 | www.iaa.es
--------------------------------------------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170210/89f95d50/attachment.html>
More information about the Mesa-users
mailing list