Oscillation Parameters

The &osc namelist group defines oscillation parameters; the input file can contain one or more, but only the last (matching) one is used. Allowable parameters are:

inner_bound (default 'REGULAR')

Inner boundary conditions; one of:

  • 'REGULAR' : Regularity-enforcing (only valid when inner grid point is at \(x = 0\))
  • 'ZERO_R' : Zero radial displacement (only valid when inner grid point is at \(x \neq 0\))
  • 'ZERO_H' : Zero horizontal displacement (only valid when inner grid point is at \(x \neq 0\))
outer_bound (default 'VACUUM')

Outer boundary conditions; one of:

  • 'VACUUM' : Zero surface pressure
  • 'DZIEM' : Formulation following Dziembowski (1971)
  • 'UNNO' : Formulation following Unno et al. (1989)
  • 'JCD' : Formulation following Jörgen Christensen-Dalsgaard (ADIPLS)
variables_set (default 'GYRE')

Dependent variables in oscillation equations; one of:

  • 'GYRE' : GYRE formulation, as desciribed in the equations.pdf document
  • 'DZIEM' : Formulation following Dziembowski (1971)
  • 'JCD' : Formulation following Jörgen Christensen-Dalsgaard (ADIPLS)
  • 'MIX' : Mixed formulation ('JCD' for gravitational components, 'DZIEM' for mechanical components)
  • 'LAGP' : Lagrangian pressure perturbation formulation
alpha_gr (default 1.)
Scaling factor for gravitational potential perturbations (see \(\alpha_{\rm gr}\) term in the osc-equation section). Set to 0. to implement the Cowling approximation
alpha_th (defaualt 1.)
Scaling factor for the thermal timescale (see \(\alpha_{\rm th}\) term in the osc-equations section). Set to 0. to implement the non-adiabatic reversible (NAR) approximation, and to large values to approach the adiabatic limit
alpha_hf (defaualt 1.)
Scaling factor for horizontal flux perturbations (see \(\alpha_{\rm hf}\) term in the osc-equations section). Set to 0. to implement the non-adiabatic radial flux (NARF) approximation
inertia_norm (default 'BOTH')

Inertia normalization factor; one of

  • 'RADIAL' : Radial amplitude squared, \(|\xi_{\rm r}|^{2}\), evaluated at x_ref
  • 'HORIZ' : Horizontal amplitude squared, \(|\lambda| |\xi_{\rm h}|^{2}\), evaluated at x_ref
  • 'BOTH' : Overall amplitude squared, \(|\xi_{\rm r}|^{2} + |\lambda| |\xi_{\rm h}|^{2}\), evaluated at x_ref
time_factor (default 'OSC')

Time-dependence factor in pulsation equations; one of:

  • OSC : Oscillatory, \(\propto \exp(-{\rm i} \omega t)\)
  • EXP : Exponential, \(\propto \exp(-\omega t)\)
conv_scheme (default 'FROZEN_PESNELL_1'`)

Scheme for treating convection; one of:

  • 'FROZEN_PESNELL_1' : Freeze convective heating altogether; case 1 described by Pesnell (1990)
  • 'FROZEN_PESNELL_4' : Freeze Lagrangian perturbation of convective luminosity; case 4 described by Pesnell (1990)
zeta_scheme (default PESNELL)

Scheme for evaluating frequency weight function \(\sderiv{\zeta}{x}\) and integral dimensionless eigenfrequency \(\omega_{\rm int}\); one of:

deps_scheme (default 'MODEL')

Scheme for calculating burning partial derivatives \((\partial\ln\epsilon/\partial\ln T)_{\rho}\) and \((\partial\ln\epsilon/\partial\ln\rho)_{T}\); one of

  • 'MODEL' : Use values from model
  • 'FILE' : Use complex (phase-lagged) values from separate file
deps_file (default '')
Name of epsilon partial derivatives file, when deps_scheme is 'FILE'
deps_file_format (default 'WOLF')

Format of epsilon partial derivative file, when deps_scheme is 'FILE'; one of:

x_ref (default 1 or outer grid point, whichever is smaller)
Reference fractional radius for photosphere, normalizations etc.
nonadiabatic (default .FALSE.)
Flag to include non-adiabatic effects
quasiad_eigfuncs (default .FALSE.)
Flag to calculate quasi-adiabatic entropy/luminosity eigenfunctions during adiabatic calculations
eddington_approx (default .FALSE.)
Flag to use the Eddington approximation
reduce_order (default .TRUE.)
Flag to reduce the order of the adiabatic radial-pulsation equations from 4 to 2
tag_list (default '', which matches all)
Comma-separated list of &mode tags to match