Output Diagnostics
vel_energy.dat
| Column |
Quantity |
Description |
More information |
| 1 |
time |
magnetic diffusion time |
|
| 2 |
Kinetic Energy |
|
Definition |
| 3 |
Total Toroidal |
|
|
| 4 |
E_eq_sym |
|
|
| 5 |
E_l=1 |
|
|
| 6 |
E_tor_m=0 |
|
|
| 7 |
Edis |
|
Definition |
| 8 |
Axisymmetric KE |
|
|
mag_energy.dat
| Column |
Quantity |
Description |
More information |
| 1 |
time |
magnetic diffusion time |
|
| 2 |
Magnetic Energy |
|
Definition |
| 3 |
Total Toroidal |
|
|
| 4 |
E_eq_sym |
|
|
| 5 |
E_l=1 |
|
|
| 6 |
E_tor_m=0 |
|
|
| 7 |
Mag. Dissipation |
|
Definition |
| 8 |
Axial Dipole g10 |
|
|
| 9 |
Axial Quadrupole g20 |
|
|
| 10 |
Axisymmetric ME |
|
|
vel_partition.dat
| Column |
Quantity |
Description |
| 1 |
time |
|
| 2 |
Total Toroidal |
|
| 3 |
Total Poloidal |
|
| 4 |
axisym. tor. |
|
| 5 |
axisym. pol. |
|
| 6 |
equat.sym.tor. |
|
| 7 |
equat.sym.pol. |
|
| 8 |
equat+axisym.tor. |
\(E_\text{kin}^\text{Tor,eq+ax} = \sum_{l \text{ even}} E^\text{Tor}_\text{kin}(l, m=0)\) |
| 9 |
equat.axisym.pol. |
\(E_\text{kin}^\text{Pol,eq+ax} = \sum_{l \text{ even}} E^\text{Pol}_\text{kin}(l, m=0)\) |
| 10 |
sum(l x u(l).u(l)) |
|
| 11 |
lpol |
|
mag_partition.dat
| Column |
Quantity |
Description |
| 1 |
time |
|
| 2 |
Total Toroidal |
|
| 3 |
Total Poloidal |
|
| 4 |
axisym. tor. |
|
| 5 |
axisym. pol. |
|
| 6 |
equat.sym.tor. |
\(E_\text{mag}^\text{Tor,eq+ax} = \sum_{l \text{ even}} E^\text{Tor}_\text{mag}(l, m=0)\) |
| 7 |
equat.sym.pol. |
\(E_\text{mag}^\text{Pol,eq+ax} = \sum_{l \text{ even}} E^\text{Pol}_\text{mag}(l, m=0)\) |
| 8 |
equat+axisym.tor. |
|
| 9 |
equat.axisym.pol. |
|
mag_cmb.dat
| Column |
Quantity |
Description |
More information |
| 1 |
time |
|
|
| 2 |
Surf Mag Energy |
|
|
| 3 |
Even mag Energy |
|
|
| 4 |
l=1 surf mag Energy |
|
|
| 5 |
N Mag pole theta |
|
|
| 6 |
N Mag pole phi |
|
|
| 7 |
Axisym Surf mag en. |
|
|
| 8 |
l<=12 surf mag Ene |
|
|
dynamo_outputs.dat
cod_nusselt.dat
Nusselt number calculations performed assuming conductive profile of
\[
T^\text{cond} = \frac{-A}{r} + B,~~\frac{\partial T^\text{cond}}{\partial r} = \frac{A}{r^2},
\]
where \(A=(r_i\times r_o)/(r_i-r_o)\), \(~B=-r_i~\) for fixed temperature (\(\Delta T=1\)), and
\(~A = -1\), \(~B=0\) for fixed flux.
| Column |
Quantity |
Description |
| 1 |
time |
|
| 2 |
Nusselt |
if FF: Nu_flux, if FT: Nu_top |
| 3 |
T_top |
\(T\rvert_{r_o}\) |
| 4 |
T_bot |
\(T\rvert_{r_i}\) |
| 5 |
dCdr_ri |
\(\partial T/\partial r \rvert_{r_i}\) |
| 6 |
dCdr_ro |
\(\partial T/\partial r \rvert_{r_o}\) |
| 7 |
Nu_top |
\((\partial T/\partial r \rvert_{r_o})/ (\partial T^\text{cond}/\partial r\rvert_{r_o})\) |
| 8 |
Nu_bot |
\((\partial T/\partial r \rvert_{r_i})/ (\partial T^\text{cond}/\partial r\rvert_{r_i})\) |
| 9 |
Nu_flux |
\(\Delta T^\text{cond}/(T\rvert_{r_i}-T\rvert_{r_o})\) |
comp_nusselt.dat
Sherwood number calculations performed assuming conductive profile of
\[
\frac{\partial \xi^\text{cond}}{\partial r} = \frac{A}{r^2}, ~~ A=\frac{r_i\times r_o}{r_i-r_o}.
\]
Column headers remain as for Thermal nusselt number for compatibility.
| Column |
Quantity |
Description |
| 1 |
time |
|
| 2 |
Nusselt |
if FF: Nu_flux, if FT: Nu_top |
| 3 |
T_top |
\(\xi\rvert_{r_o}\) |
| 4 |
T_bot |
\(\xi\rvert_{r_i}\) |
| 5 |
dCdr_ri |
\(\partial \xi/\partial r \rvert_{r_i}\) |
| 6 |
dCdr_ro |
\(\partial \xi/\partial r \rvert_{r_o}\) |
| 7 |
Nu_top |
\((\partial \xi/\partial r \rvert_{r_o})/ (\partial \xi^\text{cond}/\partial r\rvert_{r_o})\) |
| 8 |
Nu_bot |
\((\partial \xi/\partial r \rvert_{r_i})/ (\partial \xi^\text{cond}/\partial r\rvert_{r_i})\) |
| 9 |
Nu_flux |
\(\Delta \xi^\text{cond}/(\xi\rvert_{r_i}-\xi\rvert_{r_o})\) |
rot_torque.dat
mag_spectrum
One table produced every i_save_rate3 timesteps. The first half of this table contains data for harmonic degree \(l\). A two-line space at the end is followed by the data for harmonic order \(m\).
| Column |
Quantity |
Description |
More information |
| 1 |
l |
|
|
| 2 |
E(l) |
|
|
| 3 |
Epol(l) |
|
|
| 4 |
Etor(m) |
|
|
| --- |
--- |
--- |
--- |
| 1 |
m |
|
|
| 2 |
E(m) |
|
|
| 3 |
Epol(m) |
|
|
| 4 |
Etor(m) |
|
|
vel_spectrum
One table produced every i_save_rate3 timesteps.
The first half of this table contains data for harmonic degree \(l\).
A two-line space at the end is followed by the data for harmonic order \(m\).
| Column |
Quantity |
Description |
More information |
| 1 |
l |
|
|
| 2 |
E(l) |
|
|
| 3 |
Epol(l) |
|
|
| 4 |
Etor(m) |
|
|
| --- |
--- |
--- |
--- |
| 1 |
m |
|
|
| 2 |
E(m) |
|
|
| 3 |
Epol(m) |
|
|
| 4 |
Etor(m) |
|
|
cod_spectrum
Before the dynamo code was capable of doubly-diffusive convection, only one buoyancy source was present, called the codensity.
This buoyancy field is set as the thermal buoyancy source for double diffusive convection.
The first half of this table contains data for harmonic degree \(l\).
A two-line space at the end is followed by the data for harmonic order \(m\).
| Column |
Quantity |
Description |
More information |
| 1 |
l |
|
|
| 2 |
cod(l) |
|
|
| --- |
--- |
--- |
--- |
| 1 |
m |
|
|
| 2 |
cod(m) |
|
|
com_spectrum
The first half of this table contains data for harmonic degree \(l\).
A two-line space at the end is followed by the data for harmonic order \(m\).
| Column |
Quantity |
Description |
More information |
| 1 |
l |
|
|
| 2 |
com(l) |
|
|
| --- |
--- |
--- |
--- |
| 1 |
m |
|
|
| 2 |
com(m) |
|
|
tim_step.dat
temp_profile
| Column |
Quantity |
Description |
More information |
| 1 |
r |
|
|
| 2 |
C |
|
|
| 3 |
dCdr |
|
|
comp_profile
| Column |
Quantity |
Description |
More information |
| 1 |
r |
|
|
| 2 |
C |
|
|
| 3 |
dCdr |
|
|