Evolutionary Models

Setting the model_type parameter of the &model namelist group to 'EVOL' tells the frontend to read the equilibrium stellar model from a file created by a stellar evolution code (e.g., MESA).

Supported Formats

The format of the model file is specified by the file_format parameter of the &model namelist group (see the Stellar Model Parameters section). Possible choices are summarized in the table below.

`file_format`	Description
`'AMDL'`	Binary file describing an evolutionary model in AMDL format, as reverse engineered from the ADIPLS stellar oscillation code (Christensen-Dalsgaard, 2008)
`'B3'`	HDF5 file describing an evolutionary model in B3 format. This format is for testing purposes only, and will eventually be superseded and/or removed
`'FAMDL'`	Text file describing an evolutionary model in FAMDL format, as specified in the `CoRoT/ESTA File Formats` document
`'FGONG'`	Text file describing an evolutionary model in FGONG format, as specified in the updated `FGONG Format` document
`'GSM'`	HDF5 file describing an evolutionary model in GYRE Stellar Model (GSM) format, as specified in the GSM File Format section
`'MESA'`	Text file describing an evolutionary model in MESA format, as specified in the MESA File Format section
`'LOSC'`	Text file describing an evolutionary model in the revised LOSC format
`'OSC'`	Text file describing an evolutionary model in OSC format, as specified in the `CoRoT/ESTA File Formats` document)
`'WDEC'`	Text file describing an evolutionary model in WDEC format, as specified in Bischoff-Kim & Montgomery (2018)

Interpolation

Cubic spline interpolation is used to evaluate data between model grid points. The deriv_type parameter in the &model namelist group controls how the spline derivatives are set up.

Double Points

If a model contains a pair of adjacent points with the same radial coordinate \(r\), this pair is treated as a double point representing a discontinuity in the density and some other thermodynamic quantities (but not the pressure or temperature). GYRE does not attempt to interpolate across double points, but instead handles them properly when solving equations through the use of internal boundary conditions.