Legacy Equation Format¶

In this format, equations consists of series of commands that build the numerical system.

Overview¶

In order to write oscillation equations with this format, one need to write them projected onto the spherical harmonics basis.

Equation example:

$\lambda b_m^l = \sum_{l'=|m|}^{\infty} - \iint_{4\pi} \{Y_l^m\}^*Y_{l'}^m d\Omega \partial_\zeta u_m^{l'} - \iint_{4\pi} \frac{2 \zeta H + \zeta^2 N H_\zeta}{r^2 r_\zeta} \{Y_l^m\}^*Y_{l'}^m d\Omega u_m^{l'} + ...$

In order to write such an equation in TOP, one need to split it into terms and add them incrementally in the numerical system. The purpose of TOP, in particular top-build is to ease writing such equations. Several commands can help inserting such terms in the numerical system, see section Commands for a list of commands available in TOP’s language.

Few other features can help understanding how to write equations in TOP: for a description of the type of terms TOP can handle, see section Type of Terms. And see section String pre-processing for a description of string pre-processing and special variables.

Type of Terms¶

For each term to be added in the numerical system, one need to provide TOP with its type. The type of a term specify whether it depends on nothing (scalar), the radial coordinate $$r$$, $$l$$ or $$l'$$.

The type of terms currently supported by TOP are:

s: scalar term term only depending on $$r$$ (or $$\zeta$$) (only available for 1D equations) term depending on $$l$$ and $$l'$$ (only available in 2D equations) term depending on $$r$$ (or $$\zeta$$) and $$l$$ (only in 2D equations) term depending on $$r$$ (or $$\zeta$$), $$l$$ and $$l'$$ (only available in 2D equations)

String pre-processing¶

$a: stands for the current coupling matrix. This variable can only be used in a term definition the previous coupling matrix. It correspond to variable$a of the last term definition. the array containing $$l$$ values of the current equation. the array containing $$l$$ values of the current variable. is the $$x^{th}$$ $$l$$ value of the current equation. is the $$x^{th}$$ $$l$$ value of the current variable. index of the current equation. index of the current variable. radial resolution of the problem. indices of the radial coordinate (this will result in the generation of a FORTRAN loop over all radial points). indices of the horizontal coordinate $$l$$ (this will result in the generation of a FORTRAN loop over all values of $$l$$). indices of the horizontal coordinate $$l'$$ (this will result in the generation of a FORTRAN loop over all values of $$l'$$).

Commands¶

input¶

This command allow the user to define a parameter for the set of equation to be written.

Syntax: input parameter format parameter: name of the parameter. format: FORTRAN format specifier. This is used to read the parameter input file. input mass 0pf5.2

stamp¶

Used to define a string to appear in the output files.

Syntax: stamp string string: string to appear in output files. stamp eq_ESTER_all_lagrange

definition¶

Used to define named constants.

Syntax: definition type name value type: type of the constant (integer, double_precision, complex) name: name of the constant value: value of the constant definition double_precision gamma_p 1d0 + 1d0/pindex will define gamma_p with a value of $$1+\frac{1}{pindex}$$ (where $$pindex$$ has to be another variable (definition or input)

eqlist¶

Defines the name of equation to be defined in the equation file.

Syntax: eqlist eq1 eq2 eq3 ... # and so on eq1: name of the equation eqlist eqEr eqdP eqPhi eqPhiP defines 4 equations named eqEr, eqdP, eqPhi and epPhiP

varlist¶

Defines variables of the equation set.

Syntax: varlist var1 var2 var3 ... # and so on var1: name of the variable varlist Er dP Phi PhiP defines 4 variables named Er, dP, Phi and PhiP

leq¶

In TOP equation are projected into the spherical harmonic basis. This command is use to define the starting $$l$$ for this projection.

Syntax: leq eqName value eqName: name of the equation value: starting value of $$l$$ leq eqEr abs(m)+iparity

lvar¶

In TOP variables are projected into the spherical harmonic basis. This command is use to define the starting $$l$$ for this projection.

Syntax: lvar varName value varName: name of the variable value: starting value of $$l$$ lvar Er abs(m)+iparity

equation¶

This command is used to start defining an equation. This means that further command in the equation file will apply to the current equation.

Syntax: equation eqName eqName: name of the equation equation eqEr

sub¶

This is use to insert a term in the current equation: this term will be computed by calling a FORTRAN subroutine.

term¶

Used to insert a term in the equation.

Syntax: term type power expression  variable type: the type of term see type of terms in TOP. power: the power of the eigenvalue preceded by a w. expression: the mathematical expression of the term to be inserted. variable: name of the variable involved in the coupling. Further characters can be used indicate radial derives. For instance, Er' mean $$\frac{\partial Er}{\partial r}$$. Higher derivative order can be achieved either by chaining the ' character or with the ^ character followed by the derivative order: Er^2 is equivalent to Er''. term s w0 -2d0 Pi'': this would insert the term $$-2\frac{\partial^2 Pi}{\partial r^2}$$ in the current equation.

termbc¶

Used to insert a term in a boundary condition of the system.

Syntax: termbc type location power expression  variable(index) type: the type of term see type of terms in TOP. power: the power of the eigenvalue preceded by a w. location: the location where the boundary condition should be inserted in the numerical system. This is basically tells the line in the matrix to be replaced with the boundary condition. expression: the mathematical expression of the term to be inserted. variable: name of the variable involved in the coupling. Further characters can be used indicate radial derives. For instance, Er' mean $$\frac{\partial Er}{\partial r}$$. Higher derivative order can be achieved either by chaining the ' character or with the ^ character followed by the derivative order: Er^2 is equivalent to Er''. index' radial coordinate of the boundary condition. termbc t $nr w0 1d0 Phi($nr)': this would insert the term $$\Phi(r=surf)$$ in the boundary condition. (The last line of the matrix would be replaced with this boundary condition).

instruction¶

Used to add ad-hoc FORTRAN instruction in the module responsible for computing coupling integrals.

Syntax: instruction fortran fortran: the FORTRAN instruction to be inserted. instruction call modify_l0($prev,$nr, abs(m)+iparity): will insert the code call modify_l0(dm(1)%artt(:, :, :), grd(1)%nr, abs(m)+iparity). See String pre-processing.