===================================
List of parameters for ``print.in``
===================================
This page lists 3051 variables distributed into 93 files. Of these:

- 2719 (90%) are documented;
- 332 (10%) are undocumented.

Some of the variable names are shared amongst modules, so there are 1703 unique names:

- 912 variables appear 1 time;
- 373 variables appear 2 times;
- 312 variables appear 3 times;
- 79 variables appear 4 times;
- 21 variables appear 5 times;
- 6 variables appear 6 times.

.. raw:: html

   <div>Filter: <input type='text' id='customvarsearch' /></div><br/>


Module *bfield.f90*
~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *bmax*
     - :math:`\max B`
   * - *bmin*
     - :math:`\min B`
   * - *brms*
     - :math:`\langle B^2\rangle^{1/2}`
   * - *bm*
     - :math:`\langle B\rangle`
   * - *b2m*
     - :math:`\langle B^2\rangle`
   * - *bxmax*
     - :math:`\max|B_x|`
   * - *bymax*
     - :math:`\max|B_y|`
   * - *bzmax*
     - :math:`\max|B_z|`
   * - *bxm*
     - :math:`\langle B_x\rangle`
   * - *bym*
     - :math:`\langle B_y\rangle`
   * - *bzm*
     - :math:`\langle B_z\rangle`
   * - *bx2m*
     - :math:`\langle B_x^2\rangle`
   * - *by2m*
     - :math:`\langle B_y^2\rangle`
   * - *bz2m*
     - :math:`\langle B_z^2\rangle`
   * - *bxbym*
     - :math:`\langle B_x B_y\rangle`
   * - *bxbzm*
     - :math:`\langle B_x B_z\rangle`
   * - *bybzm*
     - :math:`\langle B_y B_z\rangle`
   * - *dbxmax*
     - :math:`\max|B_x - B_{\mathrm{ext,}x}|`
   * - *dbymax*
     - :math:`\max|B_y - B_{\mathrm{ext,}y}|`
   * - *dbzmax*
     - :math:`\max|B_z - B_{\mathrm{ext,}z}|`
   * - *dbxm*
     - :math:`\langle B_x - B_{\mathrm{ext,}x}\rangle`
   * - *dbym*
     - :math:`\langle B_y - B_{\mathrm{ext,}y}\rangle`
   * - *dbzm*
     - :math:`\langle B_z - B_{\mathrm{ext,}z}\rangle`
   * - *dbx2m*
     - :math:`\langle\left(B_x - B_{\mathrm{ext,}x}\right)^2\rangle`
   * - *dby2m*
     - :math:`\langle\left(B_y - B_{\mathrm{ext,}y}\right)^2\rangle`
   * - *dbz2m*
     - :math:`\langle\left(B_z - B_{\mathrm{ext,}z}\right)^2\rangle`
   * - *jmax*
     - :math:`\max J`
   * - *jmin*
     - :math:`\min J`
   * - *jrms*
     - :math:`\langle J^2\rangle^{1/2}`
   * - *jm*
     - :math:`\langle J\rangle`
   * - *j2m*
     - :math:`\langle J^2\rangle`
   * - *jxmax*
     - :math:`\max|J_x|`
   * - *jymax*
     - :math:`\max|J_y|`
   * - *jzmax*
     - :math:`\max|J_z|`
   * - *jxm*
     - :math:`\langle J_x\rangle`
   * - *jym*
     - :math:`\langle J_y\rangle`
   * - *jzm*
     - :math:`\langle J_z\rangle`
   * - *jx2m*
     - :math:`\langle J_x^2\rangle`
   * - *jy2m*
     - :math:`\langle J_y^2\rangle`
   * - *jz2m*
     - :math:`\langle J_z^2\rangle`
   * - *divbmax*
     - :math:`\max|\nabla\cdot\vec{B}|`
   * - *divbrms*
     - :math:`\langle\left(\nabla\cdot\vec{B}\right)^2\rangle^{1/2}`
   * - *betamax*
     - :math:`\max\beta`
   * - *betamin*
     - :math:`\min\beta`
   * - *betam*
     - :math:`\langle\beta\rangle`
   * - *vAmax*
     - :math:`\max v_A`
   * - *vAmin*
     - :math:`\min v_A`
   * - *vAm*
     - :math:`\langle v_A\rangle`

Module *cdata.f90*
~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *it*
     - number of time step :math:`\quad`  (since beginning of job only)
   * - *t*
     - time :math:`t`  :math:`\quad`  (since start.csh)
   * - *dt*
     - time step :math:`\delta t`
   * - *walltime*
     - wall clock time since start of run.x, in seconds
   * - *timeperstep*
     -
   * - *dtv*
     - advective timestep as a fraction of the actual one
   * - *dtdiffus*
     - diffusive timestep as a fraction of the actual one
   * - *dtdiffus2*
     - hyperdiffusive (hyper2)  timestep as a fraction of the actual one
   * - *dtdiffus3*
     - hyperdiffusive (hyper3)  timestep as a fraction of the actual one
   * - *Rmesh*
     - :math:`R_{\mathrm mesh}`
   * - *Rmesh3*
     - :math:`R_{\mathrm mesh}^{(3)}`
   * - *maxadvec*
     - maxadvec
   * - *eps_rkf*
     - time step accuracy threshold

Module *chemistry.f90*
~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Ym*
     - :math:`\left<Y_x\right>`
   * - *rhoYm*
     - :math:`\left<\rho Y_x\right>`
   * - *TYm*
     - :math:`\left<Y_{\mathrm thresh}-Y_x\right>`
   * - *dYm*
     - :math:`\delta\left<Y_x\right>/\delta t`
   * - *dYmax*
     - :math:`max\delta\left<Y_x\right>/\delta t`
   * - *Ymax*
     - :math:`\left<Y_{x,max}\right>`
   * - *Ymin*
     - :math:`\left<Y_{x,min}\right>`
   * - *hm*
     - :math:`\left<H_{x,max}\right>`
   * - *cpm*
     - :math:`\left<c_{p,x}\right>`
   * - *diffm*
     - :math:`\left<D_{x}\right>`
   * - *diffmax*
     - :math:`\left<D_{x,max}\right>`
   * - *diffmin*
     - :math:`\left<D_{x,min}\right>`
   * - *Ymz*
     - :math:`\left<Y_x\right>_{xy}(z)`
   * - *dtchem*
     - :math:`dt_{chem}`
   * - *nuclrmin*
     - :math:`\left< r_{\min} \right>`
   * - *nuclrate*
     - :math:`\left< J \right>`

Module *chemistry_simple.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Ym*
     - :math:`\left<Y_x\right>`
   * - *dYm*
     - :math:`\delta\left<Y_x\right>/\delta t`
   * - *dYmax*
     - :math:`max\delta\left<Y_x\right>/\delta t`
   * - *Ymax*
     - :math:`\left<Y_{x,max}\right>`
   * - *Ymin*
     - :math:`\left<Y_{x,min}\right>`
   * - *hm*
     - :math:`\left<H_{x,max}\right>`
   * - *cpm*
     - :math:`\left<c_{p,x}\right>`
   * - *diffm*
     - :math:`\left<D_{x}\right>`
   * - *Ymz*
     - :math:`\left<Y_x\right>_{xy}(z)`
   * - *dtchem*
     - :math:`dt_{chem}`

Module *density.f90*
~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rhom*
     - :math:`\left<\varrho\right>`  :math:`\quad`  (mean density)
   * - *rhomxmask*
     - :math:`\left<\varrho\right>`  for the density_xaver_range
   * - *rhomzmask*
     - :math:`\left<\varrho\right>`  for the density_zaver_range
   * - *rho2m*
     - :math:`\left<\varrho^2\right>`
   * - *rho4m*
     - :math:`\left<\varrho^4\right>`
   * - *rho6m*
     - :math:`\left<\varrho^6\right>`
   * - *rho8m*
     - :math:`\left<\varrho^8\right>`
   * - *rho12m*
     - :math:`\left<\varrho^{12}\right>`
   * - *rhof2m*
     - :math:`\left<\varrho'^2\right>`
   * - *lnrho2m*
     -
   * - *drho2m*
     - :math:`<(\varrho-\varrho_0)^2>`
   * - *drhom*
     - :math:`<\varrho-\varrho_0>`
   * - *rhomin*
     - :math:`\min(\rho)`
   * - *rhomax*
     - :math:`\max(\rho)`
   * - *rhominloc*
     - :math:`location of \min(\rho)`
   * - *rhomaxloc*
     - :math:`location of \max(\rho)`
   * - *lnrhomin*
     - :math:`\min(\log\rho)`
   * - *lnrhomax*
     - :math:`\max(\log\rho)`
   * - *rhorms*
     - :math:`\sqrt{<\varrho^2>}`
   * - *lnrhorms*
     - :math:`\sqrt{<(\ln\varrho)^2>}`
   * - *ugrhom*
     - :math:`\left<\uv\cdot\nabla\varrho\right>`
   * - *uglnrhom*
     - :math:`\left<\uv\cdot\nabla\ln\varrho\right>`
   * - *dtd*
     -
   * - *dtd3*
     -
   * - *rhomr*
     -
   * - *totmass*
     - :math:`\int\varrho\,dV`
   * - *mass*
     - :math:`\int\varrho\,dV`
   * - *sphmass*
     - :math:`\int\varrho\,dV`  inside :math:`r < r_{\mathrm diag}` .
   * - *inertiaxx*
     - :math:`xx`  component of the inertia tensor (spherical coordinates)
   * - *inertiayy*
     - :math:`yy`  component of the inertia tensor (spherical coordinates)
   * - *inertiazz*
     - :math:`zz`  component of the inertia tensor (spherical coordinates)
   * - *inertiaxx_car*
     - :math:`xx`  component of the inertia tensor (Cartesian coordinates)
   * - *inertiayy_car*
     - :math:`xx`  component of the inertia tensor (Cartesian coordinates)
   * - *inertiazz_car*
     - :math:`xx`  component of the inertia tensor (Cartesian coordinates)
   * - *vol*
     - :math:`\int\,dV`  (volume)
   * - *grhomax*
     - :math:`\max (|\nabla \varrho|)`
   * - *kap_tdep*
     - time-dependent diffusivity

Module *density_stratified.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *mass*
     - :math:`\int\rho\,d^3x`
   * - *rhomin*
     - :math:`\min\left|\rho\right|`
   * - *rhomax*
     - :math:`\max\left|\rho\right|`
   * - *drhom*
     - :math:`\langle\Delta\rho/\rho_0\rangle`
   * - *drho2m*
     - :math:`\langle\left(\Delta\rho/\rho_0\right)^2\rangle`
   * - *drhorms*
     - :math:`\langle\Delta\rho/\rho_0\rangle_{rms}`
   * - *drhomax*
     - :math:`\max\left|\Delta\rho/\rho_0\right|`

Module *detonate.f90*
~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *detn*
     - Number of detonated sites (summed over time steps between adjacent outputs)
   * - *dettot*
     - Total energy input (summed over time steps between adjacent outputs)

Module *dustdensity.f90*
~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *KKm*
     - :math:`\sum {\cal T}_k^{\mathrm coag}`
   * - *KK2m*
     - :math:`\sum {\cal T}_k^{\mathrm coag}`
   * - *MMxm*
     - :math:`\sum {\cal M}^x_{k,{\mathrm coag}}`
   * - *MMym*
     - :math:`\sum {\cal M}^y_{k,{\mathrm coag}}`
   * - *MMzm*
     - :math:`\sum {\cal M}^z_{k,{\mathrm coag}}`

Module *entropy.f90*
~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtc*
     - :math:`\delta t/[c_{\delta t}\,\delta_x /\max c_{\mathrm s}]`  :math:`\quad`  (time step relative to acoustic time step; see \S~\ref{time-step})
   * - *ethm*
     - :math:`\left<\varrho e\right>`  :math:`\quad`  (mean thermal [=internal] energy)
   * - *ethdivum*
     -
   * - *ssruzm*
     - :math:`\left<s \varrho u_z/c_p\right>`
   * - *ssuzm*
     - :math:`\left<s u_z/c_p\right>`
   * - *ssm*
     - :math:`\left<s/c_p\right>`  :math:`\quad`  (mean entropy)
   * - *ssbycpm*
     - :math:`\left<s/c_p\right>`  :math:`\quad`  (mean entropy)
   * - *ss2m*
     - :math:`\left<(s/c_p)^2\right>`  :math:`\quad`  (mean squared entropy) :math:`\quad`  (mean of the running average of entropy)
   * - *eem*
     - :math:`\left<e\right>`
   * - *ppm*
     - :math:`\left<p\right>`
   * - *ppmax*
     - :math:`\max(p)`
   * - *ppmin*
     - :math:`\min(p)`
   * - *csm*
     - :math:`\left<c_{\mathrm s}\right>`
   * - *csmax*
     - :math:`\max (c_{\mathrm s})`
   * - *cgam*
     - :math:`\left<c_{\gamma}\right>`
   * - *pdivum*
     - :math:`\left<p\nabla\cdot\uv\right>`
   * - *heatm*
     -
   * - *ugradpm*
     -
   * - *fradbot*
     - :math:`\int F_{\mathrm bot}\cdot d\vec{S}`
   * - *fradtop*
     - :math:`\int F_{\mathrm top}\cdot d\vec{S}`
   * - *TTtop*
     - :math:`\int T_{\mathrm top} d\vec{S}`
   * - *ethtot*
     - :math:`\int_V\varrho e\,dV`  :math:`\quad`  (total thermal [=internal] energy)
   * - *dtchi*
     - :math:`\delta t / [c_{\delta t,{\mathrm v}}\, \delta x^2/\chi_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat conductivity; see \S~\ref{time-step})
   * - *Hmax*
     - :math:`H_{\mathrm max}` :math:`\quad`  (net heat sources summed see \S~\ref{time-step})
   * - *tauhmin*
     - :math:`\min(\tau_{\mathrm heat})`
   * - *dtH*
     - :math:`\delta t / [c_{\delta t,{\mathrm s}}\, c_{\mathrm v}T /H_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat sources; see \S~\ref{time-step})
   * - *yHm*
     - mean hydrogen ionization
   * - *yHmax*
     - max of hydrogen ionization
   * - *TTm*
     - :math:`\left<T\right>`
   * - *TTmax*
     - :math:`T_{\max}`
   * - *TTmin*
     - :math:`T_{\min}`
   * - *TTmaxloc*
     - location of :math:`T_{\max}`
   * - *TTminloc*
     - location of :math:`T_{\min}`
   * - *gTmax*
     - :math:`\max (|\nabla T|)`
   * - *ssmax*
     - :math:`s_{\max}`
   * - *ssmin*
     - :math:`s_{\min}`
   * - *gTrms*
     - :math:`(\nabla T)_{\mathrm rms}`
   * - *gsrms*
     - :math:`(\nabla s)_{\mathrm rms}`
   * - *gTxgsrms*
     - :math:`(\nabla T\times\nabla s)_{\mathrm rms}`
   * - *gTxgsom*
     - :math:`\left<  (\nabla T\times\nabla s)\cdot\boldsymbol{\omega} \right>`
   * - *fconvm*
     - :math:`\left< c_p \varrho u_z T \right>`
   * - *TTp*
     -
   * - *ssmr*
     -
   * - *TTmr*
     -
   * - *ufpresm*
     - :math:`\left< -u/\rho\nabla p \right>`
   * - *Kkramersm*
     - :math:`\left< K_{\mathrm kramers} \right>`
   * - *chikrammin*
     - :math:`\min (\chi_{\mathrm kramers})`
   * - *chikrammax*
     - :math:`\max (\chi_{\mathrm kramers})`
   * - *TT2m*
     - :math:`\left<(T)^2\right>`  :math:`\quad`  (mean squared temperature)

Module *forcing.f90*
~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *bfm*
     - :math:`\left<\Bv\cdot\fv\right>`
   * - *jfm*
     - :math:`\left<\Jv\cdot\fv\right>`
   * - *rufm*
     - :math:`\left<\rho\fv\cdot\uv\right>`
   * - *ufm*
     - :math:`\left<\fv\cdot\uv\right>`
   * - *ofm*
     - :math:`\left<\ov\cdot\fv\right>`

Module *gravity_simple.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *epot*
     - :math:`\left<\varrho \Phi_{\mathrm grav} \right>`  :math:`\quad`  (mean potential energy)
   * - *epottot*
     - :math:`\int_V\varrho \Phi_{\mathrm grav} dV`  :math:`\quad`  (total potential energy)
   * - *ugm*
     - :math:`\left<\uv \cdot \gv\right>`
   * - *rugm*
     - :math:`\left<\varrho \uv \cdot \gv\right>`
   * - *rgxm*
     - :math:`\left<\varrho g_x \right>`
   * - *Wgrav*
     - :math:`\int\varrho \uv \cdot \gv \, dV`
   * - *Fgravx*
     - :math:`\int\varrho g_x \, dV`

Module *heatflux.f90*
~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtspitzer*
     - Spitzer heat conduction time step
   * - *dtq*
     - heatflux time step
   * - *dtq2*
     - heatflux time step due to tau
   * - *qmax*
     - :math:`\max(|\qv|)`
   * - *tauqmax*
     - :math:`\max(|\tau_{\mathrm Spitzer}|)`
   * - *qxmin*
     - :math:`\min(|q_x|)`
   * - *qymin*
     - :math:`\min(|q_y|)`
   * - *qzmin*
     - :math:`\min(|q_z|)`
   * - *qxmax*
     - :math:`\max(|q_x|)`
   * - *qymax*
     - :math:`\max(|q_y|)`
   * - *qzmax*
     - :math:`\max(|q_z|)`
   * - *qrms*
     - :math:`\sqrt(|\qv|^2)`
   * - *qsatmin*
     - minimum of qsat/qabs
   * - *qsatrms*
     - rms of qsat/abs

Module *hydro.f90*
~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *u2tm*
     - :math:`\left<\uv(t)\cdot\int_0^t\uv(t') dt'\right>`
   * - *uotm*
     - :math:`\left<\uv(t)\cdot\int_0^t\omv(t') dt'\right>`
   * - *outm*
     - :math:`\left<\omv(t)\cdot\int_0^t\uv(t') dt'\right>`
   * - *fkinzm*
     - :math:`\left<{1\over2} \varrho\uv^2 u_z\right>`
   * - *gamm*
     - :math:`\left<gamma\right>`
   * - *gamrms*
     - :math:`\left<\gamma^2\right>^{1/2}`
   * - *gammax*
     - :math:`\max(\gamma)`
   * - *u2m*
     - :math:`\left<\uv^2\right>`
   * - *u2sphm*
     - :math:`\int_{r=0}^{r=r_{\mathrm diag}} \uv^2 dV` , where :math:`r=\sqrt{x^2+y^2+z^2}`
   * - *um2*
     -
   * - *uxpt*
     - :math:`u_x(x_1,y_1,z_1,t)`
   * - *uypt*
     - :math:`u_y(x_1,y_1,z_1,t)`
   * - *uzpt*
     - :math:`u_z(x_1,y_1,z_1,t)`
   * - *uxp2*
     - :math:`u_x(x_2,y_2,z_2,t)`
   * - *uyp2*
     - :math:`u_y(x_2,y_2,z_2,t)`
   * - *uzp2*
     - :math:`u_z(x_2,y_2,z_2,t)`
   * - *uxuypt*
     - :math:`(u_x u_y)(x_1,y_1,z_1,t)`
   * - *uyuzpt*
     - :math:`(u_y u_z)(x_1,y_1,z_1,t)`
   * - *uzuxpt*
     - :math:`(u_z u_x)(x_1,y_1,z_1,t)`
   * - *urms*
     - :math:`\left<\uv^2\right>^{1/2}`
   * - *urmsx*
     - :math:`\left<\uv^2\right>^{1/2}`  for the hydro_xaver_range
   * - *urmsz*
     - :math:`\left<\uv^2\right>^{1/2}`  for the hydro_zaver_range
   * - *durms*
     - :math:`\left<\delta\uv^2\right>^{1/2}`
   * - *umax*
     - :math:`\max(|\uv|)`
   * - *umin*
     - :math:`\min(|\uv|)`
   * - *uxrms*
     - :math:`\left<u_x^2\right>^{1/2}`
   * - *uyrms*
     - :math:`\left<u_y^2\right>^{1/2}`
   * - *uzrms*
     - :math:`\left<u_z^2\right>^{1/2}`
   * - *uzrmaxs*
     -
   * - *uxmin*
     - :math:`\min(|u_x|)`
   * - *uymin*
     - :math:`\min(|u_y|)`
   * - *uzmin*
     - :math:`\min(|u_z|)`
   * - *uxmax*
     - :math:`\max(|u_x|)`
   * - *uymax*
     - :math:`\max(|u_y|)`
   * - *uzmax*
     - :math:`\max(|u_z|)`
   * - *uxm*
     - :math:`\left<u_x\right>`
   * - *uym*
     - :math:`\left<u_y\right>`
   * - *uzm*
     - :math:`\left<u_z\right>`
   * - *uzcx10m*
     - :math:`\left<u_z\cos10x\right>`
   * - *uzsx10m*
     - :math:`\left<u_z\sin10x\right>`
   * - *uduum*
     - :math:`\left<\boldsymbol{u}\cdot\boldsymbol{u}\right>`
   * - *ux2m*
     - :math:`\left<u_x^2\right>`
   * - *uy2m*
     - :math:`\left<u_y^2\right>`
   * - *uz2m*
     - :math:`\left<u_z^2\right>`
   * - *ux3m*
     - :math:`\left<u_x^3\right>`
   * - *uy3m*
     - :math:`\left<u_y^3\right>`
   * - *uz3m*
     - :math:`\left<u_z^3\right>`
   * - *ux4m*
     - :math:`\left<u_x^4\right>`
   * - *uy4m*
     - :math:`\left<u_y^4\right>`
   * - *uz4m*
     - :math:`\left<u_z^4\right>`
   * - *u4m*
     - :math:`\left<u^4\right>`
   * - *u6m*
     - :math:`\left<u^6\right>`
   * - *u8m*
     - :math:`\left<u^8\right>`
   * - *uxuy2m*
     - :math:`\left<u_x^2u_y^2\right>`
   * - *uyuz2m*
     - :math:`\left<u_y^2u_z^2\right>`
   * - *uzux2m*
     - :math:`\left<u_z^2u_x^2\right>`
   * - *velxx2m*
     - :math:`\left<w \gamma^2 u_x^2\right>`
   * - *velxy2m*
     - :math:`\left<w \gamma^2 u_y^2\right>`
   * - *velxz2m*
     - :math:`\left<w \gamma^2 u_z^2\right>`
   * - *velxrms*
     - :math:`\left<\sqrt{w} \gamma^2 u^2 \right>^{1/2}`
   * - *T00m*
     - :math:`\left< T_{00} \right>`
   * - *Txxm*
     - :math:`\left< T_{xx} \right>`
   * - *Tyym*
     - :math:`\left< T_{yy} \right>`
   * - *Tzzm*
     - :math:`\left< T_{zz} \right>`
   * - *Txym*
     - :math:`\left< T_{xy} \right>`
   * - *Tyzm*
     - :math:`\left< T_{yz} \right>`
   * - *Tzxm*
     - :math:`\left< T_{zx} \right>`
   * - *T0x2m*
     - :math:`\left< T_{0x}^2 \right>`
   * - *T0y2m*
     - :math:`\left< T_{0y}^2 \right>`
   * - *T0z2m*
     - :math:`\left< T_{0z}^2 \right>`
   * - *T0irms*
     - :math:`\left< T_{0i}^2 \right>^{1/2}`
   * - *ux2ccm*
     - :math:`\left<u_x^2\cos^2kz\right>`
   * - *ux2ssm*
     - :math:`\left<u_x^2\sin^2kz\right>`
   * - *uy2ccm*
     - :math:`\left<u_y^2\cos^2kz\right>`
   * - *uy2ssm*
     - :math:`\left<u_y^2\sin^2kz\right>`
   * - *uxuycsm*
     - :math:`\left<u_xu_y\cos kz\sin kz\right>`
   * - *uxuym*
     - :math:`\left<u_x u_y\right>`
   * - *uxuzm*
     - :math:`\left<u_x u_z\right>`
   * - *uyuzm*
     - :math:`\left<u_y u_z\right>`
   * - *umx*
     - :math:`\left< u_x \right>`
   * - *umy*
     - :math:`\left< u_y \right>`
   * - *umz*
     - :math:`\left< u_z \right>`
   * - *omumz*
     - :math:`\left<\left<\Wv\right>_{xy} \cdot\left<\Uv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean cross helicity production)
   * - *umamz*
     - :math:`\left<\left<\uv\right>_{xy}\cdot\left<\Av\right>_{xy}\right>`
   * - *umbmz*
     - :math:`\left<\left<\Uv\right>_{xy} \cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean cross helicity production)
   * - *umxbmz*
     - :math:`\left<\left<\Uv\right>_{xy} \times\left<\Bv\right>_{xy} \right>_z`  :math:`\quad`  (:math:`xy` -averaged mean emf)
   * - *rux2m*
     - :math:`\left<\rho u_x^2\right>`
   * - *ruy2m*
     - :math:`\left<\rho u_y^2\right>`
   * - *ruz2m*
     - :math:`\left<\rho u_z^2\right>`
   * - *divum*
     - :math:`\left<{\mathrm div}\uv)\right>`
   * - *rdivum*
     - :math:`\left<\varrho{\mathrm div}\uv)\right>`
   * - *divu2m*
     - :math:`\left<({\mathrm div}\uv)^2\right>`
   * - *gdivu2m*
     - :math:`\left<({\mathrm grad\,div}\uv)^2\right>`
   * - *u3u21m*
     - :math:`\left<u_3 u_{2,1}\right>`
   * - *u1u32m*
     - :math:`\left<u_1 u_{3,2}\right>`
   * - *u2u13m*
     - :math:`\left<u_2 u_{1,3}\right>`
   * - *u2u31m*
     - :math:`\left<u_2 u_{3,1}\right>`
   * - *u3u12m*
     - :math:`\left<u_3 u_{1,2}\right>`
   * - *u1u23m*
     - :math:`\left<u_1 u_{2,3}\right>`
   * - *u2mr*
     -
   * - *urmr*
     -
   * - *upmr*
     -
   * - *uzmr*
     -
   * - *uxfampm*
     -
   * - *uyfampm*
     -
   * - *uzfampm*
     -
   * - *uxfampim*
     -
   * - *uyfampim*
     -
   * - *uzfampim*
     -
   * - *ruxm*
     - :math:`\left<\varrho u_x\right>`  :math:`\quad`  (mean :math:`x` -momentum density)
   * - *ruym*
     - :math:`\left<\varrho u_y\right>`  :math:`\quad`  (mean :math:`y` -momentum density)
   * - *ruzm*
     - :math:`\left<\varrho u_z\right>`  :math:`\quad`  (mean :math:`z` -momentum density)
   * - *ruxtot*
     - :math:`\left<\rho |u|\right>`  :math:`\quad`  (mean absolute :math:`x` -momentum density)
   * - *rumax*
     - :math:`\max(\varrho |\uv|)`  :math:`\quad`  (maximum modulus of momentum)
   * - *ruxuym*
     - :math:`\left<\varrho u_x u_y\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *ruxuzm*
     - :math:`\left<\varrho u_x u_z\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *ruyuzm*
     - :math:`\left<\varrho u_y u_z\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *divrhourms*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm rms}`
   * - *divrhoumax*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm max}`
   * - *rlxm*
     - :math:`\left< \rho y u_z - z u_y \right>`
   * - *rlym*
     - :math:`\left< \rho z u_x - x u_z \right>`
   * - *rlzm*
     - :math:`\left< \rho x u_y - y u_x \right>`
   * - *rlx2m*
     - :math:`\left<(\rho y u_z-z u_y)^2\right>`
   * - *rly2m*
     - :math:`\left<(\rho z u_x-x u_z)^2\right>`
   * - *rlz2m*
     - :math:`\left<(\rho x u_y-y u_x)^2\right>`
   * - *tot_ang_mom*
     - Total angular momentum in spherical coordinates about the axis.
   * - *dtu*
     - :math:`\delta t/[c_{\delta t}\,\delta x /\max|\mathbf{u}|]`  :math:`\quad`  (time step relative to advective time step; see \S~\ref{time-step})
   * - *oum*
     - :math:`\left<\boldsymbol{\omega}\cdot\uv\right>`
   * - *oxum*
     - :math:`\left<\boldsymbol{\omega}\times\uv\right>`
   * - *ourms*
     - :math:`\left<(\boldsymbol{\omega}\cdot\uv)^2\right>^{1/2}`
   * - *oxurms*
     - :math:`\left<(\boldsymbol{\omega}\times\uv)^2\right>^{1/2}`
   * - *ou_int*
     - :math:`\int_V\boldsymbol{\omega}\cdot\uv\,dV`
   * - *fum*
     - :math:`\left<\fv\cdot\uv\right>`  (continuous forcing only)
   * - *odel2um*
     - :math:`\left<\boldsymbol{\omega}\nabla^2\uv\right>`
   * - *o2m*
     - :math:`\left<\boldsymbol{\omega}^2\right> \equiv \left<(\curl\uv)^2\right>`
   * - *o2u2m*
     - :math:`\left<\uv^2\boldsymbol{\omega}^2\right>`
   * - *o2sphm*
     - :math:`\int_{r=0}^{r=r_{\mathrm diag}} \boldsymbol{\omega}^2 dV` , where :math:`r=\sqrt{x^2+y^2+z^2}`
   * - *orms*
     - :math:`\left<\boldsymbol{\omega}^2\right>^{1/2}`
   * - *omax*
     - :math:`\max(|\boldsymbol{\omega}|)`
   * - *ox2m*
     - :math:`\left<\omega_x^2\right>`
   * - *oy2m*
     - :math:`\left<\omega_y^2\right>`
   * - *oz2m*
     - :math:`\left<\omega_z^2\right>`
   * - *ox3m*
     - :math:`\left<\omega_x^3\right>`
   * - *oy3m*
     - :math:`\left<\omega_y^3\right>`
   * - *oz3m*
     - :math:`\left<\omega_z^3\right>`
   * - *ox4m*
     - :math:`\left<\omega_x^4\right>`
   * - *oy4m*
     - :math:`\left<\omega_y^4\right>`
   * - *oz4m*
     - :math:`\left<\omega_z^4\right>`
   * - *oxm*
     -
   * - *oym*
     -
   * - *ozm*
     -
   * - *oxuzxm*
     - :math:`\left<\omega_x u_{z,x} \right>`
   * - *oyuzym*
     - :math:`\left<\omega_y u_{z,y} \right>`
   * - *oxoym*
     - :math:`\left<\omega_x\omega_y\right>`
   * - *oxozm*
     - :math:`\left<\omega_x\omega_z\right>`
   * - *oyozm*
     - :math:`\left<\omega_y\omega_z\right>`
   * - *qfm*
     - :math:`\left<\qv\cdot\fv\right>`
   * - *q2m*
     - :math:`\left<\qv^2\right>`
   * - *qrms*
     - :math:`\left<\qv^2\right>^{1/2}`
   * - *qmax*
     - :math:`\max(|\qv|)`
   * - *qom*
     - :math:`\left<\qv\cdot\omv\right>`
   * - *quxom*
     - :math:`\left<\qv\cdot(\uv\times\omv)\right>`
   * - *qezxum*
     - :math:`\left< (\boldsymbol{e_z} \times \mathbf{u}) \cdot \mathbf{q} \right>`
   * - *quysm*
     - :math:`\left< \frac{1}{\tau} (u_y^S - u_y) \mathbf{\hat{y}} \cdot \mathbf{q} \right>`
   * - *jxbrqm*
     - :math:`\left<(\Jv\times\Bv/\rho)\cdot\mathbf{q}\right>`
   * - *pvzm*
     - :math:`\left<\omega_z + 2\Omega/\varrho\right>`  :math:`\quad`  (z component of potential vorticity)
   * - *oumphi*
     - :math:`\left<\omv\cdot\uv\right>_\varphi`
   * - *ormphi*
     - :math:`\left<\omega_r\right>_\varphi`
   * - *opmphi*
     - :math:`\left<\omega_\phi\right>_\varphi`
   * - *ozmphi*
     - :math:`\left<\omega_z\right>_\varphi`
   * - *ormr*
     -
   * - *opmr*
     -
   * - *ozmr*
     -
   * - *uguxm*
     -
   * - *uguym*
     -
   * - *uguzm*
     -
   * - *ugurmsx*
     - :math:`\left<\left(\uv\nabla\uv\right)^2\right>^{1/2}`  for the hydro_xaver_range
   * - *gu2m*
     - :math:`\left<(\nabla\uv)^2\right>`
   * - *ugu2m*
     - :math:`\left<\uv\nabla\uv\right>^2`
   * - *dudx*
     - :math:`\left<\frac{\delta \uv}{\delta x}\right>`
   * - *Marms*
     - :math:`\left<\uv^2/\cs^2\right>`  :math:`\quad`  (rms Mach number)
   * - *Mamax*
     - :math:`\max |\uv|/\cs`  :math:`\quad`  (maximum Mach number)
   * - *fintm*
     -
   * - *fextm*
     -
   * - *duxdzma*
     -
   * - *duydzma*
     -
   * - *EEK*
     - :math:`\left<\varrho\uv^2\right>/2`
   * - *EEK2*
     - :math:`\left<(\varrho\uv^2/2)^2\right>`
   * - *EEK3*
     - :math:`\left<(\varrho\uv^2/2)^3\right>`
   * - *EEK4*
     - :math:`\left<(\varrho\uv^2/2)^4\right>`
   * - *ekin*
     - :math:`\left<{1\over2}\varrho\uv^2\right>`
   * - *ekintot*
     - :math:`\int_V{1\over2}\varrho\uv^2\, dV`
   * - *totangmom*
     -
   * - *uxglnrym*
     - :math:`\left<u_x\partial_y\ln\varrho\right>`
   * - *uyglnrxm*
     - :math:`\left<u_y\partial_x\ln\varrho\right>`
   * - *uzdivum*
     - :math:`\left<u_z\nabla\cdot\uv\right>`
   * - *uxuydivum*
     - :math:`\left<u_x u_y\nabla\cdot\uv\right>`
   * - *divuHrms*
     - :math:`(\nabla_{\mathrm H}\cdot\uv_{\mathrm H})^{\mathrm rms}`
   * - *uxxrms*
     - :math:`u_{x,x}^{\mathrm rms}`
   * - *uyyrms*
     - :math:`u_{y,y}^{\mathrm rms}`
   * - *uzzrms*
     - :math:`u_{z,z}^{\mathrm rms}`
   * - *uxzrms*
     - :math:`u_{x,z}^{\mathrm rms}`
   * - *uyzrms*
     - :math:`u_{y,z}^{\mathrm rms}`
   * - *uzyrms*
     - :math:`u_{z,y}^{\mathrm rms}`
   * - *dtF*
     - :math:`\delta t/[c_{\delta t}\,\delta x /\max|\mathbf{F}|]`  :math:`\quad`  (time step relative to max force time step; see \S~\ref{time-step})
   * - *urlm*
     - :math:`\int u_r(\theta,\phi)Y^m_{\ell}(\theta,\phi)\sin(\theta)d\theta d\phi`
   * - *udpxxm*
     - components of symmetric tensor :math:`\left< u_i \partial_j p + u_j \partial_i p \right>`

Module *hydro_kinematic.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *ourms*
     - :math:`\left<(\boldsymbol{\omega}\cdot\uv)^2\right>^{1/2}`
   * - *oxurms*
     - :math:`\left<(\boldsymbol{\omega}\times\uv)^2\right>^{1/2}`
   * - *EEK*
     - :math:`\left<\varrho\uv^2\right>/2`

Module *hydro_potential.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *u2tm*
     - :math:`\left<\uv(t)\cdot\int_0^t\uv(t') dt'\right>`
   * - *fkinzm*
     - :math:`\left<{1\over2} \varrho\uv^2 u_z\right>`
   * - *u2m*
     - :math:`\left<\uv^2\right>`
   * - *um2*
     -
   * - *uxpt*
     - :math:`u_x(x_1,y_1,z_1,t)`
   * - *uypt*
     - :math:`u_y(x_1,y_1,z_1,t)`
   * - *uzpt*
     - :math:`u_z(x_1,y_1,z_1,t)`
   * - *uxp2*
     - :math:`u_x(x_2,y_2,z_2,t)`
   * - *uyp2*
     - :math:`u_y(x_2,y_2,z_2,t)`
   * - *uzp2*
     - :math:`u_z(x_2,y_2,z_2,t)`
   * - *urms*
     - :math:`\left<\uv^2\right>^{1/2}`
   * - *urmsx*
     - :math:`\left<\uv^2\right>^{1/2}`  for the hydro_xaver_range
   * - *urmsz*
     - :math:`\left<\uv^2\right>^{1/2}`  for the hydro_zaver_range
   * - *durms*
     - :math:`\left<\delta\uv^2\right>^{1/2}`
   * - *umax*
     - :math:`\max(|\uv|)`
   * - *umin*
     - :math:`\min(|\uv|)`
   * - *uxrms*
     - :math:`\left<u_x^2\right>^{1/2}`
   * - *uyrms*
     - :math:`\left<u_y^2\right>^{1/2}`
   * - *uzrms*
     - :math:`\left<u_z^2\right>^{1/2}`
   * - *uzrmaxs*
     -
   * - *uxmin*
     - :math:`\min(|u_x|)`
   * - *uymin*
     - :math:`\min(|u_y|)`
   * - *uzmin*
     - :math:`\min(|u_z|)`
   * - *uxmax*
     - :math:`\max(|u_x|)`
   * - *uymax*
     - :math:`\max(|u_y|)`
   * - *uzmax*
     - :math:`\max(|u_z|)`
   * - *uxm*
     - :math:`\left<u_x\right>`
   * - *uym*
     - :math:`\left<u_y\right>`
   * - *uzm*
     - :math:`\left<u_z\right>`
   * - *ux2m*
     - :math:`\left<u_x^2\right>`
   * - *uy2m*
     - :math:`\left<u_y^2\right>`
   * - *uz2m*
     - :math:`\left<u_z^2\right>`
   * - *ux2ccm*
     - :math:`\left<u_x^2\cos^2kz\right>`
   * - *ux2ssm*
     - :math:`\left<u_x^2\sin^2kz\right>`
   * - *uy2ccm*
     - :math:`\left<u_y^2\cos^2kz\right>`
   * - *uy2ssm*
     - :math:`\left<u_y^2\sin^2kz\right>`
   * - *uxuycsm*
     - :math:`\left<u_xu_y\cos kz\sin kz\right>`
   * - *uxuym*
     - :math:`\left<u_x u_y\right>`
   * - *uxuzm*
     - :math:`\left<u_x u_z\right>`
   * - *uyuzm*
     - :math:`\left<u_y u_z\right>`
   * - *umx*
     - :math:`\left< u_x \right>`
   * - *umy*
     - :math:`\left< u_y \right>`
   * - *umz*
     - :math:`\left< u_z \right>`
   * - *omumz*
     - :math:`\left<\left<\Wv\right>_{xy} \cdot\left<\Uv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean cross helicity production)
   * - *umamz*
     - :math:`\left<\left<\uv\right>_{xy}\cdot\left<\Av\right>_{xy}\right>`
   * - *umbmz*
     - :math:`\left<\left<\Uv\right>_{xy} \cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean cross helicity production)
   * - *umxbmz*
     - :math:`\left<\left<\Uv\right>_{xy} \times\left<\Bv\right>_{xy} \right>_z`  :math:`\quad`  (:math:`xy` -averaged mean emf)
   * - *rux2m*
     - :math:`\left<\rho u_x^2\right>`
   * - *ruy2m*
     - :math:`\left<\rho u_y^2\right>`
   * - *ruz2m*
     - :math:`\left<\rho u_z^2\right>`
   * - *divum*
     - :math:`\left<{\mathrm div}\uv)\right>`
   * - *rdivum*
     - :math:`\left<\varrho{\mathrm div}\uv)\right>`
   * - *divu2m*
     - :math:`\left<({\mathrm div}\uv)^2\right>`
   * - *gdivu2m*
     - :math:`\left<({\mathrm grad\,div}\uv)^2\right>`
   * - *u3u21m*
     - :math:`\left<u_3 u_{2,1}\right>`
   * - *u1u32m*
     - :math:`\left<u_1 u_{3,2}\right>`
   * - *u2u13m*
     - :math:`\left<u_2 u_{1,3}\right>`
   * - *u2u31m*
     - :math:`\left<u_2 u_{3,1}\right>`
   * - *u3u12m*
     - :math:`\left<u_3 u_{1,2}\right>`
   * - *u1u23m*
     - :math:`\left<u_1 u_{2,3}\right>`
   * - *u2mr*
     -
   * - *urmr*
     -
   * - *upmr*
     -
   * - *uzmr*
     -
   * - *uxfampm*
     -
   * - *uyfampm*
     -
   * - *uzfampm*
     -
   * - *uxfampim*
     -
   * - *uyfampim*
     -
   * - *uzfampim*
     -
   * - *ruxm*
     - :math:`\left<\varrho u_x\right>`  :math:`\quad`  (mean :math:`x` -momentum density)
   * - *ruym*
     - :math:`\left<\varrho u_y\right>`  :math:`\quad`  (mean :math:`y` -momentum density)
   * - *ruzm*
     - :math:`\left<\varrho u_z\right>`  :math:`\quad`  (mean :math:`z` -momentum density)
   * - *ruxtot*
     - :math:`\left<\rho |u|\right>`  :math:`\quad`  (mean absolute :math:`x` -momentum density)
   * - *rumax*
     - :math:`\max(\varrho |\uv|)`  :math:`\quad`  (maximum modulus of momentum)
   * - *ruxuym*
     - :math:`\left<\varrho u_x u_y\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *ruxuzm*
     - :math:`\left<\varrho u_x u_z\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *ruyuzm*
     - :math:`\left<\varrho u_y u_z\right>`  :math:`\quad`  (mean Reynolds stress)
   * - *divrhourms*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm rms}`
   * - *divrhoumax*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm max}`
   * - *rlxm*
     - :math:`\left< \rho y u_z - z u_y \right>`
   * - *rlym*
     - :math:`\left< \rho z u_x - x u_z \right>`
   * - *rlzm*
     - :math:`\left< \rho x u_y - y u_x \right>`
   * - *rlx2m*
     - :math:`\left<(\rho y u_z-z u_y)^2\right>`
   * - *rly2m*
     - :math:`\left<(\rho z u_x-x u_z)^2\right>`
   * - *rlz2m*
     - :math:`\left<(\rho x u_y-y u_x)^2\right>`
   * - *tot_ang_mom*
     - Total angular momentum in spherical coordinates about the axis.
   * - *dtu*
     - :math:`\delta t/[c_{\delta t}\,\delta x /\max|\mathbf{u}|]`  :math:`\quad`  (time step relative to advective time step; see \S~\ref{time-step})
   * - *oum*
     - :math:`\left<\boldsymbol{\omega} \cdot\uv\right>`
   * - *ou_int*
     - :math:`\int_V\boldsymbol{\omega}\cdot\uv\,dV`
   * - *fum*
     - :math:`\left<\fv\cdot\uv\right>`
   * - *odel2um*
     - :math:`\left<\boldsymbol{\omega}\nabla^2\uv\right>`
   * - *o2m*
     - :math:`\left<\boldsymbol{\omega}^2\right> \equiv \left<(\curl\uv)^2\right>`
   * - *orms*
     - :math:`\left<\boldsymbol{\omega}^2\right>^{1/2}`
   * - *omax*
     - :math:`\max(|\boldsymbol{\omega}|)`
   * - *ox2m*
     - :math:`\left<\omega_x^2\right>`
   * - *oy2m*
     - :math:`\left<\omega_y^2\right>`
   * - *oz2m*
     - :math:`\left<\omega_z^2\right>`
   * - *oxm*
     -
   * - *oym*
     -
   * - *ozm*
     -
   * - *oxuzxm*
     - :math:`\left<\omega_x u_{z,x} \right>`
   * - *oyuzym*
     - :math:`\left<\omega_y u_{z,y} \right>`
   * - *oxoym*
     - :math:`\left<\omega_x\omega_y\right>`
   * - *oxozm*
     - :math:`\left<\omega_x\omega_z\right>`
   * - *oyozm*
     - :math:`\left<\omega_y\omega_z\right>`
   * - *qfm*
     - :math:`\left<\qv\cdot\fv\right>`
   * - *q2m*
     - :math:`\left<\qv^2\right>`
   * - *qrms*
     - :math:`\left<\qv^2\right>^{1/2}`
   * - *qmax*
     - :math:`\max(|\qv|)`
   * - *qom*
     - :math:`\left<\qv\cdot\omv\right>`
   * - *quxom*
     - :math:`\left<\qv\cdot(\uv\times\omv)\right>`
   * - *oumphi*
     - :math:`\left<\omv\cdot\uv\right>_\varphi`
   * - *ozmphi*
     -
   * - *ormr*
     -
   * - *opmr*
     -
   * - *ozmr*
     -
   * - *dudx*
     - :math:`\left<\frac{\delta \uv}{\delta x}\right>`
   * - *Marms*
     - :math:`\left<\uv^2/\cs^2\right>`  :math:`\quad`  (rms Mach number)
   * - *Mamax*
     - :math:`\max |\uv|/\cs`  :math:`\quad`  (maximum Mach number)
   * - *fintm*
     -
   * - *fextm*
     -
   * - *duxdzma*
     -
   * - *duydzma*
     -
   * - *EEK*
     - :math:`\left<\varrho\uv^2\right>/2`
   * - *ekin*
     - :math:`\left<{1\over2}\varrho\uv^2\right>`
   * - *ekintot*
     - :math:`\int_V{1\over2}\varrho\uv^2\, dV`
   * - *totangmom*
     -
   * - *uxglnrym*
     - :math:`\left<u_x\partial_y\ln\varrho\right>`
   * - *uyglnrxm*
     - :math:`\left<u_y\partial_x\ln\varrho\right>`
   * - *uzdivum*
     - :math:`\left<u_z\nabla\cdot\uv\right>`
   * - *uxuydivum*
     - :math:`\left<u_x u_y\nabla\cdot\uv\right>`
   * - *divuHrms*
     - :math:`(\nabla_{\mathrm H}\cdot\uv_{\mathrm H})^{\mathrm rms}`
   * - *uxxrms*
     - :math:`u_{x,x}^{\mathrm rms}`
   * - *uyyrms*
     - :math:`u_{y,y}^{\mathrm rms}`
   * - *uxzrms*
     - :math:`u_{x,z}^{\mathrm rms}`
   * - *uyzrms*
     - :math:`u_{y,z}^{\mathrm rms}`
   * - *uzyrms*
     - :math:`u_{z,y}^{\mathrm rms}`
   * - *udpxxm*
     - components of symmetric tensor :math:`\left< u_i \partial_j p + u_j \partial_i p \right>`

Module *interstellar.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *taucmin*
     - :math:`\min(\tau_{\mathrm cool})`
   * - *Hmax_ism*
     - :math:`\max(\Gamma-\rho\Lambda)`
   * - *Lamm*
     - :math:`\left<\Lambda\right>`
   * - *nrhom*
     - TBC
   * - *rhoLm*
     - :math:`\left<\rho\Lambda\right>`
   * - *Gamm*
     - :math:`\left<\Gamma\right>`

Module *lorenz_gauge.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *phim*
     - :math:`\left<\phi\right>`
   * - *phipt*
     - :math:`\phi(x1,y1,z1)`
   * - *phip2*
     - :math:`\phi(x2,y2,z2)`
   * - *phibzm*
     - :math:`\left<\phi B_z\right>`
   * - *phibzmz*
     - :math:`\left<\phi B_z\right>_{xy}`

Module *magnetic.f90*
~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *eta_tdep*
     - :math:`t` -dependent :math:`\eta`
   * - *ab_int*
     - :math:`\int\Av\cdot\Bv\;dV`
   * - *jb_int*
     - :math:`\int\jv\cdot\Bv\;dV`
   * - *b2tm*
     - :math:`\left<\bv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *bjtm*
     - :math:`\left<\bv(t)\cdot\int_0^t\jv(t') dt'\right>`
   * - *jbtm*
     - :math:`\left<\jv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *ujtm*
     - :math:`\left<\uv(t)\cdot\int_0^t\jv(t') dt'\right>`
   * - *jutm*
     - :math:`\left<\jv(t)\cdot\int_0^t\uv(t') dt'\right>`
   * - *ubtm*
     - :math:`\left<\uv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *butm*
     - :math:`\left<\bv(t)\cdot\int_0^t\uv(t') dt'\right>`
   * - *b2ruzm*
     - :math:`\left<\Bv^2\rho u_z\right>`
   * - *b2uzm*
     - :math:`\left<\Bv^2u_z\right>`
   * - *ubbzm*
     - :math:`\left<(\uv\cdot\Bv)B_z\right>`
   * - *b1m*
     - :math:`\left<|\Bv|\right>`
   * - *b2m*
     - :math:`\left<\Bv^2\right>`
   * - *EEM*
     - :math:`\left<\Bv^2\right>/2`
   * - *EEM2*
     - :math:`\left<(\Bv^2/2)^2\right>`
   * - *EEM3*
     - :math:`\left<(\Bv^2/2)^3\right>`
   * - *EEM4*
     - :math:`\left<(\Bv^2/2)^4\right>`
   * - *b4m*
     - :math:`\log_{10}\left<\Bv^4\right>`
   * - *b6m*
     - :math:`\log_{10}\left<\Bv^6\right>`
   * - *b8m*
     - :math:`\log_{10}\left<\Bv^8\right>`
   * - *b12m*
     - :math:`\log_{10}\left<\Bv^{12}\right>`
   * - *logbm*
     - :math:`\left<\log B\right>`
   * - *bm2*
     - :math:`\max(\Bv^2)`
   * - *j2m*
     - :math:`\left<\jv^2\right>`
   * - *jm2*
     - :math:`\max(\jv^2)`
   * - *abm*
     - :math:`\left<\Av\cdot\Bv\right>`
   * - *acbm*
     - :math:`\left<\Av_\mathrm{Cou}\cdot\Bv\right>`
   * - *gLamam*
     - :math:`\left<\nabla\Lambda\cdot\Av\right>`
   * - *gLambm*
     - :math:`\left<\nabla\Lambda\cdot\Bv\right>`
   * - *abumx*
     - :math:`\left<u_x\Av\cdot\Bv\right>`
   * - *abumy*
     - :math:`\left<u_y\Av\cdot\Bv\right>`
   * - *abumz*
     - :math:`\left<u_z\Av\cdot\Bv\right>`
   * - *abmh*
     - :math:`\left<\Av\cdot\Bv\right>`  (temp)
   * - *abmn*
     - :math:`\left<\Av\cdot\Bv\right>`  (north)
   * - *abms*
     - :math:`\left<\Av\cdot\Bv\right>`  (south)
   * - *abrms*
     - :math:`\left<(\Av\cdot\Bv)^2\right>^{1/2}`
   * - *jbrms*
     - :math:`\left<(\jv\cdot\Bv)^2\right>^{1/2}`
   * - *jxbrms*
     - :math:`\left<(\jv\times\Bv)^2\right>^{1/2}`
   * - *ajm*
     - :math:`\left<\jv\cdot\Av\right>`
   * - *jbm*
     - :math:`\left<\jv\cdot\Bv\right>`
   * - *a2b2m*
     - :math:`\left<\Av^2\cdot\Bv^2\right>`
   * - *j2b2m*
     - :math:`\left<\jv^2\cdot\Bv^2\right>`
   * - *hjbm*
     -
   * - *jbmh*
     - :math:`\left<\Jv\cdot\Bv\right>`  (temp)
   * - *jbmn*
     - :math:`\left<\Jv\cdot\Bv\right>`  (north)
   * - *jbms*
     - :math:`\left<\Jv\cdot\Bv\right>`  (south)
   * - *ubm*
     - :math:`\left<\uv\cdot\Bv\right>`
   * - *dubrms*
     - :math:`\left<(\uv-\Bv)^2\right>^{1/2}`
   * - *dobrms*
     - :math:`\left<(\boldsymbol{\omega}-\Bv)^2 \right>^{1/2}`
   * - *uxbxm*
     - :math:`\left<u_xB_x\right>`
   * - *uybxm*
     - :math:`\left<u_yB_x\right>`
   * - *uzbxm*
     - :math:`\left<u_zB_x\right>`
   * - *uxbym*
     - :math:`\left<u_xB_y\right>`
   * - *uybym*
     - :math:`\left<u_yB_y\right>`
   * - *uzbym*
     - :math:`\left<u_zB_y\right>`
   * - *uxbzm*
     - :math:`\left<u_xB_z\right>`
   * - *uybzm*
     - :math:`\left<u_yB_z\right>`
   * - *uzbzm*
     - :math:`\left<u_zB_z\right>`
   * - *uxjxm*
     - :math:`\left<u_xJ_x\right>`
   * - *uxjym*
     - :math:`\left<u_xJ_y\right>`
   * - *uxjzm*
     - :math:`\left<u_xJ_z\right>`
   * - *uyjxm*
     - :math:`\left<u_yJ_x\right>`
   * - *uyjym*
     - :math:`\left<u_yJ_y\right>`
   * - *uyjzm*
     - :math:`\left<u_yJ_z\right>`
   * - *uzjxm*
     - :math:`\left<u_zJ_x\right>`
   * - *uzjym*
     - :math:`\left<u_zJ_y\right>`
   * - *uzjzm*
     - :math:`\left<u_zJ_z\right>`
   * - *cosubm*
     - :math:`\left<\Uv\cdot\Bv/(|\Uv|\,|\Bv|) \right>`
   * - *jxbxm*
     - :math:`\left<j_xB_x\right>`
   * - *jybxm*
     - :math:`\left<j_yB_x\right>`
   * - *jzbxm*
     - :math:`\left<j_zB_x\right>`
   * - *jxbym*
     - :math:`\left<j_xB_y\right>`
   * - *jybym*
     - :math:`\left<j_yB_y\right>`
   * - *jzbym*
     - :math:`\left<j_zB_y\right>`
   * - *jxbzm*
     - :math:`\left<j_xB_z\right>`
   * - *jybzm*
     - :math:`\left<j_yB_z\right>`
   * - *jzbzm*
     - :math:`\left<j_zB_z\right>`
   * - *uam*
     - :math:`\left<\uv\cdot\Av\right>`
   * - *obm*
     - :math:`\left<\ov\cdot\Bv\right>`
   * - *ujm*
     - :math:`\left<\uv\cdot\Jv\right>`
   * - *fbm*
     - :math:`\left<\fv\cdot\Bv\right>`
   * - *fxbxm*
     - :math:`\left<f_x B_x\right>`
   * - *epsM*
     - :math:`\left<\eta\mu_0\jv^2\right>`
   * - *epsM2*
     - :math:`\left<(\eta\mu_0\jv^2)^2\right>`
   * - *epsM3*
     - :math:`\left<(\eta\mu_0\jv^2)^3\right>`
   * - *epsM4*
     - :math:`\left<(\eta\mu_0\jv^2)^4\right>`
   * - *epsAD*
     - :math:`\left<\rho^{-1} t_{\mathrm AD} (\vec{J}\times\vec{B})^2\right>`  (heating by ion-neutrals friction)
   * - *bxpt*
     - :math:`B_x(x_1,y_1,z_1,t)`
   * - *bypt*
     - :math:`B_y(x_1,y_1,z_1,t)`
   * - *bzpt*
     - :math:`B_z(x_1,y_1,z_1,t)`
   * - *bxbypt*
     - :math:`(B_x B_y)(x_1,y_1,z_1,t)`
   * - *bybzpt*
     - :math:`(B_y B_z)(x_1,y_1,z_1,t)`
   * - *bzbxpt*
     - :math:`(B_z B_x)(x_1,y_1,z_1,t)`
   * - *jxpt*
     - :math:`J_x(x_1,y_1,z_1,t)`
   * - *jypt*
     - :math:`J_y(x_1,y_1,z_1,t)`
   * - *jzpt*
     - :math:`J_z(x_1,y_1,z_1,t)`
   * - *Expt*
     - :math:`{\cal E}_x(x_1,y_1,z_1,t)`
   * - *Eypt*
     - :math:`{\cal E}_y(x_1,y_1,z_1,t)`
   * - *Ezpt*
     - :math:`{\cal E}_z(x_1,y_1,z_1,t)`
   * - *axpt*
     - :math:`A_x(x_1,y_1,z_1,t)`
   * - *aypt*
     - :math:`A_y(x_1,y_1,z_1,t)`
   * - *azpt*
     - :math:`A_z(x_1,y_1,z_1,t)`
   * - *bxp2*
     - :math:`B_x(x_2,y_2,z_2,t)`
   * - *byp2*
     - :math:`B_y(x_2,y_2,z_2,t)`
   * - *bzp2*
     - :math:`B_z(x_2,y_2,z_2,t)`
   * - *jxp2*
     - :math:`J_x(x_2,y_2,z_2,t)`
   * - *jyp2*
     - :math:`J_y(x_2,y_2,z_2,t)`
   * - *jzp2*
     - :math:`J_z(x_2,y_2,z_2,t)`
   * - *Exp2*
     - :math:`{\cal E}_x(x_2,y_2,z_2,t)`
   * - *Eyp2*
     - :math:`{\cal E}_y(x_2,y_2,z_2,t)`
   * - *Ezp2*
     - :math:`{\cal E}_z(x_2,y_2,z_2,t)`
   * - *axp2*
     - :math:`A_x(x_2,y_2,z_2,t)`
   * - *ayp2*
     - :math:`A_y(x_2,y_2,z_2,t)`
   * - *azp2*
     - :math:`A_z(x_2,y_2,z_2,t)`
   * - *epsM_LES*
     -
   * - *aybym2*
     -
   * - *exaym2*
     -
   * - *exabot*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm bot}`
   * - *exatop*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm top}`
   * - *exjm2*
     -
   * - *emag*
     - :math:`\int_V{1\over2\mu_0}\Bv^2\, dV`
   * - *km0EM*
     - :math:`\int E_M(k)\,dk`
   * - *km1EM*
     - :math:`\int k^{-1} E_M(k)\,dk`
   * - *brms*
     - :math:`\left<\Bv^2\right>^{1/2}`
   * - *bfrms*
     - :math:`\left<{\Bv'}^2\right>^{1/2}`
   * - *bf2m*
     - :math:`\left<{\Bv'}^2\right>`
   * - *bf4m*
     - :math:`\left<{\Bv'}^4\right>`
   * - *bmax*
     - :math:`\max(|\Bv|)`
   * - *bxmin*
     - :math:`\min(|B_x|)`
   * - *bymin*
     - :math:`\min(|B_y|)`
   * - *bzmin*
     - :math:`\min(|B_z|)`
   * - *bxmax*
     - :math:`\max(|B_x|)`
   * - *bymax*
     - :math:`\max(|B_y|)`
   * - *bzmax*
     - :math:`\max(|B_z|)`
   * - *bbxmax*
     - :math:`\max(|B_x|) excluding Bv_{ext}`
   * - *bbymax*
     - :math:`\max(|B_y|) excluding Bv_{ext}`
   * - *bbzmax*
     - :math:`\max(|B_z|) excluding Bv_{ext}`
   * - *jxmax*
     - :math:`\max(|jv_x|)`
   * - *jymax*
     - :math:`\max(|jv_y|)`
   * - *jzmax*
     - :math:`\max(|jv_z|)`
   * - *jrms*
     - :math:`\left<\jv^2\right>^{1/2}`
   * - *hjrms*
     - :math:`\left<\jv^2\right>^{1/2}`
   * - *jmax*
     - :math:`\max(|\jv|)`
   * - *vA23rms*
     - :math:`\left<\Bv^2/\varrho^{4/3}\right>^{1/2}`
   * - *vArms*
     - :math:`\left<\Bv^2/\varrho\right>^{1/2}`
   * - *vAmax*
     - :math:`\max(\Bv^2/\varrho)^{1/2}`
   * - *dtb*
     - :math:`\delta t / [c_{\delta t}\,\delta x /v_{\mathrm A,max}]`  :math:`\quad`  (time step relative to Alfv{\'e}n time step; see \S~\ref{time-step})
   * - *dteta*
     - :math:`\delta t/[c_{\delta t,{\mathrm v}}\, \delta x^2/\eta_{\mathrm max}]`  :math:`\quad`  (time step relative to resistive time step; see \S~\ref{time-step})
   * - *dteta3*
     - :math:`\delta t/[c_{\delta t,{\mathrm v3}}\, \delta x^6/\eta^{\mathrm hyper}_{\mathrm max}]`  :math:`\quad`  (time step relative to hyper resistive time step; see \S~\ref{time-step})
   * - *dtHr*
     -
   * - *dtFr*
     -
   * - *dtBr*
     -
   * - *axm*
     -
   * - *aym*
     -
   * - *azm*
     -
   * - *a2m*
     - :math:`\left<\Av^2\right>`
   * - *arms*
     - :math:`\left<\Av^2\right>^{1/2}`
   * - *az2m*
     - :math:`\left<A_z^2\right>`
   * - *amax*
     - :math:`\max(|\Av|)`
   * - *divarms*
     - :math:`\langle(\nabla\cdot\Av)^2\rangle^{1/2}`
   * - *beta1m*
     - :math:`\left<\Bv^2/(2\mu_0 p)\right>`  :math:`\quad`  (mean inverse plasma beta)
   * - *beta1max*
     - :math:`\max[\Bv^2/(2\mu_0 p)]`  :math:`\quad`  (maximum inverse plasma beta)
   * - *betam*
     - :math:`\langle\beta\rangle`
   * - *betamax*
     - :math:`\max\beta`
   * - *betamin*
     - :math:`\min\beta`
   * - *Azmid_min*
     - :math:`\min A_z^{\mathrm mid}(x)`
   * - *Azmid_max*
     - :math:`\max A_z^{\mathrm mid}(x)`
   * - *Azmid_ymin*
     - :math:`\min A_z^{\mathrm mid}(y)`
   * - *Azmid_ymax*
     - :math:`\max A_z^{\mathrm mid}(y)`
   * - *bxm*
     - :math:`\left<B_x\right>`
   * - *bym*
     - :math:`\left<B_y\right>`
   * - *bzm*
     - :math:`\left<B_z\right>`
   * - *jxm*
     - :math:`\left<J_x\right>`
   * - *jym*
     - :math:`\left<J_y\right>`
   * - *jzm*
     - :math:`\left<J_z\right>`
   * - *bxbym*
     - :math:`\left<B_x B_y\right>`
   * - *bxbzm*
     - :math:`\left<B_x B_z\right>`
   * - *bybzm*
     - :math:`\left<B_y B_z\right>`
   * - *djuidjbim*
     -
   * - *bij_cov_diffmax*
     - difference between two implementations of covariant derivatives
   * - *bmx*
     - :math:`\left<\left<\Bv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean field)
   * - *bmy*
     - :math:`\left<\left<\Bv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean field)
   * - *bmz*
     - :math:`\left<\left<\Bv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean field)
   * - *bmzS2*
     - :math:`\left<\left<\Bv_S\right>_{xy}^2\right>`
   * - *bmzA2*
     - :math:`\left<\left<\Bv_A\right>_{xy}^2\right>`
   * - *jmx*
     - :math:`\left<\left<\Jv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean current density)
   * - *jmy*
     - :math:`\left<\left<\Jv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean current density)
   * - *jmz*
     - :math:`\left<\left<\Jv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean current density)
   * - *bmzph*
     - Phase of a Beltrami field
   * - *bmzphe*
     - Error of phase of a Beltrami field
   * - *bsinphz*
     - sine of phase of a Beltrami field
   * - *bcosphz*
     - cosine of phase of a Beltrami field
   * - *emxamz3*
     - :math:`\left<\left<\Ev\right>_{xy}\times\left<\Av\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity flux)
   * - *embmz*
     - :math:`\left<\left<\Ev\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity production )
   * - *ambmz*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field)
   * - *ambmzh*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, temp)
   * - *ambmzn*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, north)
   * - *ambmzs*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, south)
   * - *jmbmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (current helicity of :math:`xy` -averaged mean field)
   * - *Rmmz*
     - :math:`\left<\frac{|\uv\times\Bv|}{|\eta\Jv|} \right>_{xy}`
   * - *kx_aa*
     - :math:`k_x`
   * - *kmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>/ \left<\left<\Bv\right>_{xy}^2\right>`
   * - *bx2m*
     - :math:`\left< B_x^2 \right>`
   * - *by2m*
     - :math:`\left< B_y^2 \right>`
   * - *bz2m*
     - :math:`\left< B_z^2 \right>`
   * - *bx3m*
     - :math:`\left< B_x^3 \right>`
   * - *by3m*
     - :math:`\left< B_y^3 \right>`
   * - *bz3m*
     - :math:`\left< B_z^3 \right>`
   * - *bx4m*
     - :math:`\left< B_x^4 \right>`
   * - *by4m*
     - :math:`\left< B_y^4 \right>`
   * - *bz4m*
     - :math:`\left< B_z^4 \right>`
   * - *jx2m*
     - :math:`\left< J_x^2 \right>`
   * - *jy2m*
     - :math:`\left< J_y^2 \right>`
   * - *jz2m*
     - :math:`\left< J_z^2 \right>`
   * - *jx4m*
     - :math:`\left< J_x^4 \right>`
   * - *jy4m*
     - :math:`\left< J_y^4 \right>`
   * - *jz4m*
     - :math:`\left< J_z^4 \right>`
   * - *jz3m*
     - :math:`\left< J_z^3 \right>`
   * - *jh2m1*
     - :math:`\left< J_\perp^2 \right>^{I}`
   * - *jx2m1*
     - :math:`\left< J_x^2 \right>^{I}`
   * - *jy2m1*
     - :math:`\left< J_y^2 \right>^{I}`
   * - *jx2m2*
     - :math:`\left< J_x^2 \right>^{II}`
   * - *jy2m2*
     - :math:`\left< J_y^2 \right>^{II}`
   * - *jx2m3*
     - :math:`\left< J_x^2 \right>^{III}`
   * - *jy2m3*
     - :math:`\left< J_y^2 \right>^{III}`
   * - *uxbm*
     - :math:`\left<\uv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *jxbm*
     - :math:`\left<\jv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *vmagfricmax*
     - :math:`\max(1/\nu_{\mathrm mag}|\jv\times\Bv/\Bv^2|)`
   * - *vmagfricrms*
     - :math:`\left<1/\nu_{\mathrm mag}|\jv\times\Bv/\Bv^2|^2\right>^{1/2}`
   * - *oxuxbm*
     -
   * - *jxbxbm*
     -
   * - *gpxbm*
     -
   * - *uxDxuxbm*
     -
   * - *b3b21m*
     - :math:`\left<B_3 B_{2,1} \right>`
   * - *b3b12m*
     - :math:`\left<B_3 B_{1,2} \right>`
   * - *b1b32m*
     - :math:`\left<B_1 B_{3,2} \right>`
   * - *b1b23m*
     - :math:`\left<B_1 B_{2,3} \right>`
   * - *b2b13m*
     - :math:`\left<B_2 B_{1,3} \right>`
   * - *b2b31m*
     - :math:`\left<B_2 B_{3,1} \right>`
   * - *udotxbm*
     -
   * - *uxbdotm*
     -
   * - *uxbmx*
     - :math:`\left<(\uv\times\Bv)_x\right>`
   * - *uxbmy*
     - :math:`\left<(\uv\times\Bv)_y\right>`
   * - *uxbmz*
     - :math:`\left<(\uv\times\Bv)_z\right>`
   * - *jxbmx*
     - :math:`\left<(\jv\times\Bv)_x\right>`
   * - *jxbmy*
     - :math:`\left<(\jv\times\Bv)_y\right>`
   * - *jxbmz*
     - :math:`\left<(\jv\times\Bv)_z\right>`
   * - *uxbcmx*
     -
   * - *uxbcmy*
     -
   * - *uxbsmx*
     -
   * - *uxbsmy*
     -
   * - *examx*
     - :math:`\left<\Ev\times\Av\right>|_x`
   * - *examy*
     - :math:`\left<\Ev\times\Av\right>|_y`
   * - *examz*
     - :math:`\left<\Ev\times\Av\right>|_z`
   * - *exatotalmx*
     - :math:`\left<\Ev\times\Av\right>|_x`
   * - *exatotalmy*
     - :math:`\left<\Ev\times\Av\right>|_y`
   * - *exatotalmz*
     - :math:`\left<\Ev\times\Av\right>|_z`
   * - *exjmx*
     - :math:`\left<\Ev\times\Jv\right>|_x`
   * - *exjmy*
     - :math:`\left<\Ev\times\Jv\right>|_y`
   * - *exjmz*
     - :math:`\left<\Ev\times\Jv\right>|_z`
   * - *dexbmx*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_x`
   * - *dexbmy*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_y`
   * - *dexbmz*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_z`
   * - *phibmx*
     - :math:`\left<\phi\Bv\right>|_x`
   * - *phibmy*
     - :math:`\left<\phi\Bv\right>|_y`
   * - *phibmz*
     - :math:`\left<\phi\Bv\right>|_z`
   * - *uxjm*
     -
   * - *b2divum*
     - :math:`\left<\Bv^2\nabla\cdot\uv\right>`
   * - *jdel2am*
     - :math:`\left<\Jv\cdot\nabla^2\Av)\right>`
   * - *jem*
     - :math:`\left<\jv\cdot\Ev\right>`
   * - *aem*
     - :math:`\left<\Av\cdot\Ev\right>`
   * - *ujxbm*
     - :math:`\left<\uv\cdot(\Jv\times\Bv)\right>`
   * - *WL2D*
     - :math:`\left<J_i u_j A_{i,j} \right>`
   * - *WL3D*
     - :math:`-\left<J_i u_j A_{j,i} \right>`
   * - *WL3D2*
     - :math:`\left<J_i A_j u_{j,i} \right>`
   * - *gb2m*
     - :math:`\left<|\hat{B}_{i,j}|^2\right>`
   * - *bij2m*
     - :math:`\left<|\hat{B}_{i,j}|^2\right>`
   * - *sijbibjm*
     - :math:`\left<S_{i,j} B_i B_j\right>`
   * - *ubgbpm*
     - :math:`\left<\uv\cdot(\Bv\cdot\nabla\Bv)\right>`
   * - *ugb22m*
     - :math:`\left<\uv\cdot\nabla\Bv^2/2)\right>`
   * - *jxbrxm*
     -
   * - *jxbrym*
     -
   * - *jxbrzm*
     -
   * - *jxbrmax*
     - :math:`\max(|\Jv\times\Bv/\rho|)`
   * - *jxbr2m*
     - :math:`\left<(\Jv\times\Bv/\rho)^2\right>`
   * - *jxbrqm*
     - :math:`\left<(\Jv\times\Bv/\rho)\cdot\mathbf{q}\right>`
   * - *uxBrms*
     -
   * - *Bresrms*
     -
   * - *Rmrms*
     -
   * - *jfm*
     -
   * - *brbpmr*
     -
   * - *vA2m*
     -
   * - *b2mr*
     -
   * - *brmr*
     -
   * - *bpmr*
     -
   * - *bzmr*
     -
   * - *armr*
     -
   * - *apmr*
     -
   * - *azmr*
     -
   * - *mflux_x*
     -
   * - *mflux_y*
     -
   * - *mflux_z*
     -
   * - *bmxy_rms*
     - :math:`\sqrt{[\left<b_x\right>_z(x,y)]^2 + [\left<b_y\right>_z(x,y)]^2 + [\left<b_z\right>_z(x,y)]^2}`
   * - *etasmagm*
     - Mean of Smagorinsky resistivity
   * - *etasmagmin*
     - Min of Smagorinsky resistivity
   * - *etasmagmax*
     - Max of Smagorinsky resistivity
   * - *etavamax*
     - Max of artificial resistivity :math:`\eta\sim v_A`
   * - *etajmax*
     - Max of artificial resistivity :math:`\eta\sim J / \sqrt{\rho}`
   * - *etaj2max*
     - Max of artificial resistivity :math:`\eta\sim J^2 / \rho`
   * - *etajrhomax*
     - Max of artificial resistivity :math:`\eta\sim J / \rho`
   * - *etaaniso*
     - :math:`\eta_1`
   * - *etaanisoBB*
     - :math:`\eta_{BB}`
   * - *cosjbm*
     - :math:`\left<\Jv\cdot\Bv/(|\Jv|\,|\Bv|)\right>`
   * - *coshjbm*
     -
   * - *jparallelm*
     - Mean value of the component of J parallel to B
   * - *jperpm*
     - Mean value of the component of J perpendicular to B
   * - *hjparallelm*
     - Mean value of the component of :math:`J_{\mathrm hyper}`  parallel to B
   * - *hjperpm*
     - Mean value of the component of :math:`J_{\mathrm hyper}`  perpendicular to B
   * - *b2sphm*
     - :math:`\int_{r=0}^{r=r_{\mathrm diag}} \Bv^2 dV` , where :math:`r=\sqrt{x^2+y^2+z^2}`
   * - *brmsx*
     - :math:`\left<\Bv^2\right>^{1/2}`  for the magnetic_xaver_range
   * - *brmsz*
     - :math:`\left<\Bv^2\right>^{1/2}`  for the magnetic_zaver_range
   * - *Exmxy*
     - :math:`\left<{\cal E}_x\right>_{z}`
   * - *Eymxy*
     - :math:`\left<{\cal E}_y\right>_{z}`
   * - *Ezmxy*
     - :math:`\left<{\cal E}_z\right>_{z}`

Module *neutralvelocity.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *epsKn*
     - :math:`\left<2\nu_n\varrho_n\Strain_n^2\right>`

Module *noentropy.f90*
~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtc*
     - :math:`\delta t/[c_{\delta t}\,\delta_x /\max c_{\mathrm s}]`  :math:`\quad`  (time step relative to acoustic time step; see \S~\ref{time-step})
   * - *ethm*
     - :math:`\left<\varrho e\right>`  :math:`\quad`  (mean thermal [=internal] energy)
   * - *pdivum*
     - :math:`\left<p\nabla\uv\right>`
   * - *csm*
     - :math:`\left<c_{\mathrm s}\right>`

Module *particles_caustics.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *TrSigmapm*
     - :math:`\langle{\mathrm Tr}\left[\sigma\right]\rangle`
   * - *blowupm*
     - Mean no. of times :math:`\sigma`  falls below cutoff
   * - *lnVpm*
     - Mean of (logarithm of) Volume around an inertial particle

Module *particles_chemistry.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Shchm*
     - :math:`\mathrm{mean particle Sherwood number}`

Module *particles_dust.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *xpm*
     - :math:`x_{part}`
   * - *xpmin*
     - :math:`x_{part}`
   * - *xpmax*
     - :math:`x_{part}`
   * - *xp2m*
     - :math:`x^2_{part}`
   * - *vrelpabsm*
     - :math:`\mathrm{Absolute value of mean relative velocity}`
   * - *vpxm*
     - :math:`u_{part}`
   * - *vpx2m*
     - :math:`u^2_{part}`
   * - *vpcoalx2m*
     - :math:`v^2_{coal,part}`  (thermal/coalescence velocity entering the Brownian coagulation kernel)
   * - *ekinp*
     - :math:`E_{kin,part}`
   * - *vpxmax*
     - :math:`MAX(u_{part})`
   * - *vpxmin*
     - :math:`MIN(u_{part})`
   * - *npm*
     - :math:`\mathrm{mean particle number density}`

Module *particles_lagrangian.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *xpm*
     - :math:`x_{part}`
   * - *xp2m*
     - :math:`x^2_{part}`
   * - *vrelpabsm*
     - :math:`\mathrm{Absolute value of mean relative velocity}`
   * - *vpxm*
     - :math:`u_{part}`
   * - *vpx2m*
     - :math:`u^2_{part}`
   * - *ekinp*
     - :math:`E_{kin,part}`
   * - *vpxmax*
     - :math:`MAX(u_{part})`
   * - *vpxmin*
     - :math:`MIN(u_{part})`
   * - *npm*
     - :math:`\mathrm{mean particle number density}`

Module *particles_mass_swarm.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *mpm*
     - :math:`\overline{m_p}`
   * - *mpmin*
     - :math:`\min_j m_{p,j}`
   * - *mpmax*
     - :math:`\max_j m_{p,j}`

Module *particles_surfspec.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtpchem*
     - :math:`dt_{particle,chemistry}`

Module *polymer.f90*
~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *polytrm*
     - :math:`\left\langle Tr[C_{ij}]\right\rangle`
   * - *frmax*
     - :math:`\max(f(r))`

Module *pscalar.f90*
~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rhoccm*
     - :math:`\left<\varrho c\right>`
   * - *ccmax*
     - :math:`\max(c)`
   * - *ccmin*
     -
   * - *lnccm*
     -
   * - *mcct*
     -
   * - *gcc5m*
     -
   * - *gcc10m*
     -
   * - *ucm*
     -
   * - *uudcm*
     -
   * - *Cz2m*
     -
   * - *Cz4m*
     -
   * - *cc1m*
     -
   * - *cc2m*
     -
   * - *cc3m*
     -
   * - *cc4m*
     -
   * - *cc5m*
     -
   * - *cc6m*
     -
   * - *cc7m*
     -
   * - *cc8m*
     -
   * - *cc9m*
     -
   * - *cc10m*
     -
   * - *Crmsm*
     -
   * - *gcc1m*
     -
   * - *gcc2m*
     -
   * - *gcc3m*
     -
   * - *gcc4m*
     -
   * - *gcc6m*
     -
   * - *gcc7m*
     -
   * - *gcc8m*
     -
   * - *gcc9m*
     -
   * - *lnccmz*
     -
   * - *lnccmy*
     -
   * - *lnccmx*
     -
   * - *ccmz*
     -
   * - *ccmy*
     -
   * - *ccmx*
     -
   * - *ccglnrm*
     - :math:`\left<c\nabla_z\varrho\right>`

Module *pscalar_nolog.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rhoccm*
     - :math:`\left< \rho c_1 \right>`
   * - *ccmax*
     - :math:`\max[ c_1 ]`
   * - *ccmin*
     - :math:`\min[ c_1 ]`
   * - *mrclncm*
     - :math:`\left< \rho c_1 \log(c_1) \right>`
   * - *rhoccmax*
     - :math:`\max[ \rho c_1 ]`
   * - *rhoc2m*
     - :math:`\left< \rho c_2 \right>`
   * - *rhoc3m*
     - :math:`\left< \rho c_3 \right>`
   * - *ucm*
     - :math:`\left< u_z c_1 \right>`  (with an extra factor of :math:`2 \cos(z)`  if lgradC_profile=T)
   * - *uudcm*
     - :math:`\left< u_z \vec{u} \cdot \grad c_1 \right>`
   * - *Cz2m*
     - :math:`\left< \rho c_1 z^2 \right>`
   * - *Cz4m*
     - :math:`\left< \rho c_1 z^4 \right>`
   * - *Crmsm*
     - :math:`\sqrt{\left< \rho c_1^2 \right>}`
   * - *ccrms*
     - :math:`\sqrt{\left< c_1^2 \right>}`
   * - *uxcm*
     - :math:`\left< u_x c_1 \right>`
   * - *uycm*
     - :math:`\left< u_y c_1 \right>`
   * - *uzcm*
     - :math:`\left< u_z c_1 \right>`
   * - *cc1m*
     - :math:`\left< c_1^{-1} \right>`
   * - *cc2m*
     - :math:`\left< c_1^{-2} \right>`
   * - *cc3m*
     - :math:`\left< c_1^{-3} \right>`
   * - *cc4m*
     - :math:`\left< c_1^{-4} \right>`
   * - *cc5m*
     - :math:`\left< c_1^{-5} \right>`
   * - *cc6m*
     - :math:`\left< c_1^{-6} \right>`
   * - *cc7m*
     - :math:`\left< c_1^{-7} \right>`
   * - *cc8m*
     - :math:`\left< c_1^{-8} \right>`
   * - *cc9m*
     - :math:`\left< c_1^{-9} \right>`
   * - *cc10m*
     - :math:`\left< c_1^{-10} \right>`
   * - *gcc1m*
     - :math:`\left< \left|\grad c_1 \right| \right>`
   * - *gcc2m*
     - :math:`\left< \left|\grad c_1 \right|^{2} \right>`
   * - *gcc3m*
     - :math:`\left< \left|\grad c_1 \right|^{3} \right>`
   * - *gcc4m*
     - :math:`\left< \left|\grad c_1 \right|^{4} \right>`
   * - *gcc5m*
     - :math:`\left< \left|\grad c_1 \right|^{5} \right>`
   * - *gcc6m*
     - :math:`\left< \left|\grad c_1 \right|^{6} \right>`
   * - *gcc7m*
     - :math:`\left< \left|\grad c_1 \right|^{7} \right>`
   * - *gcc8m*
     - :math:`\left< \left|\grad c_1 \right|^{8} \right>`
   * - *gcc9m*
     - :math:`\left< \left|\grad c_1 \right|^{9} \right>`
   * - *gcc10m*
     - :math:`\left< \left|\grad c_1 \right|^{10} \right>`
   * - *cugccm*
     - :math:`\left< c_1 \vec{u} \cdot \grad c_1 \right>`
   * - *ccugum*
     - :math:`\left< c_1 \vec{u} \cdot \grad u_z \right>`
   * - *ccglnrm*
     - :math:`\left< c_1 \partial_z \log\rho \right>`

Module *selfgravity.f90*
~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rugpotselfm*
     - :math:`\left<\rho\uv\cdot\nabla\Phi\right>`
   * - *gpotself2m*
     - :math:`\left<(\nabla\Phi)^2\right>`

Module *shear.f90*
~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtshear*
     - advec\_shear/cdt
   * - *deltay*
     - deltay

Module *shock.f90*
~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *shockmax*
     - Max shock number

Module *shock_highorder.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *shockm*
     -
   * - *shockmin*
     -
   * - *shockmax*
     -
   * - *gshockmax*
     - :math:`\max\left|\nabla\nu_{shock}\right|`

Module *solid_cells_ogrid_chemistry.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Ym*
     - :math:`\left<Y_x\right>`
   * - *dYm*
     - :math:`\delta\left<Y_x\right>/\delta t`
   * - *dYmax*
     - :math:`max\delta\left<Y_x\right>/\delta t`
   * - *Ymax*
     - :math:`\left<Y_{x,max}\right>`
   * - *Ymin*
     - :math:`\left<Y_{x,min}\right>`
   * - *hm*
     - :math:`\left<H_{x,max}\right>`
   * - *cpm*
     - :math:`\left<c_{p,x}\right>`
   * - *diffm*
     - :math:`\left<D_{x}\right>`
   * - *Ymz*
     - :math:`\left<Y_x\right>_{xy}(z)`
   * - *dtchem*
     - :math:`dt_{chem}`

Module *temperature_idealgas.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *TTmax*
     - :math:`\max (T)`
   * - *gTmax*
     - :math:`\max (|\nabla T|)`
   * - *TTmin*
     - :math:`\min (T)`
   * - *TTm*
     - :math:`\left< T \right>`
   * - *TTzmask*
     - :math:`\left< T \right>`  for the temp_zaver_range
   * - *TT2m*
     - :math:`\left< T^2 \right>`
   * - *TugTm*
     - :math:`\left< T\uv\cdot\nabla T \right>`
   * - *Trms*
     - :math:`\sqrt{\left< T^2 \right>}`
   * - *uxTm*
     - :math:`\left< u_x T \right>`
   * - *uyTm*
     - :math:`\left< u_y T \right>`
   * - *uzTm*
     - :math:`\left< u_z T \right>`
   * - *gT2m*
     - :math:`\left< (\nabla T)^2 \right>`
   * - *guxgTm*
     - :math:`\left< \nabla u_x \cdot \nabla T \right>`
   * - *guygTm*
     - :math:`\left< \nabla u_y \cdot \nabla T \right>`
   * - *guzgTm*
     - :math:`\left< \nabla u_z \cdot \nabla T \right>`
   * - *Tugux_uxugTm*
     - :math:`\left< T \uv\cdot\nabla u_x + u_x \uv\cdot\nabla T \right> =\left< \uv\cdot\nabla(u_x T) \right>`
   * - *Tuguy_uyugTm*
     - :math:`\left< T \uv\cdot\nabla u_y + u_y \uv\cdot\nabla T \right> =\left< \uv\cdot\nabla(u_y T) \right>`
   * - *Tuguz_uzugTm*
     - :math:`\left< T \uv\cdot\nabla u_z + u_z \uv\cdot\nabla T \right> =\left< \uv\cdot\nabla(u_z T) \right>`
   * - *Tdxpm*
     - :math:`\left< T dp/dx \right>`
   * - *Tdypm*
     - :math:`\left< T dp/dy \right>`
   * - *Tdzpm*
     - :math:`\left< T dp/dz \right>`
   * - *fradtop*
     - :math:`<-K{dT\over dz}>_{\text{top}}`  :math:`\quad`  (top radiative flux)
   * - *fradbot*
     - :math:`<-K{dT\over dz}>_{\text{bot}}`  :math:`\quad`  (bottom radiative flux)
   * - *yHmax*
     - DOCUMENT ME
   * - *yHmin*
     - DOCUMENT ME
   * - *yHm*
     - DOCUMENT ME
   * - *ethm*
     - :math:`\left< e_{\text{th}}\right> = \left< c_v \rho T \right> `  :math:`\quad`  (mean thermal energy)
   * - *eem*
     - :math:`\left< e \right> = \left< c_v T \right>`  :math:`\quad`  (mean internal energy)
   * - *ethtot*
     - :math:`\int_V\varrho e\,dV`  :math:`\quad`  (total thermal energy)
   * - *ssm*
     - :math:`\overline{S}`
   * - *thcool*
     - :math:`\tau_{\mathrm cool}`
   * - *ppm*
     - :math:`\overline{P}`
   * - *csm*
     - :math:`\overline{c}_{\mathrm s}`
   * - *csmax*
     - :math:`\max (c_{\mathrm s})`
   * - *dtc*
     - :math:`\delta t/[c_{\delta t}\,\delta_x /\max c_{\mathrm s}]`  :math:`\quad`  (time step relative to acoustic time step; see \S~\ref{time-step})
   * - *dtchi*
     - :math:`\delta t / [c_{\delta t,{\mathrm v}}\, \delta x^2/\chi_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat conductivity; see \S~\ref{time-step})

Module *temperature_ionization.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *TTmax*
     - :math:`\max (T)`
   * - *TTmin*
     - :math:`\min (T)`
   * - *TTm*
     - :math:`\left< T \right>`
   * - *yHmax*
     -
   * - *yHmin*
     -
   * - *yHm*
     -
   * - *ethm*
     - :math:`\left< e_{\text{th}}\right> = \left< c_v \rho T \right> `  :math:`\quad`  (mean thermal energy)
   * - *ssm*
     -
   * - *dtc*
     -
   * - *eem*
     - :math:`\left< e \right> `  :math:`\quad`  (mean internal energy)
   * - *ppm*
     - :math:`\left< p \right> `
   * - *Tppm*
     - :math:`\left<\max(p_{\mathrm thresh}-p,0)_{\mathrm norm}\right> `
   * - *mum*
     -
   * - *heatThm*
     - :math:`\alpha_\mathrm{Th}`
   * - *TTref*
     - :math:`T_\mathrm{ref}`
   * - *ppmax*
     -
   * - *ppmin*
     -

Module *testfield_axisym.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alpPERP*
     - :math:`\alpha_\perp`
   * - *alpPARA*
     - :math:`\alpha_\perp`
   * - *gam*
     - :math:`\gamma`
   * - *betPERP*
     - :math:`\beta_\perp`
   * - *betPARA*
     - :math:`\beta_\perp`
   * - *del*
     - :math:`\delta`
   * - *kapPERP*
     - :math:`\kappa_\perp`
   * - *kapPARA*
     - :math:`\kappa_\perp`
   * - *mu*
     - :math:`\mu`
   * - *alpPERPz*
     - :math:`\alpha_\perp(z)`
   * - *alpPARAz*
     - :math:`\alpha_\perp(z)`
   * - *gamz*
     - :math:`\gamma(z)`
   * - *betPERPz*
     - :math:`\beta_\perp(z)`
   * - *betPARAz*
     - :math:`\beta_\perp(z)`
   * - *delz*
     - :math:`\delta(z)`
   * - *kapPERPz*
     - :math:`\kappa_\perp(z)`
   * - *kapPARAz*
     - :math:`\kappa_\perp(z)`
   * - *muz*
     - :math:`\mu(z)`
   * - *bx1pt*
     - :math:`b_x^{1}`
   * - *bx2pt*
     - :math:`b_x^{2}`
   * - *bx3pt*
     - :math:`b_x^{3}`
   * - *b1rms*
     - :math:`\left<b_{1}^2\right>^{1/2}`
   * - *b2rms*
     - :math:`\left<b_{2}^2\right>^{1/2}`
   * - *b3rms*
     - :math:`\left<b_{3}^2\right>^{1/2}`

Module *testfield_axisym2.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alpPERP*
     - :math:`\alpha_\perp`
   * - *alpPARA*
     - :math:`\alpha_\perp`
   * - *gam*
     - :math:`\gamma`
   * - *betPERP*
     - :math:`\beta_\perp`
   * - *betPARA*
     - :math:`\beta_\perp`
   * - *del*
     - :math:`\delta`
   * - *kapPERP*
     - :math:`\kappa_\perp`
   * - *kapPARA*
     - :math:`\kappa_\perp`
   * - *mu*
     - :math:`\mu`
   * - *bx1pt*
     - :math:`b_x^{1}`
   * - *bx2pt*
     - :math:`b_x^{2}`
   * - *bx3pt*
     - :math:`b_x^{3}`
   * - *b1rms*
     - :math:`\left<b_{1}^2\right>^{1/2}`
   * - *b2rms*
     - :math:`\left<b_{2}^2\right>^{1/2}`
   * - *b3rms*
     - :math:`\left<b_{3}^2\right>^{1/2}`

Module *testfield_meri.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *E11xy*
     - :math:`E_{11xy}`
   * - *E12xy*
     - :math:`E_{12xy}`
   * - *E13xy*
     - :math:`E_{13xy}`
   * - *E21xy*
     - :math:`E_{21xy}`
   * - *E22xy*
     - :math:`E_{22xy}`
   * - *E23xy*
     - :math:`E_{23xy}`
   * - *E31xy*
     - :math:`E_{31xy}`
   * - *E32xy*
     - :math:`E_{32xy}`
   * - *E33xy*
     - :math:`E_{33xy}`
   * - *E41xy*
     - :math:`E_{41xy}`
   * - *E42xy*
     - :math:`E_{42xy}`
   * - *E43xy*
     - :math:`E_{43xy}`
   * - *E51xy*
     - :math:`E_{51xy}`
   * - *E52xy*
     - :math:`E_{52xy}`
   * - *E53xy*
     - :math:`E_{53xy}`
   * - *E61xy*
     - :math:`E_{61xy}`
   * - *E62xy*
     - :math:`E_{62xy}`
   * - *E63xy*
     - :math:`E_{63xy}`
   * - *E71xy*
     - :math:`E_{71xy}`
   * - *E72xy*
     - :math:`E_{72xy}`
   * - *E73xy*
     - :math:`E_{73xy}`
   * - *E81xy*
     - :math:`E_{81}`
   * - *E82xy*
     - :math:`E_{82}`
   * - *E83xy*
     - :math:`E_{83}`
   * - *E91xy*
     - :math:`E_{91}`
   * - *E92xy*
     - :math:`E_{92}`
   * - *E93xy*
     - :math:`E_{93}`
   * - *a11xy*
     - :math:`\alpha_{11}`
   * - *a12xy*
     - :math:`\alpha_{12}`
   * - *a13xy*
     - :math:`\alpha_{13}`
   * - *a21xy*
     - :math:`\alpha_{21}`
   * - *a22xy*
     - :math:`\alpha_{22}`
   * - *a23xy*
     - :math:`\alpha_{23}`
   * - *a31xy*
     - :math:`\alpha_{31}`
   * - *a32xy*
     - :math:`\alpha_{32}`
   * - *a33xy*
     - :math:`\alpha_{33}`
   * - *b111xy*
     - :math:`\b_{111}`
   * - *b121xy*
     - :math:`\b_{121}`
   * - *b131xy*
     - :math:`\b_{131}`
   * - *b211xy*
     - :math:`\b_{211}`
   * - *b221xy*
     - :math:`\b_{221}`
   * - *b231xy*
     - :math:`\b_{231}`
   * - *b311xy*
     - :math:`\b_{311}`
   * - *b321xy*
     - :math:`\b_{321}`
   * - *b331xy*
     - :math:`\b_{331}`
   * - *b112xy*
     - :math:`\b_{112}`
   * - *b122xy*
     - :math:`\b_{122}`
   * - *b132xy*
     - :math:`\b_{132}`
   * - *b212xy*
     - :math:`\b_{212}`
   * - *b222xy*
     - :math:`\b_{222}`
   * - *b232xy*
     - :math:`\b_{232}`
   * - *b312xy*
     - :math:`\b_{312}`
   * - *b322xy*
     - :math:`\b_{322}`
   * - *b332xy*
     - :math:`\b_{332}`

Module *testflow_z.f90*
~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *gal*
     - GAL-coefficients,     couple  :math:`\overline F`  and :math:`\overline U`
   * - *aklam*
     - AKA-:math:`\lambda` -tensor, couples :math:`\overline F`  and :math:`\overline W = \nabla\times{\overline U}`
   * - *gamma*
     - :math:`\gamma` -vector,      couples :math:`\overline F`  and :math:`\nabla\cdot{\overline U}`
   * - *nu*
     - :math:`\nu` -tensor,         couples :math:`\overline F`  and :math:`\partial^2 {\overline U}/\partial z^2`
   * - *zeta*
     - :math:`\zeta` -vector,       couples :math:`\overline F`  and :math:`{\overline G}_z = \nabla_z {\overline H}`
   * - *xi*
     - :math:`\xi` -vector,         couples :math:`\overline F`  and :math:`\partial^2 {\overline H}/\partial z^2`
   * - *aklamQ*
     - :math:`aklam^Q` -vector,     couples :math:`\overline Q`  and :math:`\overline W`
   * - *gammaQ*
     - :math:`\gamma^Q` -scalar,    couples :math:`\overline Q`  and :math:`\nabla\cdot{\overline U}=dU_z/dz`
   * - *nuQ*
     - :math:`\nu^Q` -vector,       couples :math:`\overline Q`  and :math:`\partial^2 \overline U/\partial z^2`
   * - *zetaQ*
     - :math:`\zeta^Q` -scalar,      couples :math:`\overline Q`  and :math:`{\overline G}_z`
   * - *xiQ*
     - :math:`\xi^Q` -scalar,        couples :math:`\overline Q`  and :math:`\partial^2 {\overline H}/\partial z^2`
   * - *galij*
     -
   * - *aklamij*
     - :math:`\alpha_{K,ij}`
   * - *gammai*
     - :math:`\gamma_i`
   * - *nuij*
     - :math:`\nu_{ij}`
   * - *zetai*
     - :math:`\zeta_i`  :math:`\xi_i`  :math:`\nu^Q_i`
   * - *aklamQi*
     - :math:`aklam^Q_i`
   * - *Fipq*
     - :math:`{\cal F}_i^{pq}`
   * - *Qpq*
     - :math:`{\cal Q}^{pq}`
   * - *upqrms*
     - :math:`\left<{u^{pq}}^2\right>`
   * - *hpqrms*
     - :math:`\left<{h^{pq}}^2\right>`
   * - *ux0mz*
     - :math:`\left<u_{x}\right>_{xy}`
   * - *uy0mz*
     - :math:`\left<u_{y}\right>_{xy}`
   * - *uz0mz*
     - :math:`\left<u_{z}\right>_{xy}`

Module *thermal_energy.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *TTmax*
     - :math:`\max (T)`
   * - *TTmin*
     - :math:`\min (T)`
   * - *ppm*
     - :math:`\left< p \right>`
   * - *TTm*
     - :math:`\left<T\right>`
   * - *ethm*
     - :math:`\left< e_{\text{th}}\right> = \left< c_v \rho T \right> `  :math:`\quad`  (mean thermal energy)
   * - *ethtot*
     - :math:`\int_V e_{\text{th}}\,dV`  :math:`\quad`  (total thermal energy)
   * - *ethmin*
     - :math:`\mathrm{min} e_\text{th}`
   * - *ethmax*
     - :math:`\mathrm{max} e_\text{th}`
   * - *eem*
     - :math:`\left< e \right> = \left< c_v T \right>`  :math:`\quad`  (mean internal energy)
   * - *etot*
     - :math:`\langle e_\textrm{th} + \rho u^2 / 2\rangle`

Module *training_torchfort.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *loss*
     - torchfort training loss
   * - *tauerror*
     - :math:`\sqrt{\left<(\sum_{i,j} u_i*u_j - tau_{ij})^2\right>}`

Module *viscosity.f90*
~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *nu_tdep*
     - time-dependent viscosity
   * - *ell_gam*
     - time-dependent photon mean free path
   * - *fviscm*
     - Mean value of viscous acceleration
   * - *fviscmin*
     - Min value of viscous acceleration (redundant)
   * - *fviscmax*
     - Max absolute viscous acceleration
   * - *fviscrmsx*
     - Rms value of viscous acceleration for the vis_xaver_range
   * - *num*
     - Mean value of viscosity
   * - *numax*
     - Max value of viscosity
   * - *numin*
     - Min value of viscosity
   * - *nusmagm*
     - Mean value of Smagorinsky viscosity
   * - *nusmagmin*
     - Min value of Smagorinsky viscosity
   * - *nusmagmax*
     - Max value of Smagorinsky viscosity
   * - *nu_LES*
     - Mean value of Smagorinsky viscosity
   * - *visc_heatm*
     - Mean value of viscous heating
   * - *qfviscm*
     - :math:`\left<\qv\cdot \fv_{\mathrm visc}\right>`
   * - *ufviscm*
     - :math:`\left<\uv\cdot \fv_{\mathrm visc}\right>`
   * - *Sij2m*
     - :math:`\left<\Strain^2\right>`
   * - *epsK*
     - :math:`\left<2\nu\varrho\Strain^2\right>`
   * - *epsK2*
     - :math:`\left<(2\nu\varrho\Strain^2)^2\right>`
   * - *epsK3*
     - :math:`\left<(2\nu\varrho\Strain^2)^3\right>`
   * - *epsK4*
     - :math:`\left<(2\nu\varrho\Strain^2)^4\right>`
   * - *epsKint*
     - :math:`\int(2\nu\varrho\Strain^2)\,dV`
   * - *epsK_LES*
     -
   * - *sijoiojm*
     - :math:`\left<S_{i,j} \omega_i \omega_j\right>`
   * - *dtnu*
     - :math:`\delta t/[c_{\delta t,{\mathrm v}}\, \delta x^2/\nu_{\mathrm max}]`  :math:`\quad`  (time step relative to viscous time step; see \S~\ref{time-step})
   * - *meshRemax*
     - Max mesh Reynolds number
   * - *mesh3Remax*
     - Max hyper3 mesh Reynolds number
   * - *Reshock*
     - Mesh Reynolds number at shock
   * - *nuD2uxbxm*
     -
   * - *nuD2uxbym*
     -
   * - *nuD2uxbzm*
     -

Module *experimental/anelastic.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rhom*
     - :math:`\left<\varrho\right>`  :math:`\quad`  (mean density)
   * - *rho2m*
     -
   * - *lnrho2m*
     -
   * - *drho2m*
     -
   * - *drhom*
     -
   * - *rhomin*
     -
   * - *rhomax*
     -
   * - *ugrhom*
     - :math:`\left<\uv\cdot\nabla\varrho\right>`
   * - *uglnrhom*
     -
   * - *dtd*
     -
   * - *rhomz*
     -
   * - *rhomy*
     -
   * - *rhomx*
     -
   * - *rhomxy*
     -
   * - *rhomxz*
     -
   * - *rhomr*
     -
   * - *totmass*
     -
   * - *mass*
     - :math:`\int\varrho\,dV`
   * - *divrhoum*
     - :math:`\left<\nabla\cdot(\varrho\uv)\right>`
   * - *divrhourms*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm rms}`
   * - *divrhoumax*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm max}`

Module *experimental/cosmicray_current.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *ekincr*
     - :math:`\left<{1\over2}\varrho\uv_{\mathrm cr}^2\right>`
   * - *ethmcr*
     - :math:`\left<\varrho_{\mathrm cr} e_{\mathrm cr}\right>`

Module *experimental/entropy_anelastic.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtc*
     - :math:`\delta t/[c_{\delta t}\,\delta_x /\max c_{\mathrm s}]`  :math:`\quad`  (time step relative to acoustic time step; see \S~\ref{time-step})
   * - *ethm*
     - :math:`\left<\varrho e\right>`  :math:`\quad`  (mean thermal [=internal] energy)
   * - *ethdivum*
     -
   * - *ssm*
     - :math:`\left<s/c_p\right>`  :math:`\quad`  (mean entropy)
   * - *ss2m*
     - :math:`\left<(s/c_p)^2\right>`  :math:`\quad`  (mean squared entropy)
   * - *eem*
     - :math:`\left<e\right>`
   * - *ppm*
     - :math:`\left<p\right>`
   * - *csm*
     - :math:`\left<c_{\mathrm s}\right>`
   * - *pdivum*
     - :math:`\left<p\nabla\uv\right>`
   * - *heatm*
     -
   * - *ugradpm*
     -
   * - *fradbot*
     - :math:`\int F_{\mathrm bot}\cdot d\vec{S}`
   * - *fradtop*
     - :math:`\int F_{\mathrm top}\cdot d\vec{S}`
   * - *TTtop*
     - :math:`\int T_{\mathrm top} d\vec{S}`
   * - *ethtot*
     - :math:`\int_V\varrho e\,dV`  :math:`\quad`  (total thermal [=internal] energy)
   * - *dtchi*
     - :math:`\delta t / [c_{\delta t,{\mathrm v}}\, \delta x^2/\chi_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat conductivity; see \S~\ref{time-step})
   * - *yHm*
     -
   * - *yHmax*
     -
   * - *TTm*
     -
   * - *TTmax*
     -
   * - *TTmin*
     -
   * - *fconvm*
     -
   * - *fconvz*
     -
   * - *dcoolz*
     -
   * - *fradz*
     -
   * - *fturbz*
     -
   * - *ssmx*
     -
   * - *ssmy*
     -
   * - *ssmz*
     -
   * - *TTp*
     -
   * - *ssmr*
     -
   * - *TTmx*
     -
   * - *TTmy*
     -
   * - *TTmz*
     -
   * - *TTmxy*
     -
   * - *TTmxz*
     -
   * - *TTmr*
     -
   * - *uxTTmz*
     -
   * - *uyTTmz*
     -
   * - *uzTTmz*
     -
   * - *uxTTmxy*
     -
   * - *ssmxy*
     - :math:`\left< s \right>_{z}`
   * - *ssmxz*
     - :math:`\left< s \right>_{y}`

Module *experimental/magnetic_shearboxJ.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *ab_int*
     - :math:`\int\Av\cdot\Bv\;dV`
   * - *jb_int*
     - :math:`\int\jv\cdot\Bv\;dV`
   * - *b2tm*
     - :math:`\left<\bv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *bjtm*
     - :math:`\left<\bv(t)\cdot\int_0^t\jv(t') dt'\right>`
   * - *jbtm*
     - :math:`\left<\jv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *b2ruzm*
     - :math:`\left<\Bv^2\rho u_z\right>`
   * - *b2uzm*
     - :math:`\left<\Bv^2u_z\right>`
   * - *ubbzm*
     - :math:`\left<(\uv\cdot\Bv)B_z\right>`
   * - *b1m*
     - :math:`\left<|\Bv|\right>`
   * - *b2m*
     - :math:`\left<\Bv^2\right>`
   * - *bm2*
     - :math:`\max(\Bv^2)`
   * - *j2m*
     - :math:`\left<\jv^2\right>`
   * - *jm2*
     - :math:`\max(\jv^2)`
   * - *abm*
     - :math:`\left<\Av\cdot\Bv\right>`
   * - *abumx*
     - :math:`\left<u_x\Av\cdot\Bv\right>`
   * - *abumy*
     - :math:`\left<u_y\Av\cdot\Bv\right>`
   * - *abumz*
     - :math:`\left<u_z\Av\cdot\Bv\right>`
   * - *abmh*
     - :math:`\left<\Av\cdot\Bv\right>`  (temp)
   * - *abmn*
     - :math:`\left<\Av\cdot\Bv\right>`  (north)
   * - *abms*
     - :math:`\left<\Av\cdot\Bv\right>`  (south)
   * - *abrms*
     - :math:`\left<(\Av\cdot\Bv)^2\right>^{1/2}`
   * - *jbrms*
     - :math:`\left<(\jv\cdot\Bv)^2\right>^{1/2}`
   * - *ajm*
     - :math:`\left<\jv\cdot\Av\right>`
   * - *jbm*
     - :math:`\left<\jv\cdot\Bv\right>`
   * - *hjbm*
     -
   * - *jbmh*
     - :math:`\left<\Jv\cdot\Bv\right>`  (temp)
   * - *jbmn*
     - :math:`\left<\Jv\cdot\Bv\right>`  (north)
   * - *jbms*
     - :math:`\left<\Jv\cdot\Bv\right>`  (south)
   * - *ubm*
     - :math:`\left<\uv\cdot\Bv\right>`
   * - *dubrms*
     - :math:`\left<(\uv-\Bv)^2\right>^{1/2}`
   * - *dobrms*
     - :math:`\left<(\boldsymbol{\omega}-\Bv)^2 \right>^{1/2}`
   * - *uxbxm*
     - :math:`\left<u_xB_x\right>`
   * - *uybxm*
     - :math:`\left<u_yB_x\right>`
   * - *uzbxm*
     - :math:`\left<u_zB_x\right>`
   * - *uxbym*
     - :math:`\left<u_xB_y\right>`
   * - *uybym*
     - :math:`\left<u_yB_y\right>`
   * - *uzbym*
     - :math:`\left<u_zB_y\right>`
   * - *uxbzm*
     - :math:`\left<u_xB_z\right>`
   * - *uybzm*
     - :math:`\left<u_yB_z\right>`
   * - *uzbzm*
     - :math:`\left<u_zB_z\right>`
   * - *cosubm*
     - :math:`\left<\Uv\cdot\Bv/(|\Uv|\,|\Bv|) \right>`
   * - *jxbxm*
     - :math:`\left<j_xB_x\right>`
   * - *jybxm*
     - :math:`\left<j_yB_x\right>`
   * - *jzbxm*
     - :math:`\left<j_zB_x\right>`
   * - *jxbym*
     - :math:`\left<j_xB_y\right>`
   * - *jybym*
     - :math:`\left<j_yB_y\right>`
   * - *jzbym*
     - :math:`\left<j_zB_y\right>`
   * - *jxbzm*
     - :math:`\left<j_xB_z\right>`
   * - *jybzm*
     - :math:`\left<j_yB_z\right>`
   * - *jzbzm*
     - :math:`\left<j_zB_z\right>`
   * - *uam*
     - :math:`\left<\uv\cdot\Av\right>`
   * - *ujm*
     - :math:`\left<\uv\cdot\Jv\right>`
   * - *fbm*
     - :math:`\left<\fv\cdot\Bv\right>`
   * - *fxbxm*
     - :math:`\left<f_x B_x\right>`
   * - *epsM*
     - :math:`\left<\eta\mu_0\jv^2\right>`
   * - *epsAD*
     - :math:`\left<\rho^{-1} t_{\mathrm AD} (\vec{J}\times\vec{B})^2\right>`  (heating by ion-neutrals friction)
   * - *bxpt*
     - :math:`B_x(x_1,y_1,z_1,t)`
   * - *bypt*
     - :math:`B_y(x_1,y_1,z_1,t)`
   * - *bzpt*
     - :math:`B_z(x_1,y_1,z_1,t)`
   * - *jxpt*
     - :math:`J_x(x_1,y_1,z_1,t)`
   * - *jypt*
     - :math:`J_y(x_1,y_1,z_1,t)`
   * - *jzpt*
     - :math:`J_z(x_1,y_1,z_1,t)`
   * - *Expt*
     - :math:`{\cal E}_x(x_1,y_1,z_1,t)`
   * - *Eypt*
     - :math:`{\cal E}_y(x_1,y_1,z_1,t)`
   * - *Ezpt*
     - :math:`{\cal E}_z(x_1,y_1,z_1,t)`
   * - *axpt*
     - :math:`A_x(x_1,y_1,z_1,t)`
   * - *aypt*
     - :math:`A_y(x_1,y_1,z_1,t)`
   * - *azpt*
     - :math:`A_z(x_1,y_1,z_1,t)`
   * - *bxp2*
     - :math:`B_x(x_2,y_2,z_2,t)`
   * - *byp2*
     - :math:`B_y(x_2,y_2,z_2,t)`
   * - *bzp2*
     - :math:`B_z(x_2,y_2,z_2,t)`
   * - *jxp2*
     - :math:`J_x(x_2,y_2,z_2,t)`
   * - *jyp2*
     - :math:`J_y(x_2,y_2,z_2,t)`
   * - *jzp2*
     - :math:`J_z(x_2,y_2,z_2,t)`
   * - *Exp2*
     - :math:`{\cal E}_x(x_2,y_2,z_2,t)`
   * - *Eyp2*
     - :math:`{\cal E}_y(x_2,y_2,z_2,t)`
   * - *Ezp2*
     - :math:`{\cal E}_z(x_2,y_2,z_2,t)`
   * - *axp2*
     - :math:`A_x(x_2,y_2,z_2,t)`
   * - *ayp2*
     - :math:`A_y(x_2,y_2,z_2,t)`
   * - *azp2*
     - :math:`A_z(x_2,y_2,z_2,t)`
   * - *epsM_LES*
     -
   * - *aybym2*
     -
   * - *exaym2*
     -
   * - *exabot*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm bot}`
   * - *exatop*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm top}`
   * - *exjm2*
     -
   * - *emag*
     - :math:`\int_V{1\over2\mu_0}\Bv^2\, dV`
   * - *brms*
     - :math:`\left<\Bv^2\right>^{1/2}`
   * - *bfrms*
     - :math:`\left<{\Bv'}^2\right>^{1/2}`
   * - *bmax*
     - :math:`\max(|\Bv|)`
   * - *bxmin*
     - :math:`\min(|B_x|)`
   * - *bymin*
     - :math:`\min(|B_y|)`
   * - *bzmin*
     - :math:`\min(|B_z|)`
   * - *bxmax*
     - :math:`\max(|B_x|)`
   * - *bymax*
     - :math:`\max(|B_y|)`
   * - *bzmax*
     - :math:`\max(|B_z|)`
   * - *bbxmax*
     - :math:`\max(|B_x|) excluding Bv_{ext}`
   * - *bbymax*
     - :math:`\max(|B_y|) excluding Bv_{ext}`
   * - *bbzmax*
     - :math:`\max(|B_z|) excluding Bv_{ext}`
   * - *jxmax*
     - :math:`\max(|jv_x|)`
   * - *jymax*
     - :math:`\max(|jv_y|)`
   * - *jzmax*
     - :math:`\max(|jv_z|)`
   * - *jrms*
     - :math:`\left<\jv^2\right>^{1/2}`
   * - *hjrms*
     - :math:`\left<\jv^2\right>^{1/2}`
   * - *jmax*
     - :math:`\max(|\jv|)`
   * - *vArms*
     - :math:`\left<\Bv^2/\varrho\right>^{1/2}`
   * - *vAmax*
     - :math:`\max(\Bv^2/\varrho)^{1/2}`
   * - *dtb*
     - :math:`\delta t / [c_{\delta t}\,\delta x /v_{\mathrm A,max}]`  :math:`\quad`  (time step relative to Alfv{\'e}n time step; see \S~\ref{time-step})
   * - *dteta*
     - :math:`\delta t/[c_{\delta t,{\mathrm v}}\, \delta x^2/\eta_{\mathrm max}]`  :math:`\quad`  (time step relative to resistive time step; see \S~\ref{time-step})
   * - *axm*
     -
   * - *aym*
     -
   * - *azm*
     -
   * - *a2m*
     - :math:`\left<\Av^2\right>`
   * - *arms*
     - :math:`\left<\Av^2\right>^{1/2}`
   * - *amax*
     - :math:`\max(|\Av|)`
   * - *divarms*
     - :math:`\langle(\nabla\cdot\Av)^2\rangle^{1/2}`
   * - *beta1m*
     - :math:`\left<\Bv^2/(2\mu_0 p)\right>`  :math:`\quad`  (mean inverse plasma beta)
   * - *beta1max*
     - :math:`\max[\Bv^2/(2\mu_0 p)]`  :math:`\quad`  (maximum inverse plasma beta)
   * - *betam*
     - :math:`\langle\beta\rangle`
   * - *betamax*
     - :math:`\max\beta`
   * - *betamin*
     - :math:`\min\beta`
   * - *bxm*
     - :math:`\left<B_x\right>`
   * - *bym*
     - :math:`\left<B_y\right>`
   * - *bzm*
     - :math:`\left<B_z\right>`
   * - *bxbym*
     - :math:`\left<B_x B_y\right>`
   * - *bxbzm*
     -
   * - *bybzm*
     -
   * - *djuidjbim*
     -
   * - *bmx*
     - :math:`\left<\left<\Bv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean field)
   * - *bmy*
     - :math:`\left<\left<\Bv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean field)
   * - *bmz*
     - :math:`\left<\left<\Bv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean field)
   * - *bmzS2*
     - :math:`\left<\left<\Bv_S\right>_{xy}^2\right>`
   * - *bmzA2*
     - :math:`\left<\left<\Bv_A\right>_{xy}^2\right>`
   * - *jmx*
     - :math:`\left<\left<\Jv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean current density)
   * - *jmy*
     - :math:`\left<\left<\Jv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean current density)
   * - *jmz*
     - :math:`\left<\left<\Jv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean current density)
   * - *bmzph*
     - Phase of a Beltrami field
   * - *bmzphe*
     - Error of phase of a Beltrami field
   * - *bsinphz*
     - sine of phase of a Beltrami field
   * - *bcosphz*
     - cosine of phase of a Beltrami field
   * - *emxamz3*
     - :math:`\left<\left<\Ev\right>_{xy}\times\left<\Av\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity flux)
   * - *embmz*
     - :math:`\left<\left<\Ev\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity production )
   * - *ambmz*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field)
   * - *ambmzh*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, temp)
   * - *ambmzn*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, north)
   * - *ambmzs*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, south)
   * - *jmbmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (current helicity of :math:`xy` -averaged mean field)
   * - *kx_aa*
     - :math:`k_x`
   * - *kmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>/ \left<\left<\Bv\right>_{xy}^2\right>`
   * - *bx2m*
     - :math:`\left< B_x^2 \right>`
   * - *by2m*
     - :math:`\left< B_y^2 \right>`
   * - *bz2m*
     - :math:`\left< B_z^2 \right>`
   * - *uxbm*
     - :math:`\left<\uv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *jxbm*
     - :math:`\left<\jv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *magfricmax*
     - Magneto-Frictional velocity :math:`\left<\jv\times\Bv\right>\cdot\Bv^2`
   * - *oxuxbm*
     -
   * - *jxbxbm*
     -
   * - *gpxbm*
     -
   * - *uxDxuxbm*
     -
   * - *b3b21m*
     - :math:`\left<B_3 B_{2,1} \right>`
   * - *b3b12m*
     - :math:`\left<B_3 B_{1,2} \right>`
   * - *b1b32m*
     - :math:`\left<B_1 B_{3,2} \right>`
   * - *b1b23m*
     - :math:`\left<B_1 B_{2,3} \right>`
   * - *b2b13m*
     - :math:`\left<B_2 B_{1,3} \right>`
   * - *b2b31m*
     - :math:`\left<B_2 B_{3,1} \right>`
   * - *udotxbm*
     -
   * - *uxbdotm*
     -
   * - *uxbmx*
     - :math:`\left<(\uv\times\Bv)_x\right>`
   * - *uxbmy*
     - :math:`\left<(\uv\times\Bv)_y\right>`
   * - *uxbmz*
     - :math:`\left<(\uv\times\Bv)_z\right>`
   * - *jxbmx*
     - :math:`\left<(\jv\times\Bv)_x\right>`
   * - *jxbmy*
     - :math:`\left<(\jv\times\Bv)_y\right>`
   * - *jxbmz*
     - :math:`\left<(\jv\times\Bv)_z\right>`
   * - *uxbcmx*
     -
   * - *uxbcmy*
     -
   * - *uxbsmx*
     -
   * - *uxbsmy*
     -
   * - *examx*
     - :math:`\left<\Ev\times\Av\right>|_x`
   * - *examy*
     - :math:`\left<\Ev\times\Av\right>|_y`
   * - *examz*
     - :math:`\left<\Ev\times\Av\right>|_z`
   * - *exjmx*
     - :math:`\left<\Ev\times\Jv\right>|_x`
   * - *exjmy*
     - :math:`\left<\Ev\times\Jv\right>|_y`
   * - *exjmz*
     - :math:`\left<\Ev\times\Jv\right>|_z`
   * - *dexbmx*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_x`
   * - *dexbmy*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_y`
   * - *dexbmz*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_z`
   * - *phibmx*
     - :math:`\left<\phi\Bv\right>|_x`
   * - *phibmy*
     - :math:`\left<\phi\Bv\right>|_y`
   * - *phibmz*
     - :math:`\left<\phi\Bv\right>|_z`
   * - *uxjm*
     -
   * - *b2divum*
     - :math:`\left<\Bv^2\nabla\cdot\uv\right>`
   * - *ujxbm*
     - :math:`\left<\uv\cdot(\Jv\times\Bv)\right>`
   * - *jxbrxm*
     -
   * - *jxbrym*
     -
   * - *jxbrzm*
     -
   * - *jxbrmax*
     - :math:`\max(|\Jv\times\Bv/\rho|)`
   * - *jxbr2m*
     - :math:`\left<(\Jv\times\Bv/\rho)^2\right>`
   * - *uxBrms*
     -
   * - *Bresrms*
     -
   * - *Rmrms*
     -
   * - *jfm*
     -
   * - *brbpmr*
     -
   * - *vA2m*
     -
   * - *b2mr*
     -
   * - *brmr*
     -
   * - *bpmr*
     -
   * - *bzmr*
     -
   * - *armr*
     -
   * - *apmr*
     -
   * - *azmr*
     -
   * - *mflux_x*
     -
   * - *mflux_y*
     -
   * - *mflux_z*
     -
   * - *bmxy_rms*
     - :math:`\sqrt{[\left<b_x\right>_z(x,y)]^2 + [\left<b_y\right>_z(x,y)]^2 + [\left<b_z\right>_z(x,y)]^2}`
   * - *etasmagm*
     - Mean of Smagorinsky resistivity
   * - *etasmagmin*
     - Min of Smagorinsky resistivity
   * - *etasmagmax*
     - Max of Smagorinsky resistivity
   * - *etavamax*
     - Max of artificial resistivity :math:`\eta\sim v_A`
   * - *etajmax*
     - Max of artificial resistivity :math:`\eta\sim J / \sqrt{\rho}`
   * - *etaj2max*
     - Max of artificial resistivity :math:`\eta\sim J^2 / \rho`
   * - *etajrhomax*
     - Max of artificial resistivity :math:`\eta\sim J / \rho`
   * - *cosjbm*
     - :math:`\left<\Jv\cdot\Bv/(|\Jv|\,|\Bv|)\right>`
   * - *coshjbm*
     -
   * - *jparallelm*
     - Mean value of the component of J parallel to B
   * - *jperpm*
     - Mean value of the component of J perpendicular to B
   * - *hjparallelm*
     - Mean value of the component of :math:`J_{\mathrm hyper}`  parallel to B
   * - *hjperpm*
     - Mean value of the component of :math:`J_{\mathrm hyper}`  perpendicular to B
   * - *brmsx*
     - :math:`\left<\Bv^2\right>^{1/2}`  for the magnetic_xaver_range
   * - *brmsz*
     - :math:`\left<\Bv^2\right>^{1/2}`  for the magnetic_zaver_range
   * - *Exmxy*
     - :math:`\left<{\cal E}_x\right>_{z}`
   * - *Eymxy*
     - :math:`\left<{\cal E}_y\right>_{z}`
   * - *Ezmxy*
     - :math:`\left<{\cal E}_z\right>_{z}`

Module *experimental/particles_dust_brdeplete.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *xpm*
     - :math:`x_{part}`
   * - *xp2m*
     - :math:`x^2_{part}`
   * - *vrelpabsm*
     - :math:`\mathrm{Absolute value of mean relative velocity}`
   * - *vpxm*
     - :math:`u_{part}`
   * - *vpx2m*
     - :math:`u^2_{part}`
   * - *ekinp*
     - :math:`E_{kin,part}`
   * - *vpxmax*
     - :math:`MAX(u_{part})`
   * - *vpxmin*
     - :math:`MIN(u_{part})`
   * - *npm*
     - :math:`\mathrm{mean particle number density}`

Module *experimental/solid_cells_CGEO.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *c_dragx*
     -
   * - *c_dragy*
     -
   * - *c_dragz*
     -
   * - *c_dragx_p*
     -
   * - *c_dragy_p*
     -
   * - *c_dragz_p*
     -
   * - *Nusselt*
     -

Module *experimental/solid_cells_reactive.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *c_dragx*
     -
   * - *c_dragy*
     -
   * - *c_dragz*
     -
   * - *c_dragx_p*
     -
   * - *c_dragy_p*
     -
   * - *c_dragz_p*
     -
   * - *Nusselt*
     -

Module *experimental/test_chemistry.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Ym*
     - :math:`\left<Y_x\right>`
   * - *dYm*
     - :math:`\delta\left<Y_x\right>/\delta t`
   * - *dYmax*
     - :math:`max\delta\left<Y_x\right>/\delta t`
   * - *Ymax*
     - :math:`\left<Y_{x,max}\right>`
   * - *Ymin*
     - :math:`\left<Y_{x,min}\right>`
   * - *hm*
     - :math:`\left<H_{x,max}\right>`
   * - *cpm*
     - :math:`\left<c_{p,x}\right>`
   * - *diffm*
     - :math:`\left<D_{x}\right>`
   * - *Ymz*
     - :math:`\left<Y_x\right>_{xy}(z)`
   * - *dtchem*
     - :math:`dt_{chem}`

Module *experimental/testperturb.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *alp31*
     - :math:`\alpha_{31}`
   * - *alp12*
     - :math:`\alpha_{12}`
   * - *alp22*
     - :math:`\alpha_{22}`
   * - *alp32*
     - :math:`\alpha_{32}`
   * - *eta11*
     - :math:`\eta_{113}k`
   * - *eta21*
     - :math:`\eta_{213}k`
   * - *eta31*
     - :math:`\eta_{313}k`
   * - *eta12*
     - :math:`\eta_{123}k`
   * - *eta22*
     - :math:`\eta_{223}k`
   * - *eta32*
     - :math:`\eta_{323}k`

Module *inactive/gravitational_waves_hTXk_no_xpara.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *g11pt*
     - :math:`g_{11}(x_1,y_1,z_1,t)`
   * - *g22pt*
     - :math:`g_{22}(x_1,y_1,z_1,t)`
   * - *g33pt*
     - :math:`g_{33}(x_1,y_1,z_1,t)`
   * - *g12pt*
     - :math:`g_{12}(x_1,y_1,z_1,t)`
   * - *g23pt*
     - :math:`g_{23}(x_1,y_1,z_1,t)`
   * - *g31pt*
     - :math:`g_{31}(x_1,y_1,z_1,t)`
   * - *hhTpt*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXpt*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTpt*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXpt*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hhTp2*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXp2*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTp2*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXp2*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hrms*
     - :math:`\bra{h_T^2+h_X^2}^{1/2}`
   * - *EEGW*
     - :math:`\bra{g_T^2+g_X^2}\,c^2/(32\pi G)`
   * - *gg2m*
     - :math:`\bra{g_T^2+g_X^2}`
   * - *hhT2m*
     - :math:`\bra{h_T^2}`
   * - *hhX2m*
     - :math:`\bra{h_X^2}`
   * - *hhTXm*
     - :math:`\bra{h_T h_X}`
   * - *ggT2m*
     - :math:`\bra{g_T^2}`
   * - *ggX2m*
     - :math:`\bra{g_X^2}`
   * - *ggTXm*
     - :math:`\bra{g_T g_X}`

Module *inactive/testfield_xz.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *E111z*
     - :math:`{\cal E}_1^{11}`
   * - *E211z*
     - :math:`{\cal E}_2^{11}`
   * - *E311z*
     - :math:`{\cal E}_3^{11}`
   * - *E121z*
     - :math:`{\cal E}_1^{21}`
   * - *E221z*
     - :math:`{\cal E}_2^{21}`
   * - *E321z*
     - :math:`{\cal E}_3^{21}`
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *eta11*
     - :math:`\eta_{113}k`
   * - *eta21*
     - :math:`\eta_{213}k`
   * - *b11rms*
     - :math:`\left<b_{11}^2\right>`
   * - *b21rms*
     - :math:`\left<b_{21}^2\right>`

Module *magnetic/maxwell.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *aa2m*
     - :math:`\langle A^2\rangle`
   * - *ee2m*
     - :math:`\langle E^2\rangle`
   * - *EEEM*
     - :math:`\langle(E^2+B^2)/2\rangle`
   * - *akxpt*
     - :math:`Akx^{pt}`
   * - *ekxpt*
     - :math:`Ekx^{pt}`
   * - *sigma*
     - :math:`\sigma`
   * - *emag*
     - :math:`\int_V{1\over2\mu_0}\Bv^2\, dV`
   * - *bmax*
     - :math:`\max(|\Bv|)`
   * - *brms*
     - :math:`\left<\Bv^2\right>^{1/2}`
   * - *arms*
     - :math:`\left<\Av^2\right>^{1/2}`
   * - *erms*
     - :math:`\left<\Ev^2\right>^{1/2}`
   * - *bfrms*
     - :math:`\left<{\Bv'}^2\right>^{1/2}`

Module *magnetic/meanfield.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *qsm*
     - :math:`\left<q_p(\overline{B})\right>`
   * - *qpm*
     - :math:`\left<q_p(\overline{B})\right>`
   * - *qem*
     - :math:`\left<q_e(\overline{B})\right>` , in the paper referred to as :math:`\left<q_g(\overline{B})\right>`
   * - *qam*
     - :math:`\left<q_a(\overline{B})\right>`
   * - *alpm*
     - :math:`\left<\alpha\right>`   !(where is this implemented?)
   * - *etatm*
     - :math:`\left<\eta_{\mathrm t}\right>`
   * - *EMFmz1*
     - :math:`\left<{\cal E}\right>_{xy}|_x`
   * - *EMFmz2*
     - :math:`\left<{\cal E}\right>_{xy}|_y`
   * - *EMFmz3*
     - :math:`\left<{\cal E}\right>_{xy}|_z`
   * - *EMFdotBm*
     - :math:`\left<{\cal E}\cdot\Bv \right>`
   * - *EMFdotB_int*
     - :math:`\int{\cal E}\cdot\Bv dV`
   * - *alpKjbm*
     - :math:`\left<\alpha_\mathrm{K}\overline{\Bv}\cdot\overline{\Jv}\right>`
   * - *alpKm*
     - :math:`\left<\alpha_\mathrm{K}\right>`

Module *magnetic/meanfield_demfdt.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *EMFrms*
     - :math:`(\left<{\cal E}\right>)_{\mathrm rms}`
   * - *EMFmax*
     - :math:`\max(\left<{\cal E}\right>)`
   * - *EMFmin*
     - :math:`\min(\left<{\cal E}\right>)`

Module *obsolete/density_anelastic.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rhom*
     - :math:`\left<\varrho\right>`  :math:`\quad`  (mean density)
   * - *rho2m*
     -
   * - *lnrho2m*
     -
   * - *drho2m*
     -
   * - *drhom*
     -
   * - *rhomin*
     -
   * - *rhomax*
     -
   * - *ugrhom*
     - :math:`\left<\uv\cdot\nabla\varrho\right>`
   * - *uglnrhom*
     -
   * - *dtd*
     -
   * - *rhomz*
     -
   * - *rhomy*
     -
   * - *rhomx*
     -
   * - *rhomxy*
     -
   * - *rhomxz*
     -
   * - *rhomr*
     -
   * - *totmass*
     -
   * - *mass*
     - :math:`\int\varrho\,dV`
   * - *divrhoum*
     - :math:`\left<\nabla\cdot(\varrho\uv)\right>`
   * - *divrhourms*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm rms}`
   * - *divrhoumax*
     - :math:`\left|\nabla\cdot(\varrho\uv)\right|_{\mathrm max}`

Module *obsolete/magnetic_axisym.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *aphi2m*
     - :math:`\left<A_\phi^2\right>`
   * - *bphi2m*
     - :math:`\left<B_\phi^2\right>`
   * - *bpol2m*
     - :math:`\left<B_p^2\right>`

Module *obsolete/magnetic_old.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *ab_int*
     - :math:`\int\Av\cdot\Bv\;dV`
   * - *jb_int*
     - :math:`\int\jv\cdot\Bv\;dV`
   * - *b2tm*
     - :math:`\left<\bv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *bjtm*
     - :math:`\left<\bv(t)\cdot\int_0^t\jv(t') dt'\right>`
   * - *jbtm*
     - :math:`\left<\jv(t)\cdot\int_0^t\bv(t') dt'\right>`
   * - *b2ruzm*
     - :math:`\left<\Bv^2\rho u_z\right>`
   * - *b2uzm*
     - :math:`\left<\Bv^2u_z\right>`
   * - *ubbzm*
     - :math:`\left<(\uv\cdot\Bv)B_z\right>`
   * - *b1m*
     - :math:`\left<|\Bv|\right>`
   * - *b2m*
     - :math:`\left<\Bv^2\right>`
   * - *bm2*
     - :math:`\max(\Bv^2)`
   * - *j2m*
     - :math:`\left<\jv^2\right>`
   * - *jm2*
     - :math:`\max(\jv^2)`
   * - *abm*
     - :math:`\left<\Av\cdot\Bv\right>`
   * - *abumx*
     - :math:`\left<u_x\Av\cdot\Bv\right>`
   * - *abumy*
     - :math:`\left<u_y\Av\cdot\Bv\right>`
   * - *abumz*
     - :math:`\left<u_z\Av\cdot\Bv\right>`
   * - *abmh*
     - :math:`\left<\Av\cdot\Bv\right>`  (temp)
   * - *abmn*
     - :math:`\left<\Av\cdot\Bv\right>`  (north)
   * - *abms*
     - :math:`\left<\Av\cdot\Bv\right>`  (south)
   * - *abrms*
     - :math:`\left<(\Av\cdot\Bv)^2\right>^{1/2}`
   * - *ajm*
     - :math:`\left<\jv\cdot\Av\right>`
   * - *jbm*
     - :math:`\left<\jv\cdot\Bv\right>`
   * - *jbmh*
     - :math:`\left<\Av\cdot\Bv\right>`  (temp)
   * - *jbmn*
     - :math:`\left<\Av\cdot\Bv\right>`  (north)
   * - *jbms*
     - :math:`\left<\Av\cdot\Bv\right>`  (south)
   * - *ubm*
     - :math:`\left<\uv\cdot\Bv\right>`
   * - *dubrms*
     - :math:`\left<(\uv-\Bv)^2\right>^{1/2}`
   * - *dobrms*
     - :math:`\left<(\boldsymbol{\omega}-\Bv)^2\right>^{1/2}`
   * - *uxbxm*
     - :math:`\left<u_xB_x\right>`
   * - *uybym*
     - :math:`\left<u_yB_y\right>`
   * - *uzbzm*
     - :math:`\left<u_zB_z\right>`
   * - *cosubm*
     - :math:`\left<\Uv\cdot\Bv/(|\Uv|\,|\Bv|)\right>`
   * - *uam*
     - :math:`\left<\uv\cdot\Av\right>`
   * - *ujm*
     - :math:`\left<\uv\cdot\Jv\right>`
   * - *fbm*
     - :math:`\left<\fv\cdot\Bv\right>`
   * - *fxbxm*
     - :math:`\left<f_x B_x\right>`
   * - *epsM*
     - :math:`\left<2\eta\mu_0\jv^2\right>`
   * - *epsAD*
     - :math:`\left<\rho^{-1} t_{\mathrm AD} (\vec{J}\times\vec{B})^2\right>`  (heating by ion-neutrals friction)
   * - *bxpt*
     - :math:`B_x(x_0,y_0,z_0,t)`
   * - *bypt*
     - :math:`B_y(x_0,y_0,z_0,t)`
   * - *bzpt*
     - :math:`B_z(x_0,y_0,z_0,t)`
   * - *Expt*
     - :math:`{\cal E}_x(x_0,y_0,z_0,t)`
   * - *Eypt*
     - :math:`{\cal E}_y(x_0,y_0,z_0,t)`
   * - *Ezpt*
     - :math:`{\cal E}_z(x_0,y_0,z_0,t)`
   * - *epsM_LES*
     -
   * - *aybym2*
     -
   * - *exaym2*
     -
   * - *exabot*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm bot}`
   * - *exatop*
     - :math:`\int\Ev\times\Av\,dS|_{\mathrm top}`
   * - *exjm2*
     -
   * - *emag*
     - :math:`\int_V{1\over2\mu_0}\Bv^2\, dV`
   * - *brms*
     - :math:`\left<\Bv^2\right>^{1/2}`
   * - *bmax*
     - :math:`\max(|\Bv|)`
   * - *bxmin*
     - :math:`\min(|B_x|)`
   * - *bymin*
     - :math:`\min(|B_y|)`
   * - *bzmin*
     - :math:`\min(|B_z|)`
   * - *bxmax*
     - :math:`\max(|B_x|)`
   * - *bymax*
     - :math:`\max(|B_y|)`
   * - *bzmax*
     - :math:`\max(|B_z|)`
   * - *jrms*
     - :math:`\left<\jv^2\right>^{1/2}`
   * - *jmax*
     - :math:`\max(|\jv|)`
   * - *vArms*
     - :math:`\left<\Bv^2/\varrho\right>^{1/2}`
   * - *vAmax*
     - :math:`\max(\Bv^2/\varrho)^{1/2}`
   * - *dtb*
     - :math:`\delta t / [c_{\delta t}\,\delta x /v_{\mathrm A,max}]`  :math:`\quad`  (time step relative to Alfv{\'e}n time step; see \S~\ref{time-step})
   * - *dteta*
     - :math:`\delta t/[c_{\delta t,{\mathrm v}}\, \delta x^2/\eta_{\mathrm max}]`  :math:`\quad`  (time step relative to resistive time step; see \S~\ref{time-step})
   * - *axm*
     -
   * - *aym*
     -
   * - *azm*
     -
   * - *arms*
     -
   * - *amax*
     -
   * - *beta1m*
     - :math:`\left<\Bv^2/(2\mu_0 p)\right>`  :math:`\quad`  (mean inverse plasma beta)
   * - *beta1max*
     - :math:`\max[\Bv^2/(2\mu_0 p)]`  :math:`\quad`  (maximum inverse plasma beta)
   * - *bxm*
     - :math:`\left<\left<B\right>_{yz}^2\right>_x^{1/2}`
   * - *bym*
     - :math:`\left<\left<B\right>_{zx}^2\right>_y^{1/2}`
   * - *bzm*
     - :math:`\left<\left<B\right>_{xy}^2\right>_z^{1/2}`
   * - *bx2m*
     - :math:`\left<B_x^2\right>`
   * - *by2m*
     - :math:`\left<B_y^2\right>`
   * - *bz2m*
     - :math:`\left<B_z^2\right>`
   * - *bx2mx*
     - :math:`\left<B_x^2\right>_{yz}`
   * - *by2mx*
     - :math:`\left<B_y^2\right>_{yz}`
   * - *bz2mx*
     - :math:`\left<B_z^2\right>_{yz}`
   * - *bxbym*
     - :math:`\left<B_x B_y\right>`
   * - *bxbymx*
     - :math:`\left<B_x B_y\right>_{yz}`
   * - *bxbzm*
     -
   * - *bybzm*
     -
   * - *djuidjbim*
     -
   * - *bxbymy*
     -
   * - *bxbzmy*
     -
   * - *bybzmy*
     -
   * - *bxbymz*
     -
   * - *bxbzmz*
     -
   * - *bybzmz*
     -
   * - *b2mz*
     -
   * - *axmz*
     - :math:`\left<{\cal A}_x\right>_{xy}`
   * - *aymz*
     - :math:`\left<{\cal A}_y\right>_{xy}`
   * - *azmz*
     - :math:`\left<{\cal A}_z\right>_{xy}`
   * - *abuxmz*
     - :math:`\left<(\Av \cdot \Bv) u_x \right>_{xy}`
   * - *abuymz*
     - :math:`\left<(\Av \cdot \Bv) u_y \right>_{xy}`
   * - *abuzmz*
     - :math:`\left<(\Av \cdot \Bv) u_z \right>_{xy}`
   * - *uabxmz*
     - :math:`\left<(\uv \cdot \Av) B_x \right>_{xy}`
   * - *uabymz*
     - :math:`\left<(\uv \cdot \Av) B_y \right>_{xy}`
   * - *uabzmz*
     - :math:`\left<(\uv \cdot \Av) B_z \right>_{xy}`
   * - *bbxmz*
     - :math:`\left<{\cal B}'_x\right>_{xy}`
   * - *bbymz*
     - :math:`\left<{\cal B}'_y\right>_{xy}`
   * - *bbzmz*
     - :math:`\left<{\cal B}'_z\right>_{xy}`
   * - *bxmz*
     - :math:`\left<{\cal B}_x\right>_{xy}`
   * - *bymz*
     - :math:`\left<{\cal B}_y\right>_{xy}`
   * - *bzmz*
     - :math:`\left<{\cal B}_z\right>_{xy}`
   * - *jxmz*
     - :math:`\left<{\cal J}_x\right>_{xy}`
   * - *jymz*
     - :math:`\left<{\cal J}_y\right>_{xy}`
   * - *jzmz*
     - :math:`\left<{\cal J}_z\right>_{xy}`
   * - *Exmz*
     - :math:`\left<{\cal E}_x\right>_{xy}`
   * - *Eymz*
     - :math:`\left<{\cal E}_y\right>_{xy}`
   * - *Ezmz*
     - :math:`\left<{\cal E}_z\right>_{xy}`
   * - *bmx*
     - :math:`\left<\left<\Bv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean field)
   * - *bmy*
     - :math:`\left<\left<\Bv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean field)
   * - *bmz*
     - :math:`\left<\left<\Bv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean field)
   * - *bmzS2*
     - :math:`\left<\left<\Bv_S\right>_{xy}^2\right>`
   * - *bmzA2*
     - :math:`\left<\left<\Bv_A\right>_{xy}^2\right>`
   * - *jmx*
     - :math:`\left<\left<\Jv\right>_{yz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`yz` -averaged mean current density)
   * - *jmy*
     - :math:`\left<\left<\Jv\right>_{xz}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xz` -averaged mean current density)
   * - *jmz*
     - :math:`\left<\left<\Jv\right>_{xy}^2 \right>^{1/2}`  :math:`\quad`  (energy of :math:`xy` -averaged mean current density)
   * - *bmzph*
     - Phase of a Beltrami field
   * - *bmzphe*
     - Error of phase of a Beltrami field
   * - *bsinphz*
     - sine of phase of a Beltrami field
   * - *bcosphz*
     - cosine of phase of a Beltrami field
   * - *emxamz3*
     - :math:`\left<\left<\Ev\right>_{xy}\times\left<\Av\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity flux)
   * - *embmz*
     - :math:`\left<\left<\Ev\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (:math:`xy` -averaged mean field helicity production )
   * - *ambmz*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field)
   * - *ambmzh*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, temp)
   * - *ambmzn*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, north)
   * - *ambmzs*
     - :math:`\left<\left<\Av\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>`  :math:`\quad`   (magnetic helicity of :math:`xy` -averaged mean field, south)
   * - *jmbmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy} \right>`  :math:`\quad`  (current helicity of :math:`xy` -averaged mean field)
   * - *kx_aa*
     - :math:`k_x`
   * - *kmz*
     - :math:`\left<\left<\Jv\right>_{xy}\cdot\left<\Bv\right>_{xy}\right>/ \left<\left<\Bv\right>_{xy}^2\right>`
   * - *bx2my*
     - :math:`\left< B_x^2 \right>_{xz}`
   * - *by2my*
     - :math:`\left< B_y^2 \right>_{xz}`
   * - *bz2my*
     - :math:`\left< B_z^2 \right>_{xz}`
   * - *bx2mz*
     - :math:`\left< B_x^2 \right>_{xy}`
   * - *by2mz*
     - :math:`\left< B_y^2 \right>_{xy}`
   * - *bz2mz*
     - :math:`\left< B_z^2 \right>_{xy}`
   * - *bx2rmz*
     - :math:`\left< B_x^2/\varrho \right>_{xy}`
   * - *by2rmz*
     - :math:`\left< B_y^2/\varrho \right>_{xy}`
   * - *bz2rmz*
     - :math:`\left< B_z^2/\varrho \right>_{xy}`
   * - *jbmz*
     - :math:`\left<\Jv\cdot\Bv\right>|_{xy}`
   * - *abmz*
     - :math:`\left<\Av\cdot\Bv\right>|_{xy}`
   * - *uamz*
     - :math:`\left<\uv\cdot\Av\right>|_{xy}`
   * - *examz1*
     - :math:`\left<\Ev\times\Av\right>_{xy}|_x`
   * - *examz2*
     - :math:`\left<\Ev\times\Av\right>_{xy}|_y`
   * - *examz3*
     - :math:`\left<\Ev\times\Av\right>_{xy}|_z`
   * - *bxmxy*
     - :math:`\left< B_x \right>_{xy}`
   * - *bymxy*
     - :math:`\left< B_y \right>_{xy}`
   * - *bzmxy*
     - :math:`\left< B_z \right>_{xy}`
   * - *jxmxy*
     - :math:`\left< J_x \right>_{xy}`
   * - *jymxy*
     - :math:`\left< J_y \right>_{xy}`
   * - *jzmxy*
     - :math:`\left< J_z \right>_{xy}`
   * - *axmxz*
     - :math:`\left< A_x \right>_{xz}`
   * - *aymxz*
     - :math:`\left< A_y \right>_{xz}`
   * - *azmxz*
     - :math:`\left< A_z \right>_{xz}`
   * - *bxmxz*
     - :math:`\left< B_x \right>_{xz}`
   * - *bymxz*
     - :math:`\left< B_y \right>_{xz}`
   * - *bzmxz*
     - :math:`\left< B_z \right>_{xz}`
   * - *bx2mxz*
     - :math:`\left< B_x^2 \right>_{xz}`
   * - *by2mxz*
     - :math:`\left< B_y^2 \right>_{xz}`
   * - *bz2mxz*
     - :math:`\left< B_z^2 \right>_{xz}`
   * - *bx2mxy*
     - :math:`\left< B_x^2 \right>_{z}`
   * - *by2mxy*
     - :math:`\left< B_y^2 \right>_{z}`
   * - *bz2mxy*
     - :math:`\left< B_z^2 \right>_{z}`
   * - *bxbymxy*
     - :math:`\left< B_x B_y \right>_{z}`
   * - *bxbzmxy*
     - :math:`\left< B_x B_z \right>_{z}`
   * - *bybzmxy*
     - :math:`\left< B_y B_z \right>_{z}`
   * - *bxbymxz*
     - :math:`\left< B_x B_y \right>_{y}`
   * - *bxbzmxz*
     - :math:`\left< B_x B_z \right>_{y}`
   * - *bybzmxz*
     - :math:`\left< B_y B_z \right>_{y}`
   * - *uxbm*
     - :math:`\left<\uv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *jxbm*
     - :math:`\left<\jv\times\Bv\right>\cdot\Bv_0/B_0^2`
   * - *oxuxbm*
     -
   * - *jxbxbm*
     -
   * - *gpxbm*
     -
   * - *uxDxuxbm*
     -
   * - *b3b21m*
     - :math:`\left<B_3 B_{2,1} \right>`
   * - *b1b32m*
     - :math:`\left<B_1 B_{3,2} \right>`
   * - *b2b13m*
     - :math:`\left<B_2 B_{1,3} \right>`
   * - *udotxbm*
     -
   * - *uxbdotm*
     -
   * - *uxbmx*
     - :math:`\left<(\uv\times\Bv)_x\right>`
   * - *uxbmy*
     - :math:`\left<(\uv\times\Bv)_y\right>`
   * - *uxbmz*
     - :math:`\left<(\uv\times\Bv)_z\right>`
   * - *jxbmx*
     - :math:`\left<(\jv\times\Bv)_x\right>`
   * - *jxbmy*
     - :math:`\left<(\jv\times\Bv)_y\right>`
   * - *jxbmz*
     - :math:`\left<(\jv\times\Bv)_z\right>`
   * - *uxbcmx*
     -
   * - *uxbcmy*
     -
   * - *uxbsmx*
     -
   * - *uxbsmy*
     -
   * - *examx*
     - :math:`\left<\Ev\times\Av\right>|_x`
   * - *examy*
     - :math:`\left<\Ev\times\Av\right>|_y`
   * - *examz*
     - :math:`\left<\Ev\times\Av\right>|_z`
   * - *exjmx*
     - :math:`\left<\Ev\times\Jv\right>|_x`
   * - *exjmy*
     - :math:`\left<\Ev\times\Jv\right>|_y`
   * - *exjmz*
     - :math:`\left<\Ev\times\Jv\right>|_z`
   * - *dexbmx*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_x`
   * - *dexbmy*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_y`
   * - *dexbmz*
     - :math:`\left<\nabla\times\Ev\times\Bv\right>|_z`
   * - *phibmx*
     - :math:`\left<\phi\Bv\right>|_x`
   * - *phibmy*
     - :math:`\left<\phi\Bv\right>|_y`
   * - *phibmz*
     - :math:`\left<\phi\Bv\right>|_z`
   * - *uxjm*
     -
   * - *b2divum*
     - :math:`\left<\Bv^2\nabla\cdot\uv\right>`
   * - *ujxbm*
     - :math:`\left<\uv\cdot(\Jv\times\Bv\right>`
   * - *jxbrxm*
     -
   * - *jxbrym*
     -
   * - *jxbrzm*
     -
   * - *jxbr2m*
     - :math:`\left<(\Jv\times\Bv/\rho)^2\right>`
   * - *jxbrxmx*
     -
   * - *jxbrymx*
     -
   * - *jxbrzmx*
     -
   * - *jxbrxmy*
     -
   * - *jxbrymy*
     -
   * - *jxbrzmy*
     -
   * - *jxbrxmz*
     -
   * - *jxbrymz*
     -
   * - *jxbrzmz*
     -
   * - *uxBrms*
     -
   * - *Bresrms*
     -
   * - *Rmrms*
     -
   * - *jfm*
     -
   * - *brbpmr*
     -
   * - *vA2m*
     -
   * - *b2mr*
     -
   * - *brmr*
     -
   * - *bpmr*
     -
   * - *bzmr*
     -
   * - *armr*
     -
   * - *apmr*
     -
   * - *azmr*
     -
   * - *bxmx*
     -
   * - *bymx*
     -
   * - *bzmx*
     -
   * - *bxmy*
     -
   * - *bymy*
     -
   * - *bzmy*
     -
   * - *mflux_x*
     -
   * - *mflux_y*
     -
   * - *mflux_z*
     -
   * - *bmxy_rms*
     - :math:`\sqrt{[\left<b_x\right>_z(x,y)]^2 + [\left<b_y\right>_z(x,y)]^2 + [\left<b_z\right>_z(x,y)]^2}`
   * - *etasmagm*
     - Mean of Smagorinsky resistivity
   * - *etasmagmin*
     - Min of Smagorinsky resistivity
   * - *etasmagmax*
     - Max of Smagorinsky resistivity
   * - *etavamax*
     - Max of artificial resistivity :math:`\eta\sim v_A`
   * - *etajmax*
     - Max of artificial resistivity :math:`\eta\sim J / \sqrt{\rho}`
   * - *etaj2max*
     - Max of artificial resistivity :math:`\eta\sim J^2 / \rho`
   * - *etajrhomax*
     - Max of artificial resistivity :math:`\eta\sim J / \rho`
   * - *cosjbm*
     - :math:`\left<\Jv\cdot\Bv/(|\Jv|\,|\Bv|)\right>`
   * - *jparallelm*
     - Mean value of the component of J parallel to B
   * - *jperpm*
     - Mean value of the component of J perpendicular to B
   * - *Exmxz*
     - :math:`\left<{\cal E}_x\right>_{y}`
   * - *Eymxz*
     - :math:`\left<{\cal E}_y\right>_{y}`
   * - *Ezmxz*
     - :math:`\left<{\cal E}_z\right>_{y}`
   * - *etatotalmx*
     - :math:`\left<\eta\right>_{yz}`
   * - *etatotalmz*
     - :math:`\left<\eta\right>_{xy}`

Module *obsolete/special/axionSU2back.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *Q*
     - :math:`Q`
   * - *Qdot*
     - :math:`\dot{Q}`
   * - *Qddot*
     - :math:`\ddot{Q}`
   * - *chi*
     - :math:`\chi`
   * - *chidot*
     - :math:`\dot{\chi}`
   * - *chiddot*
     - :math:`\ddot{\chi}`
   * - *psi*
     - :math:`\psi`
   * - *psiL*
     - :math:`\psi_L`
   * - *psidot*
     - :math:`\dot\psi`
   * - *psiddot*
     - :math:`\ddot\psi`
   * - *TR*
     - :math:`T_R`
   * - *TL*
     - :math:`T_L`
   * - *TRdot*
     - :math:`\dot T_R`
   * - *TRddot*
     - :math:`\ddot T_R`
   * - *imTR*
     - :math:`\Im T_R`
   * - *psi_anal*
     - :math:`\psi^{\mathrm anal}`
   * - *TR_anal*
     - :math:`T_R^{\mathrm anal}`
   * - *TReff2m*
     - :math:`|T_R|^2_{\mathrm eff}`
   * - *TReff2km*
     - :math:`k|T_R|^2_{\mathrm eff}`
   * - *TRdoteff2m*
     - :math:`|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TRdoteff2km*
     - :math:`k|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TLeff2m*
     - :math:`|T_R|^2_{\mathrm eff}`
   * - *TLeff2km*
     - :math:`k|T_R|^2_{\mathrm eff}`
   * - *TLdoteff2m*
     - :math:`|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TLdoteff2km*
     - :math:`k|T_R\dot{T}_R|_{\mathrm eff}`
   * - *dgrant_up*
     - :math:`{\cal T}^\chi`
   * - *grand2*
     - :math:`{\cal T}^Q`  (test)
   * - *dgrant*
     - :math:`\dot{\cal T}^\chi`
   * - *fact*
     - :math:`\Theta(t)`
   * - *k0*
     - :math:`k0`
   * - *dk*
     - :math:`dk`

Module *obsolete/visc_smagorinsky.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *nu_LES*
     - Mean value of Smagorinsky viscosity

Module *special/1D_loop.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtchi2*
     - heatconduction
   * - *dtrad*
     - radiative loss from RTV
   * - *dtspitzer*
     - Spitzer heat conduction time step
   * - *qmax*
     - max of heat flux vector
   * - *qrms*
     - rms of heat flux vector

Module *special/advective_gauge.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *LamRAm*
     - :math:`\left<\Lambda_{r\to a}\right>`
   * - *LamRApt*
     - :math:`\Lambda_{r\to a}(x1,y1,z1)`
   * - *LamRAp2*
     - :math:`\Lambda_{r\to a}(x2,y2,z2)`
   * - *LamRArms*
     - :math:`\left<\Lambda_{r\to a}^2\right>^{1/2}`
   * - *LamRAbzm*
     - :math:`\left<\Lambda_{r\to a} B_z\right>`
   * - *LamRAbzmz*
     - :math:`\left<\Lambda_{r\to a} B_z\right>_{xy}`
   * - *gLamRAbm*
     - :math:`\left<\Lambda_{r\to a}\Bv\right>`
   * - *apbrms*
     - :math:`\left<(\Av'\Bv)^2\right>^{1/2}`
   * - *jxarms*
     - :math:`\left<(\Jv\times\Av)^2\right>^{1/2}`
   * - *jxaprms*
     - :math:`\left<(\Jv\times\Av')^2\right>^{1/2}`
   * - *jxgLamRArms*
     - :math:`\left<(\Jv\times\nabla\Lambda_{r\to a})^2\right>^{1/2}`
   * - *gLamRArms*
     - :math:`\left<(\nabla\Lambda_{r\to a})^2\right>^{1/2}`
   * - *divabrms*
     - :math:`\left<[(\nabla\cdot\Av)\Bv]^2\right>^{1/2}`
   * - *divapbrms*
     - :math:`\left<[(\nabla\cdot\Av')\Bv]^2\right>^{1/2}`
   * - *d2LamRAbrms*
     - :math:`\left<[(\nabla^2\Lambda_{r\to a})\Bv]^2\right>^{1/2}`
   * - *d2LamRArms*
     - :math:`\left<[\nabla^2\Lambda_{r\to a}]^2\right>^{1/2}`
   * - *aabm*
     - :math:`\left<\Av_\mathrm{Adv}\cdot\Bv\right>`

Module *special/ascale_collapse.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *redshift*
     - redshift :math:`z`
   * - *Hubble*
     - :math:`H(a)`
   * - *ascale*
     - :math:`a`
   * - *lna*
     - :math:`\ln a`
   * - *tph*
     - :math:`t_\mathrm{phys}`

Module *special/axionSU2back.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *a*
     - :math:`a`
   * - *phi*
     - :math:`phi`
   * - *phidot*
     - :math:`phidot`
   * - *H*
     - :math:`Hubble_parameter`
   * - *Q*
     - :math:`Q`
   * - *Qdot*
     - :math:`\dot{Q}`
   * - *Qddot*
     - :math:`\ddot{Q}`
   * - *chi*
     - :math:`\chi`
   * - *chidot*
     - :math:`\dot{\chi}`
   * - *chiddot*
     - :math:`\ddot{\chi}`
   * - *psi*
     - :math:`\psi`
   * - *psiL*
     - :math:`\psi_L`
   * - *psidot*
     - :math:`\dot\psi`
   * - *psiddot*
     - :math:`\ddot\psi`
   * - *TR*
     - :math:`T_R`
   * - *TL*
     - :math:`T_L`
   * - *uR*
     - :math:`u_R`
   * - *uL*
     - :math:`u_L`
   * - *TRdot*
     - :math:`\dot T_R`
   * - *TRddot*
     - :math:`\ddot T_R`
   * - *imTR*
     - :math:`\Im T_R`
   * - *psi_anal*
     - :math:`\psi^{\mathrm anal}`
   * - *TR_anal*
     - :math:`T_R^{\mathrm anal}`
   * - *TReff2m*
     - :math:`|T_R|^2_{\mathrm eff}`
   * - *TReff2km*
     - :math:`k|T_R|^2_{\mathrm eff}`
   * - *TRdoteff2m*
     - :math:`|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TRdoteff2km*
     - :math:`k|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TRpsim*
     - :math:`\langle T_R^* \psi\rangle`
   * - *TRpsikm*
     - :math:`\langle T_R^* \psi (k/a)\rangle`
   * - *TRpsidotm*
     - :math:`\langle T_R^* \psi'\rangle`
   * - *TRdotpsim*
     - :math:`\langle {T_R^*}' \psi\rangle`
   * - *TLeff2m*
     - :math:`|T_L|^2_{\mathrm eff}`
   * - *TLeff2km*
     - :math:`k|T_L|^2_{\mathrm eff}`
   * - *uReff2m*
     - :math:`|u_R|^2_{\mathrm eff}`
   * - *uReff2km*
     - :math:`k|u_R|^2_{\mathrm eff}`
   * - *uLeff2m*
     - :math:`|u_L|^2_{\mathrm eff}`
   * - *uLeff2km*
     - :math:`k|u_L|^2_{\mathrm eff}`
   * - *TLdoteff2m*
     - :math:`|T_R\dot{T}_R|_{\mathrm eff}`
   * - *TLdoteff2km*
     - :math:`k|T_R\dot{T}_R|_{\mathrm eff}`
   * - *dgrant_up*
     - :math:`{\cal T}^\chi`
   * - *grand2*
     - :math:`{\cal T}^Q`  (test)
   * - *dgrant*
     - :math:`\dot{\cal T}^\chi`
   * - *rhoT*
     - :math:`\rho_T`
   * - *JJ_R*
     - :math:`J_R`
   * - *JJ_L*
     - :math:`J_L`
   * - *JJ*
     - :math:`J`
   * - *fact*
     - :math:`\Theta(t)`
   * - *k0*
     - :math:`k0`
   * - *dk*
     - :math:`dk`

Module *special/axionU1back.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *a*
     - :math:`a`
   * - *phi*
     - :math:`phi`
   * - *phidot*
     - :math:`phidot`
   * - *H*
     - :math:`Hubble_parameter`
   * - *AL*
     - :math:`\dot\psi`
   * - *ALdot*
     - :math:`\ddot\psi`
   * - *AR*
     - :math:`T_R`
   * - *ARdot*
     - :math:`T_L`
   * - *imAR*
     - :math:`\Im T_R`
   * - *imAL*
     - :math:`T_R^{\mathrm anal}`
   * - *imARdot*
     - :math:`\Im T_R`
   * - *k0*
     - :math:`k0`
   * - *dk*
     - :math:`dk`
   * - *rhoe*
     - :math:`rho_e`
   * - *rhob*
     - :math:`rho_B`
   * - *edotb*
     - :math:`edotb`

Module *special/backreact_infl.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *phim*
     - :math:`\left<\phi\right>`
   * - *phi2m*
     - :math:`\left<\phi^2\right>`
   * - *phirms*
     - :math:`\left<\phi^2\right>^{1/2}`
   * - *dphim*
     - :math:`\left<\phi'\right>`
   * - *dphi2m*
     - :math:`\left<(\phi')^2\right>`
   * - *dphirms*
     - :math:`\left<(\phi')^2\right>^{1/2}`
   * - *dtphi*
     - :math:`dt/cdtphi`
   * - *Hscriptm*
     - :math:`\left<{\cal a*H}\right>`
   * - *ascale*
     - :math:`a`
   * - *lnam*
     - :math:`\left<\ln a\right>`
   * - *ddotam*
     - :math:`a''/a`
   * - *a2rhopm*
     - :math:`a^2 (\rho+p)`
   * - *a2rhom*
     - :math:`a^2 \rho`
   * - *a2rhophim*
     - :math:`a^2 \rho_\phi`
   * - *a4rhophim*
     - :math:`a^4 \rho_\phi`
   * - *a2rhogphim*
     - :math:`0.5 <grad \phi^2>`
   * - *rho_chi*
     - :math:`\rho_\chi`
   * - *rho_rad*
     - :math:`\rho_\mathrm{rad}`
   * - *sigEma*
     - :math:`\rho_\chi`
   * - *sigBma*
     - :math:`\rho_\chi`
   * - *count_eb0a*
     - :math:`f_\mathrm{EB0}`
   * - *heating*
     - :math:`\theta_\mathrm{heat}`
   * - *wstate*
     - :math:`w_\mathrm{state}`
   * - *wstate_aver*
     - :math:`\langle w_\mathrm{state}\rangle`
   * - *Gamma_phi*
     - :math:`\langle w_\mathrm{state}\rangle`
   * - *Gam_phi*
     - :math:`\langle w_\mathrm{state}\rangle`

Module *special/chiral_mhd.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *muSm*
     - :math:`\left<\mu_S\right>`
   * - *muSrms*
     - :math:`\left<\mu_S^2\right>^{1/2}`
   * - *muSmax*
     - :math:`\max\mu`
   * - *mu5m*
     - :math:`\left<\mu_5\right>`
   * - *mu51m*
     - :math:`\left<|\mu_5|\right>`
   * - *mu53m*
     - :math:`\left<\mu_5^3\right>`
   * - *mu54m*
     - :math:`\left<\mu_5^4\right>`
   * - *mu5rms*
     - :math:`\left<\mu_5^2\right>^{1/2}`
   * - *mu5min*
     - :math:`\min\mu_5`
   * - *mu5max*
     - :math:`\max\mu_5`
   * - *mu5abs*
     - :math:`\max|\mu_5|`
   * - *srce5m*
     - :math:`\left<S_5\right>`
   * - *gamf5m*
     - :math:`\left<\Gamma_5\right>`
   * - *gmu5rms*
     - :math:`\left<(\nabla\mu_5)^2\right>^{1/2}`
   * - *gmuSrms*
     - :math:`\left<(\nabla\mu_S)^2\right>^{1/2}`
   * - *gmu5mx*
     - :math:`\left<\nabla\mu_5\right>_x`
   * - *gmu5my*
     - :math:`\left<\nabla\mu_5\right>_y`
   * - *gmu5mz*
     - :math:`\left<\nabla\mu_5\right>_z`
   * - *bgmu5rms*
     - :math:`\left<(\Bv\cdot\nabla\mu_5)^2\right>^{1/2}`
   * - *bgmuSrms*
     - :math:`\left<(\Bv\cdot\nabla\mu_S)^2\right>^{1/2}`
   * - *mu5bjm*
     - :math:`\left<\mu_5 ((\nabla\times\Bv)\cdot\Bv) \right>`
   * - *mu5bjrms*
     - :math:`\left<(\mu_5 ((\nabla\times\Bv)\cdot\Bv))^2 \right>^{1/2}`
   * - *dt_lambda5*
     - :math:`\mathrm{min}(\mu_5/\Bv^2) \delta x/(\lambda \eta)`
   * - *dt_D5*
     - :math:`(\lambda \eta \mathrm{max}(\Bv^2))^{-1}`
   * - *dt_gammaf5*
     - :math:`1/\Gamma_\mathrm{f}`
   * - *dt_CMW*
     - :math:`\delta x/((C_\mu C_5)^{1/2} \mathrm{max}(|\Bv|))`
   * - *dt_Dmu*
     - :math:`(\lambda \eta \mathrm{min}(\Bv^2))^{-1}`
   * - *dt_vmu*
     - :math:`\delta x /(\eta \mathrm{max}(|\mu_5 |))`
   * - *dt_chiral*
     - total time-step contribution from chiral MHD
   * - *mu5bxm*
     - :math:`\left<\mu_5B_x\right>`
   * - *mu5b2m*
     - :math:`\left<\mu_5B^2\right>`
   * - *mu5jbm*
     - :math:`\left<\mu_5\Jv\cdot\Bv\right>`
   * - *jxm*
     - :math:`\langle J_x\rangle`
   * - *Dmu5_tdep*
     - :math:`D(t)`

Module *special/collapse.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *betm*
     - :math:`\left<\beta\right>`
   * - *massm*
     - :math:`\left<m\right>`

Module *special/coronae.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtchi2*
     - :math:`\delta t / [c_{\delta t,{\mathrm v}}\, \delta x^2/\chi_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat conductivity; see \S~\ref{time-step})
   * - *dtspitzer*
     - Spitzer heat conduction time step
   * - *dtrad*
     - radiative loss from RTV

Module *special/disp_current.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *EEEM*
     - :math:`\left<\Ev^2+\Bv^2\right>/2`
   * - *erms*
     - :math:`\left<\Ev^2\right>^{1/2}`
   * - *eprimerms*
     - :math:`\left<(E')^2\right>^{1/2}`
   * - *bprimerms*
     - :math:`\left<(B')^2\right>^{1/2}`
   * - *jprimerms*
     - :math:`\left<(J')^2\right>^{1/2}`
   * - *gam_EBrms*
     - :math:`\left<(\gamma')^2\right>^{1/2}`
   * - *boostprms*
     - :math:`\left<\mbox{boost}^2\right>^{1/2}`
   * - *edotrms*
     - :math:`\left<\dot{\Ev}^2\right>^{1/2}`
   * - *emax*
     - :math:`\max(|\Ev|)`
   * - *a0rms*
     - :math:`\left<A_0^2\right>^{1/2}`
   * - *grms*
     - :math:`\left<C-\nabla\cdot\Av\right>^{1/2}`
   * - *da0rms*
     - :math:`\left<C-\nabla\cdot\Av\right>^{1/2}`
   * - *BcurlEm*
     - :math:`\left<\Bv\cdot\nabla\times\Ev\right>`
   * - *divJrms*
     - :math:`\left<\nab\Jv^2\right>^{1/2}`
   * - *divErms*
     - :math:`\left<\nab\Ev^2\right>^{1/2}`
   * - *rhoerms*
     - :math:`\left<\rho_e^2\right>^{1/2}`
   * - *divJm*
     - :math:`\left<\nab\Jv\right>`
   * - *divEm*
     - :math:`\left<\nab\Ev\right>`
   * - *rhoem*
     - :math:`\left<\rho_e\right>`
   * - *count_eb0*
     - :math:`f_\mathrm{EB0}`
   * - *mfpf*
     - :math:`-f'/f`
   * - *fppf*
     - :math:`f''/f`
   * - *afact*
     - :math:`a`  (scale factor)
   * - *constrainteqn*
     - :math:`<deldotE+>`
   * - *exm*
     - :math:`\left<E_x\right>`
   * - *eym*
     - :math:`\left<E_y\right>`
   * - *ezm*
     - :math:`\left<E_z\right>`
   * - *etaSchw*
     - :math:`1/\left<\sigma_\mathrm{E}\right>`
   * - *dteta*
     - :math:`dt/cdtv`
   * - *dtsigE*
     - :math:`dt/cdt\_sigE`
   * - *sigEm*
     - :math:`\left<\sigma_\mathrm{E}\right>`
   * - *sigBm*
     - :math:`\left<\sigma_\mathrm{B}\right>`
   * - *sigErms*
     - :math:`\left<\sigma_\mathrm{E}^2\right>^{1/2}`
   * - *sigBrms*
     - :math:`\left<\sigma_\mathrm{B}^2\right>^{1/2}`
   * - *sigEE2m*
     - :math:`\left<\sigma_\mathrm{E}\Ev^2\right>`
   * - *sigBBEm*
     - :math:`\left<\sigma_\mathrm{E}\Bv\cdot\Ev\right>`
   * - *adphiBm*
     - :math:`\left<(\alpha/f)<\phi'\Bv\cdot\Ev\right>`
   * - *adphiBrms*
     - :math:`\left<[(\alpha/f)<\phi'\Bv]\right>^{1/2}`
   * - *Johmrms*
     - :math:`\left<\Jv^2\right>^{1/2}`
   * - *J2sigEm*
     - :math:`\left<\Jv^2/\sigma_E\right>`
   * - *ujxb1m*
     - :math:`\left<\uv\cdot(\Jv\times\Bv)_E\right>`
   * - *curlBrms*
     - :math:`\left<(\mathrm{curl}\Bv)^2\right>^{1/2}`
   * - *echarge*
     - :math:`\left<e_\mathrm{eff}\right>`
   * - *ebm*
     - :math:`\left<\Ev\cdot\Bv\right>`

Module *special/electroweaksu2.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *W1rms*
     - :math:`\left<{\Wv^1}^2\right>^{1/2}`
   * - *W2rms*
     - :math:`\left<{\Wv^2}^2\right>^{1/2}`
   * - *W3rms*
     - :math:`\left<{\Wv^3}^2\right>^{1/2}`
   * - *W1max*
     - :math:`\max(|\Wv^1|)`
   * - *W2max*
     - :math:`\max(|\Wv^2|)`
   * - *W3max*
     - :math:`\max(|\Wv^3|)`
   * - *dW1rms*
     - :math:`\left<{{\dot{\Wv}}_1}^2\right>^{1/2}`
   * - *dW2rms*
     - :math:`\left<{{\dot{\Wv}}_2}^2\right>^{1/2}`
   * - *dW3rms*
     - :math:`\left<{{\dot{\Wv}}_3}^2\right>^{1/2}`
   * - *dW1max*
     - :math:`\max(|\dot{\Wv}_1|)`
   * - *dW2max*
     - :math:`\max(|\dot{\Wv}_2|)`
   * - *dW3max*
     - :math:`\max(|\dot{\Wv}_3|)`
   * - *W1ddotrms*
     - :math:`\left<{\ddot{\Wv}_1}^2\right>^{1/2}`
   * - *W2ddotrms*
     - :math:`\left<{\ddot{\Wv}_2}^2\right>^{1/2}`
   * - *W3ddotrms*
     - :math:`\left<{\ddot{\Wv}_3}^2\right>^{1/2}`
   * - *divW1rms*
     - :math:`\left<{\nab{\Wv}_1}^2\right>^{1/2}`
   * - *divW2rms*
     - :math:`\left<{\nab{\Wv}_2}^2\right>^{1/2}`
   * - *divW3rms*
     - :math:`\left<{\nab{\Wv}_3}^2\right>^{1/2}`
   * - *divW1m*
     - :math:`\left<\nab{\Wv}_1\right>`
   * - *divW2m*
     - :math:`\left<\nab{\Wv}_2\right>`
   * - *divW3m*
     - :math:`\left<\nab{\Wv}_3\right>`
   * - *divdotW1rms*
     - :math:`\left<{\nab\dot{\Wv}_1}^2\right>^{1/2}`
   * - *divdotW2rms*
     - :math:`\left<{\nab\dot{\Wv}_2}^2\right>^{1/2}`
   * - *divdotW3rms*
     - :math:`\left<{\nab\dot{\Wv}_3}^2\right>^{1/2}`
   * - *rhoW1rms*
     - :math:`\left<{\rho_{\Wv_1}}^2\right>^{1/2}`
   * - *rhoW2rms*
     - :math:`\left<{\rho_{\Wv_2}}^2\right>^{1/2}`
   * - *rhoW3rms*
     - :math:`\left<{\rho_{\Wv_3}}^2\right>^{1/2}`
   * - *divdotW1m*
     - :math:`\left<\nab\dot\Wv_1\right>`
   * - *divdotW2m*
     - :math:`\left<\nab\dot\Wv_2\right>`
   * - *divdotW3m*
     - :math:`\left<\nab\dot\Wv_3\right>`
   * - *rhoW1m*
     - :math:`\left<\rho_e\Wv_1\right>`
   * - *rhoW2m*
     - :math:`\left<\rho_e\Wv_2\right>`
   * - *rhoW3m*
     - :math:`\left<\rho_e\Wv_3\right>`
   * - *constrainteqnW*
     - :math:`<deldotW+>`
   * - *W1xm*
     - :math:`\left<W_x^1\right>`
   * - *W1ym*
     - :math:`\left<W_y^1\right>`
   * - *W1zm*
     - :math:`\left<W_z^1\right>`
   * - *W2xm*
     - :math:`\left<W_x^2\right>`
   * - *W2ym*
     - :math:`\left<W_y^2\right>`
   * - *W2zm*
     - :math:`\left<W_z^2\right>`
   * - *W3xm*
     - :math:`\left<W_x^3\right>`
   * - *W3ym*
     - :math:`\left<W_y^3\right>`
   * - *W3zm*
     - :math:`\left<W_z^3\right>`
   * - *dW1xm*
     - :math:`\left<\dot{W}^1_x\right>`
   * - *dW1ym*
     - :math:`\left<\dot{W}^1_y\right>`
   * - *dW1zm*
     - :math:`\left<\dot{W}^1_z\right>`
   * - *dW2xm*
     - :math:`\left<\dot{W}^2_x\right>`
   * - *dW2ym*
     - :math:`\left<\dot{W}^2_y\right>`
   * - *dW2zm*
     - :math:`\left<\dot{W}^2_z\right>`
   * - *dW3xm*
     - :math:`\left<\dot{W}^3_x\right>`
   * - *dW3ym*
     - :math:`\left<\dot{W}^3_y\right>`
   * - *dW3zm*
     - :math:`\left<\dot{W}^3_z\right>`
   * - *W1dotW1m*
     - :math:`\left<\dot{\Wv}_1\cdot{\Wv}_1\right>`
   * - *W2dotW2m*
     - :math:`\left<\dot{\Wv}_2\cdot{\Wv}_2\right>`
   * - *W3dotW3m*
     - :math:`\left<\dot{\Wv}_3\cdot{\Wv}_3\right>`

Module *special/gravitational_waves.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *hhT2m*
     - :math:`\left<h_{\mathrm T}^2\right>`
   * - *hhX2m*
     - :math:`\left<h_{\mathrm X}^2\right>`
   * - *hhThhXm*
     - :math:`\left<h_{\mathrm T}h_{\mathrm X}\right>`
   * - *ggTpt*
     - :math:`g_{\mathrm T}(x_1,y_1,z_1,t)`
   * - *strTpt*
     - :math:`S_{\mathrm T}(x_1,y_1,z_1,t)`
   * - *strXpt*
     - :math:`S_{\mathrm X}(x_1,y_1,z_1,t)`

Module *special/gravitational_waves_hij6.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *h22rms*
     - :math:`h_{22}^{\mathrm rms}`
   * - *h33rms*
     - :math:`h_{33}^{\mathrm rms}`
   * - *h23rms*
     - :math:`h_{23}^{\mathrm rms}`
   * - *g11pt*
     - :math:`g_{11}(x_1,y_1,z_1,t)`
   * - *g22pt*
     - :math:`g_{22}(x_1,y_1,z_1,t)`
   * - *g33pt*
     - :math:`g_{33}(x_1,y_1,z_1,t)`
   * - *g12pt*
     - :math:`g_{12}(x_1,y_1,z_1,t)`
   * - *g23pt*
     - :math:`g_{23}(x_1,y_1,z_1,t)`
   * - *g31pt*
     - :math:`g_{31}(x_1,y_1,z_1,t)`
   * - *hhTpt*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXpt*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTpt*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXpt*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hhTp2*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXp2*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTp2*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXp2*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hrms*
     - :math:`\bra{h_T^2+h_X^2}^{1/2}`
   * - *EEGW*
     - :math:`\bra{g_T^2+g_X^2}\,c^2/(32\pi G)`
   * - *gg2m*
     - :math:`\bra{g_T^2+g_X^2}`
   * - *hhT2m*
     - :math:`\bra{h_T^2}`
   * - *hhX2m*
     - :math:`\bra{h_X^2}`
   * - *hhTXm*
     - :math:`\bra{h_T h_X}`
   * - *ggT2m*
     - :math:`\bra{g_T^2}`
   * - *ggX2m*
     - :math:`\bra{g_X^2}`
   * - *ggTXm*
     - :math:`\bra{g_T g_X}`
   * - *ggTm*
     - :math:`\bra{g_T}`
   * - *ggXm*
     - :math:`\bra{g_X}`
   * - *hijij2m*
     - :math:`\bra{h_{ij,ij}^2}`
   * - *gijij2m*
     - :math:`\bra{g_{ij,ij}^2}`

Module *special/gravitational_waves_hTXk.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *STrept*
     - :math:`Re S_{T}(k_1,k_1,k_1,t)`
   * - *STimpt*
     - :math:`Im S_{T}(k_1,k_1,k_1,t)`
   * - *SXrept*
     - :math:`Re S_{X}(k_1,k_1,k_1,t)`
   * - *SXimpt*
     - :math:`Im S_{X}(k_1,k_1,k_1,t)`
   * - *hTrept*
     - :math:`Re h_{T}(k_1,k_1,k_1,t)`
   * - *hTimpt*
     - :math:`Im h_{T}(k_1,k_1,k_1,t)`
   * - *hXrept*
     - :math:`Re h_{X}(k_1,k_1,k_1,t)`
   * - *hXimpt*
     - :math:`Im h_{X}(k_1,k_1,k_1,t)`
   * - *gTrept*
     - :math:`Re h_{T}(k_1,k_1,k_1,t)`
   * - *gTimpt*
     - :math:`Im h_{T}(k_1,k_1,k_1,t)`
   * - *gXrept*
     - :math:`Re h_{X}(k_1,k_1,k_1,t)`
   * - *gXimpt*
     - :math:`Im h_{X}(k_1,k_1,k_1,t)`
   * - *STrep2*
     - :math:`Re S_{T}(k_2,k_2,k_2,t)`
   * - *STimp2*
     - :math:`Im S_{T}(k_2,k_2,k_2,t)`
   * - *SXrep2*
     - :math:`Re S_{X}(k_2,k_2,k_2,t)`
   * - *SXimp2*
     - :math:`Im S_{X}(k_2,k_2,k_2,t)`
   * - *hTrep2*
     - :math:`Re h_{T}(k_2,k_2,k_2,t)`
   * - *hTimp2*
     - :math:`Im h_{T}(k_2,k_2,k_2,t)`
   * - *hXrep2*
     - :math:`Re h_{X}(k_2,k_2,k_2,t)`
   * - *hXimp2*
     - :math:`Im h_{X}(k_2,k_2,k_2,t)`
   * - *gTrep2*
     - :math:`Re g_{T}(k_2,k_2,k_2,t)`
   * - *gTimp2*
     - :math:`Im g_{T}(k_2,k_2,k_2,t)`
   * - *gXrep2*
     - :math:`Re g_{X}(k_2,k_2,k_2,t)`
   * - *gXimp2*
     - :math:`Im g_{X}(k_2,k_2,k_2,t)`
   * - *g11pt*
     - :math:`g_{11}(x_1,y_1,z_1,t)`
   * - *g22pt*
     - :math:`g_{22}(x_1,y_1,z_1,t)`
   * - *g33pt*
     - :math:`g_{33}(x_1,y_1,z_1,t)`
   * - *g12pt*
     - :math:`g_{12}(x_1,y_1,z_1,t)`
   * - *g23pt*
     - :math:`g_{23}(x_1,y_1,z_1,t)`
   * - *g31pt*
     - :math:`g_{31}(x_1,y_1,z_1,t)`
   * - *hhTpt*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXpt*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTpt*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXpt*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hhTp2*
     - :math:`h_{T}(x_1,y_1,z_1,t)`
   * - *hhXp2*
     - :math:`h_{X}(x_1,y_1,z_1,t)`
   * - *ggTp2*
     - :math:`\dot{h}_{T}(x_1,y_1,z_1,t)`
   * - *ggXp2*
     - :math:`\dot{h}_{X}(x_1,y_1,z_1,t)`
   * - *hrms*
     - :math:`\bra{h_T^2+h_X^2}^{1/2}`
   * - *EEGW*
     - :math:`\bra{g_T^2+g_X^2}\,c^2/(32\pi G)`
   * - *gg2m*
     - :math:`\bra{g_T^2+g_X^2}`
   * - *Stgm*
     - :math:`\bra{S_Tg_T+S_Xg_X}`
   * - *hhT2m*
     - :math:`\bra{h_T^2}`
   * - *hhX2m*
     - :math:`\bra{h_X^2}`
   * - *hhTXm*
     - :math:`\bra{h_T h_X}`
   * - *ggT2m*
     - :math:`\bra{g_T^2}`
   * - *ggX2m*
     - :math:`\bra{g_X^2}`
   * - *ggTXm*
     - :math:`\bra{g_T g_X}`
   * - *nlin0*
     - :math:`\bra{nlin0}`
   * - *nlin1*
     - :math:`\bra{nlin1}`
   * - *nlin2*
     - :math:`\bra{nlin2}`
   * - *h11rms*
     - :math:`\bra{h_{11}^2}^{1/2}`
   * - *h22rms*
     - :math:`\bra{h_{22}^2}^{1/2}`
   * - *h33rms*
     - :math:`\bra{h_{33}^2}^{1/2}`
   * - *h12rms*
     - :math:`\bra{h_{12}^2}^{1/2}`
   * - *h23rms*
     - :math:`\bra{h_{23}^2}^{1/2}`
   * - *h31rms*
     - :math:`\bra{h_{31}^2}^{1/2}`

Module *special/klein_gordon.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *phim*
     - :math:`\left<\phi\right>`
   * - *phi2m*
     - :math:`\left<\phi^2\right>`
   * - *phirms*
     - :math:`\left<\phi^2\right>^{1/2}`
   * - *dphim*
     - :math:`\left<\phi'\right>`
   * - *dphi2m*
     - :math:`\left<(\phi')^2\right>`
   * - *dphirms*
     - :math:`\left<(\phi')^2\right>^{1/2}`
   * - *psim*
     - :math:`\left<\psi\right>`
   * - *psi2m*
     - :math:`\left<\psi^2\right>`
   * - *psirms*
     - :math:`\left<\psi^2\right>^{1/2}`
   * - *dpsim*
     - :math:`\left<\psi'\right>`
   * - *dpsi2m*
     - :math:`\left<(\psi')^2\right>`
   * - *dpsirms*
     - :math:`\left<(\psi')^2\right>^{1/2}`
   * - *Hscriptm*
     - :math:`\left<{\cal a*H}\right>`
   * - *lnam*
     - :math:`\left<\ln a\right>`
   * - *Vprimem*
     - :math:`\left<V_{,\phi}\right>`
   * - *Vprimepsim*
     - :math:`\left<V_{,\psi}\right>`
   * - *ddotam*
     - :math:`a''/a`
   * - *a2rhopm*
     - :math:`a^2 (rho+p)`
   * - *a2rhom*
     - :math:`a^2 rho`
   * - *a2rhophim*
     - :math:`a^2 rho`
   * - *a2rhogphim*
     - :math:`0.5 <grad phi^2>`
   * - *a2rhopsim*
     - :math:`a^2 rho`
   * - *a2rhogpsim*
     - :math:`0.5 <grad psi^2>`
   * - *rho_chi*
     - :math:`\rho_\chi`
   * - *sigEma*
     - :math:`\rho_\chi`
   * - *sigBma*
     - :math:`\rho_\chi`
   * - *count_eb0a*
     - :math:`f_\mathrm{EB0}`

Module *special/Lambda_CDM.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *redshift*
     - redshift :math:`z`
   * - *Hubble*
     - :math:`H(a)`
   * - *ascale*
     - :math:`a`
   * - *lna*
     - :math:`\ln a`
   * - *tph*
     - :math:`t_\mathrm{phys}`

Module *special/lorenz_gauge.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *phim*
     - :math:`\left<\phi\right>`
   * - *phipt*
     - :math:`\phi(x1,y1,z1)`
   * - *phip2*
     - :math:`\phi(x2,y2,z2)`
   * - *phibzm*
     - :math:`\left<\phi B_z\right>`
   * - *phibzmz*
     - :math:`\left<\phi B_z\right>_{xy}`

Module *special/lucky_droplet.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rad*
     - :math:`r/r_\ast`
   * - *tauk*
     - :math:`\tau_k`
   * - *tt1m*
     - :math:`\langle T \rangle`
   * - *qq1m*
     - :math:`\langle \ln T \rangle`
   * - *qq2m*
     - :math:`\langle \ln T^2 \rangle`
   * - *qq3m*
     - :math:`\langle \ln T^3 \rangle`
   * - *qq4m*
     - :math:`\langle \ln T^4 \rangle`

Module *special/radial_dist_func.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *rad*
     - :math:`r/r_\ast`
   * - *tauk*
     - :math:`\tau_k`
   * - *tt1m*
     - :math:`\langle T \rangle`
   * - *qq1m*
     - :math:`\langle \ln T \rangle`
   * - *qq2m*
     - :math:`\langle \ln T^2 \rangle`
   * - *qq3m*
     - :math:`\langle \ln T^3 \rangle`
   * - *qq4m*
     - :math:`\langle \ln T^4 \rangle`

Module *special/reaction_0D.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *eem*
     - :math:`\langle \eta \rangle`
   * - *ee1m*
     - :math:`\langle |\eta| \rangle`
   * - *ee2m*
     - :math:`\langle \eta^2 \rangle`
   * - *ee3m*
     - :math:`\langle \eta^3 \rangle`
   * - *ee4m*
     - :math:`\langle \eta^4 \rangle`
   * - *ee10*
     - :math:`\langle \eta_{10\%} \rangle`
   * - *ee50*
     - :math:`\langle \eta_{50\%} \rangle`
   * - *ee90*
     - :math:`\langle \eta_{90\%} \rangle`
   * - *ee99*
     - :math:`\langle \eta_{99\%} \rangle`
   * - *AAm*
     - :math:`\langle [A] \rangle`
   * - *DDm*
     - :math:`\langle [D] \rangle`
   * - *LLm*
     - :math:`\langle [L] \rangle`
   * - *DLm*
     - :math:`\langle [D]+[L] \rangle`
   * - *kC*
     - :math:`k_C`
   * - *A1*
     - :math:`A_1`
   * - *A2*
     - :math:`A_2`
   * - *A3*
     - :math:`A_3`
   * - *A4*
     - :math:`A_4`
   * - *A5*
     - :math:`A_5`
   * - *D1*
     - :math:`D_1`
   * - *D2*
     - :math:`D_2`
   * - *D3*
     - :math:`D_3`
   * - *D4*
     - :math:`D_4`
   * - *D5*
     - :math:`D_5`
   * - *L1*
     - :math:`L_1`
   * - *L2*
     - :math:`L_2`
   * - *L3*
     - :math:`L_3`
   * - *L4*
     - :math:`L_4`
   * - *L5*
     - :math:`L_5`

Module *special/rel_1d.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *betm*
     - :math:`\left<\beta\right>`
   * - *betmax*
     - :math:`\beta_{\max}`

Module *special/schur_triple_decomp.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *uSH2*
     - :math:`\left<{u^\mathrm{SH}}^2\right>`
   * - *uRR2*
     - :math:`\left<{u^\mathrm{RR}}^2\right>`
   * - *uEL2*
     - :math:`\left<{u^\mathrm{EL}}^2\right>`
   * - *uRRm*
     - :math:`\left<2u^\mathrm{SH}u^\mathrm{RR}\right>`
   * - *bSH2*
     - :math:`\left<{b^\mathrm{SH}}^2\right>`
   * - *bRR2*
     - :math:`\left<{b^\mathrm{RR}}^2\right>`
   * - *bEL2*
     - :math:`\left<{b^\mathrm{EL}}^2\right>`
   * - *bRRm*
     - :math:`\left<2b^\mathrm{SH}b^\mathrm{RR}\right>`

Module *special/solar_corona.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *dtvel*
     - Velocity driver time step
   * - *dtnewt*
     - Radiative cooling time step
   * - *dtradloss*
     - Radiative losses time step
   * - *dtchi2*
     - :math:`\delta t / [c_{\delta t,{\mathrm v}}\, \delta x^2/\chi_{\mathrm max}]`  :math:`\quad`  (time step relative to time step based on heat conductivity; see \S~\ref{time-step})
   * - *dtspitzer*
     - Spitzer heat conduction time step
   * - *mag_flux*
     - Total vertical magnetic flux at

Module *test_methods/testfield_axisym4.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alpPERP*
     - :math:`\alpha_\perp`
   * - *alpPARA*
     - :math:`\alpha_\perp`
   * - *gam*
     - :math:`\gamma`
   * - *betPERP*
     - :math:`\beta_\perp`
   * - *betPERP2*
     - :math:`\beta_\perp^{(2)}`
   * - *betPARA*
     - :math:`\beta_\perp`
   * - *del*
     - :math:`\delta`
   * - *del2*
     - :math:`\delta^{(2)}`
   * - *kapPERP*
     - :math:`\kappa_\perp`
   * - *kapPERP2*
     - :math:`\kappa_\perp^{(2)}`
   * - *kapPARA*
     - :math:`\kappa_\perp`
   * - *mu*
     - :math:`\mu`
   * - *mu2*
     - :math:`\mu^{(2)}`
   * - *alpPERPz*
     - :math:`\alpha_\perp(z)`
   * - *alpPARAz*
     - :math:`\alpha_\perp(z)`
   * - *gamz*
     - :math:`\gamma(z)`
   * - *betPERPz*
     - :math:`\beta_\perp(z)`
   * - *betPARAz*
     - :math:`\beta_\perp(z)`
   * - *delz*
     - :math:`\delta(z)`
   * - *kapPERPz*
     - :math:`\kappa_\perp(z)`
   * - *kapPARAz*
     - :math:`\kappa_\perp(z)`
   * - *muz*
     - :math:`\mu(z)`
   * - *bx1pt*
     - :math:`b_x^{1}`
   * - *bx2pt*
     - :math:`b_x^{2}`
   * - *bx3pt*
     - :math:`b_x^{3}`
   * - *b1rms*
     - :math:`\left<b_{1}^2\right>^{1/2}`
   * - *b2rms*
     - :math:`\left<b_{2}^2\right>^{1/2}`
   * - *b3rms*
     - :math:`\left<b_{3}^2\right>^{1/2}`

Module *test_methods/testfield_compress_z.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *alp31*
     - :math:`\alpha_{31}`
   * - *alp12*
     - :math:`\alpha_{12}`
   * - *alp22*
     - :math:`\alpha_{22}`
   * - *alp32*
     - :math:`\alpha_{32}`
   * - *eta11*
     - :math:`\eta_{11}k`
   * - *eta21*
     - :math:`\eta_{21}k`
   * - *eta12*
     - :math:`\eta_{12}k`
   * - *eta22*
     - :math:`\eta_{22}k`
   * - *alpK*
     - :math:`\alpha^K`
   * - *alpM*
     - :math:`\alpha^M`
   * - *alpMK*
     - :math:`\alpha^{MK}`
   * - *phi11*
     - :math:`\phi_{11}`
   * - *phi21*
     - :math:`\phi_{21}`
   * - *phi12*
     - :math:`\phi_{12}`
   * - *phi22*
     - :math:`\phi_{22}`
   * - *phi32*
     - :math:`\phi_{32}`
   * - *psi11*
     - :math:`\psi_{11}k`
   * - *psi21*
     - :math:`\psi_{21}k`
   * - *psi12*
     - :math:`\psi_{12}k`
   * - *psi22*
     - :math:`\psi_{22}k`
   * - *sig1*
     - :math:`\sigma_1`
   * - *sig2*
     - :math:`\sigma_2`
   * - *sig3*
     - :math:`\sigma_3`
   * - *tau1*
     - :math:`\tau_1`
   * - *tau2*
     - :math:`\tau_2`
   * - *phiK*
     - :math:`\phi^K`
   * - *phiM*
     - :math:`\phi^M`
   * - *phiMK*
     - :math:`\phi^{MK}`
   * - *alp11cc*
     - :math:`\alpha_{11}\cos^2 kz`
   * - *alp21sc*
     - :math:`\alpha_{21}\sin kz\cos kz`
   * - *alp12cs*
     - :math:`\alpha_{12}\cos kz\sin kz`
   * - *alp22ss*
     - :math:`\alpha_{22}\sin^2 kz`
   * - *eta11cc*
     - :math:`\eta_{11}\cos^2 kz`
   * - *eta21sc*
     - :math:`\eta_{21}\sin kz\cos kz`
   * - *eta12cs*
     - :math:`\eta_{12}\cos kz\sin kz`
   * - *eta22ss*
     - :math:`\eta_{22}\sin^2 kz`
   * - *s2kzDFm*
     - :math:`\left<\sin2kz\nabla\cdot F\right>`
   * - *M11*
     - :math:`{\cal M}_{11}`
   * - *M22*
     - :math:`{\cal M}_{22}`
   * - *M33*
     - :math:`{\cal M}_{33}`
   * - *M11cc*
     - :math:`{\cal M}_{11}\cos^2 kz`
   * - *M11ss*
     - :math:`{\cal M}_{11}\sin^2 kz`
   * - *M22cc*
     - :math:`{\cal M}_{22}\cos^2 kz`
   * - *M22ss*
     - :math:`{\cal M}_{22}\sin^2 kz`
   * - *M12cs*
     - :math:`{\cal M}_{12}\cos kz\sin kz`
   * - *bx11pt*
     - :math:`b_x^{11}`
   * - *bx21pt*
     - :math:`b_x^{21}`
   * - *bx12pt*
     - :math:`b_x^{12}`
   * - *bx22pt*
     - :math:`b_x^{22}`
   * - *bx0pt*
     - :math:`b_x^{0}`
   * - *by11pt*
     - :math:`b_y^{11}`
   * - *by21pt*
     - :math:`b_y^{21}`
   * - *by12pt*
     - :math:`b_y^{12}`
   * - *by22pt*
     - :math:`b_y^{22}`
   * - *by0pt*
     - :math:`b_y^{0}`
   * - *u11rms*
     - :math:`\left<u_{11}^2\right>^{1/2}`
   * - *u21rms*
     - :math:`\left<u_{21}^2\right>^{1/2}`
   * - *u12rms*
     - :math:`\left<u_{12}^2\right>^{1/2}`
   * - *u22rms*
     - :math:`\left<u_{22}^2\right>^{1/2}`
   * - *h11rms*
     - :math:`\left<h_{11}^2\right>^{1/2}`
   * - *h21rms*
     - :math:`\left<h_{21}^2\right>^{1/2}`
   * - *h12rms*
     - :math:`\left<h_{12}^2\right>^{1/2}`
   * - *h22rms*
     - :math:`\left<h_{22}^2\right>^{1/2}`
   * - *j11rms*
     - :math:`\left<j_{11}^2\right>^{1/2}`
   * - *b11rms*
     - :math:`\left<b_{11}^2\right>^{1/2}`
   * - *b21rms*
     - :math:`\left<b_{21}^2\right>^{1/2}`
   * - *b12rms*
     - :math:`\left<b_{12}^2\right>^{1/2}`
   * - *b22rms*
     - :math:`\left<b_{22}^2\right>^{1/2}`
   * - *ux0m*
     - :math:`\left<u_{0_x}\right>`
   * - *uy0m*
     - :math:`\left<u_{0_y}\right>`
   * - *ux11m*
     - :math:`\left<u_{11_x}\right>`
   * - *uy11m*
     - :math:`\left<u_{11_y}\right>`
   * - *u0rms*
     - :math:`\left<u_{0}^2\right>^{1/2}`
   * - *b0rms*
     - :math:`\left<b_{0}^2\right>^{1/2}`
   * - *h0rms*
     - :math:`\left<h_{0}^2\right>^{1/2}`
   * - *rho0m*
     - :math:`\left<\exp h_{0}\right>`
   * - *u0max*
     - :math:`\operatorname{max}\left|\boldsymbol{u}_{0}\right|`
   * - *b0max*
     - :math:`\operatorname{max}\left|\boldsymbol{b}_{0}\right|`
   * - *h0max*
     - :math:`\operatorname{max}h_{0}`
   * - *bhrms*
     - :math:`\left<b_{h}^2\right>^{1/2}`
   * - *jb0m*
     - :math:`\left<jb_{0}\right>`
   * - *E11rms*
     - :math:`\left<{\cal E}_{11}^2\right>^{1/2}`
   * - *E21rms*
     - :math:`\left<{\cal E}_{21}^2\right>^{1/2}`
   * - *E12rms*
     - :math:`\left<{\cal E}_{12}^2\right>^{1/2}`
   * - *E22rms*
     - :math:`\left<{\cal E}_{22}^2\right>^{1/2}`
   * - *E0rms*
     - :math:`\left<{\cal E}_{0}^2\right>^{1/2}`
   * - *E0mrms*
     - :math:`\left<{\cal E}_{0}^2\right>^{1/2}`
   * - *E0xrms*
     - :math:`\left<{\cal E}_{0,x}^2\right>^{1/2}`
   * - *E0yrms*
     - :math:`\left<{\cal E}_{0,y}^2\right>^{1/2}`
   * - *Ex11pt*
     - :math:`{\cal E}_x^{11}`
   * - *Ex21pt*
     - :math:`{\cal E}_x^{21}`
   * - *Ex12pt*
     - :math:`{\cal E}_x^{12}`
   * - *Ex22pt*
     - :math:`{\cal E}_x^{22}`
   * - *Ex0pt*
     - :math:`{\cal E}_x^{0}`
   * - *Ey11pt*
     - :math:`{\cal E}_y^{11}`
   * - *Ey21pt*
     - :math:`{\cal E}_y^{21}`
   * - *Ey12pt*
     - :math:`{\cal E}_y^{12}`
   * - *Ey22pt*
     - :math:`{\cal E}_y^{22}`
   * - *Ey0pt*
     - :math:`{\cal E}_y^{0}`
   * - *bamp*
     - bamp
   * - *E111z*
     - :math:`{\cal E}_1^{11}`
   * - *E211z*
     - :math:`{\cal E}_2^{11}`
   * - *E311z*
     - :math:`{\cal E}_3^{11}`
   * - *E121z*
     - :math:`{\cal E}_1^{21}`
   * - *E221z*
     - :math:`{\cal E}_2^{21}`
   * - *E321z*
     - :math:`{\cal E}_3^{21}`
   * - *E112z*
     - :math:`{\cal E}_1^{12}`
   * - *E212z*
     - :math:`{\cal E}_2^{12}`
   * - *E312z*
     - :math:`{\cal E}_3^{12}`
   * - *E122z*
     - :math:`{\cal E}_1^{22}`
   * - *E222z*
     - :math:`{\cal E}_2^{22}`
   * - *E322z*
     - :math:`{\cal E}_3^{22}`
   * - *alp11z*
     - :math:`\alpha_{11}(z,t)`
   * - *alp21z*
     - :math:`\alpha_{21}(z,t)`
   * - *alp12z*
     - :math:`\alpha_{12}(z,t)`
   * - *alp22z*
     - :math:`\alpha_{22}(z,t)`
   * - *eta11z*
     - :math:`\eta_{11}(z,t)`
   * - *eta21z*
     - :math:`\eta_{21}(z,t)`
   * - *eta12z*
     - :math:`\eta_{12}(z,t)`
   * - *eta22z*
     - :math:`\eta_{22}(z,t)`
   * - *E10z*
     - :math:`{\cal E}_1^{0}`
   * - *E20z*
     - :math:`{\cal E}_2^{0}`
   * - *E30z*
     - :math:`{\cal E}_3^{0}`
   * - *EBpq*
     - :math:`{\cal E}\cdot\Bv^{pq}`
   * - *E0Um*
     - :math:`{\cal E}^0\cdot\Uv`
   * - *E0Wm*
     - :math:`{\cal E}^0\cdot\Wv`
   * - *bx0mz*
     - :math:`\left<b_{x}\right>_{xy}`
   * - *by0mz*
     - :math:`\left<b_{y}\right>_{xy}`
   * - *bz0mz*
     - :math:`\left<b_{z}\right>_{xy}`
   * - *M11z*
     - :math:`\left<{\cal M}_{11}\right>_{xy}`
   * - *M22z*
     - :math:`\left<{\cal M}_{22}\right>_{xy}`
   * - *M33z*
     - :math:`\left<{\cal M}_{33}\right>_{xy}`

Module *test_methods/testfield_nonlin_z.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *alp31*
     - :math:`\alpha_{31}`
   * - *alp12*
     - :math:`\alpha_{12}`
   * - *alp22*
     - :math:`\alpha_{22}`
   * - *alp32*
     - :math:`\alpha_{32}`
   * - *eta11*
     - :math:`\eta_{11}k`
   * - *eta21*
     - :math:`\eta_{21}k`
   * - *eta12*
     - :math:`\eta_{12}k`
   * - *eta22*
     - :math:`\eta_{22}k`
   * - *alpK*
     - :math:`\alpha^K`
   * - *alpM*
     - :math:`\alpha^M`
   * - *alpMK*
     - :math:`\alpha^{MK}`
   * - *phi11*
     - :math:`\phi_{11}`
   * - *phi21*
     - :math:`\phi_{21}`
   * - *phi12*
     - :math:`\phi_{12}`
   * - *phi22*
     - :math:`\phi_{22}`
   * - *phi32*
     - :math:`\phi_{32}`
   * - *psi11*
     - :math:`\psi_{11}k`
   * - *psi21*
     - :math:`\psi_{21}k`
   * - *psi12*
     - :math:`\psi_{12}k`
   * - *psi22*
     - :math:`\psi_{22}k`
   * - *phiK*
     - :math:`\phi^K`
   * - *phiM*
     - :math:`\phi^M`
   * - *phiMK*
     - :math:`\phi^{MK}`
   * - *alp11cc*
     - :math:`\alpha_{11}\cos^2 kz`
   * - *alp21sc*
     - :math:`\alpha_{21}\sin kz\cos kz`
   * - *alp12cs*
     - :math:`\alpha_{12}\cos kz\sin kz`
   * - *alp22ss*
     - :math:`\alpha_{22}\sin^2 kz`
   * - *eta11cc*
     - :math:`\eta_{11}\cos^2 kz`
   * - *eta21sc*
     - :math:`\eta_{21}\sin kz\cos kz`
   * - *eta12cs*
     - :math:`\eta_{12}\cos kz\sin kz`
   * - *eta22ss*
     - :math:`\eta_{22}\sin^2 kz`
   * - *s2kzDFm*
     - :math:`\left<\sin2kz\nabla\cdot F\right>`
   * - *M11*
     - :math:`{\cal M}_{11}`
   * - *M22*
     - :math:`{\cal M}_{22}`
   * - *M33*
     - :math:`{\cal M}_{33}`
   * - *M11cc*
     - :math:`{\cal M}_{11}\cos^2 kz`
   * - *M11ss*
     - :math:`{\cal M}_{11}\sin^2 kz`
   * - *M22cc*
     - :math:`{\cal M}_{22}\cos^2 kz`
   * - *M22ss*
     - :math:`{\cal M}_{22}\sin^2 kz`
   * - *M12cs*
     - :math:`{\cal M}_{12}\cos kz\sin kz`
   * - *bx11pt*
     - :math:`b_x^{11}`
   * - *bx21pt*
     - :math:`b_x^{21}`
   * - *bx12pt*
     - :math:`b_x^{12}`
   * - *bx22pt*
     - :math:`b_x^{22}`
   * - *bx0pt*
     - :math:`b_x^{0}`
   * - *by11pt*
     - :math:`b_y^{11}`
   * - *by21pt*
     - :math:`b_y^{21}`
   * - *by12pt*
     - :math:`b_y^{12}`
   * - *by22pt*
     - :math:`b_y^{22}`
   * - *by0pt*
     - :math:`b_y^{0}`
   * - *u11rms*
     - :math:`\left<u_{11}^2\right>^{1/2}`
   * - *u21rms*
     - :math:`\left<u_{21}^2\right>^{1/2}`
   * - *u12rms*
     - :math:`\left<u_{12}^2\right>^{1/2}`
   * - *u22rms*
     - :math:`\left<u_{22}^2\right>^{1/2}`
   * - *j11rms*
     - :math:`\left<j_{11}^2\right>^{1/2}`
   * - *b11rms*
     - :math:`\left<b_{11}^2\right>^{1/2}`
   * - *b21rms*
     - :math:`\left<b_{21}^2\right>^{1/2}`
   * - *b12rms*
     - :math:`\left<b_{12}^2\right>^{1/2}`
   * - *b22rms*
     - :math:`\left<b_{22}^2\right>^{1/2}`
   * - *ux0m*
     - :math:`\left<u_{0_x}\right>`
   * - *uy0m*
     - :math:`\left<u_{0_y}\right>`
   * - *ux11m*
     - :math:`\left<u_{11_x}\right>`
   * - *uy11m*
     - :math:`\left<u_{11_y}\right>`
   * - *u0rms*
     - :math:`\left<u_{0}^2\right>^{1/2}`
   * - *b0rms*
     - :math:`\left<b_{0}^2\right>^{1/2}`
   * - *u0max*
     - :math:`\operatorname{max}\left|\boldsymbol{u}_{0}\right|`
   * - *b0max*
     - :math:`\operatorname{max}\left|\boldsymbol{b}_{0}\right|`
   * - *jb0m*
     - :math:`\left<j_0\cdot b_0\right>`
   * - *ub0m*
     - :math:`\left<u_0\cdot b_0\right>`
   * - *uj0m*
     - :math:`\left<u_0\cdot j_0\right>`
   * - *E11rms*
     - :math:`\left<{\cal E}_{11}^2\right>^{1/2}`
   * - *E21rms*
     - :math:`\left<{\cal E}_{21}^2\right>^{1/2}`
   * - *E12rms*
     - :math:`\left<{\cal E}_{12}^2\right>^{1/2}`
   * - *E22rms*
     - :math:`\left<{\cal E}_{22}^2\right>^{1/2}`
   * - *E0rms*
     - :math:`\left<{\cal E}_{0}^2\right>^{1/2}`
   * - *Ex11pt*
     - :math:`{\cal E}_x^{11}`
   * - *Ex21pt*
     - :math:`{\cal E}_x^{21}`
   * - *Ex12pt*
     - :math:`{\cal E}_x^{12}`
   * - *Ex22pt*
     - :math:`{\cal E}_x^{22}`
   * - *Ex0pt*
     - :math:`{\cal E}_x^{0}`
   * - *Ey11pt*
     - :math:`{\cal E}_y^{11}`
   * - *Ey21pt*
     - :math:`{\cal E}_y^{21}`
   * - *Ey12pt*
     - :math:`{\cal E}_y^{12}`
   * - *Ey22pt*
     - :math:`{\cal E}_y^{22}`
   * - *Ey0pt*
     - :math:`{\cal E}_y^{0}`
   * - *bamp*
     - bamp
   * - *E111z*
     - :math:`{\cal E}_1^{11}`
   * - *E211z*
     - :math:`{\cal E}_2^{11}`
   * - *E311z*
     - :math:`{\cal E}_3^{11}`
   * - *E121z*
     - :math:`{\cal E}_1^{21}`
   * - *E221z*
     - :math:`{\cal E}_2^{21}`
   * - *E321z*
     - :math:`{\cal E}_3^{21}`
   * - *E112z*
     - :math:`{\cal E}_1^{12}`
   * - *E212z*
     - :math:`{\cal E}_2^{12}`
   * - *E312z*
     - :math:`{\cal E}_3^{12}`
   * - *E122z*
     - :math:`{\cal E}_1^{22}`
   * - *E222z*
     - :math:`{\cal E}_2^{22}`
   * - *E322z*
     - :math:`{\cal E}_3^{22}`
   * - *E10z*
     - :math:`{\cal E}_1^{0}`
   * - *E20z*
     - :math:`{\cal E}_2^{0}`
   * - *E30z*
     - :math:`{\cal E}_3^{0}`
   * - *EBpq*
     - :math:`{\cal E}\cdot\Bv^{pq}`
   * - *E0Um*
     - :math:`{\cal E}^0\cdot\Uv`
   * - *E0Wm*
     - :math:`{\cal E}^0\cdot\Wv`
   * - *bx0mz*
     - :math:`\left<b_{x}\right>_{xy}`
   * - *by0mz*
     - :math:`\left<b_{y}\right>_{xy}`
   * - *bz0mz*
     - :math:`\left<b_{z}\right>_{xy}`
   * - *M11z*
     - :math:`\left<{\cal M}_{11}\right>_{xy}`
   * - *M22z*
     - :math:`\left<{\cal M}_{22}\right>_{xy}`
   * - *M33z*
     - :math:`\left<{\cal M}_{33}\right>_{xy}`

Module *test_methods/testfield_x.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *alp31*
     - :math:`\alpha_{31}`
   * - *alp12*
     - :math:`\alpha_{12}`
   * - *alp22*
     - :math:`\alpha_{22}`
   * - *alp32*
     - :math:`\alpha_{32}`
   * - *eta11*
     - :math:`\eta_{11}k`
   * - *eta21*
     - :math:`\eta_{21}k`
   * - *eta12*
     - :math:`\eta_{12}k`
   * - *eta22*
     - :math:`\eta_{22}k`
   * - *alp11cc*
     - :math:`\alpha_{11}\cos^2 kx`
   * - *alp21sc*
     - :math:`\alpha_{21}\sin kx\cos kx`
   * - *alp12cs*
     - :math:`\alpha_{12}\cos kx\sin kx`
   * - *alp22ss*
     - :math:`\alpha_{22}\sin^2 kx`
   * - *eta11cc*
     - :math:`\eta_{11}\cos^2 kx`
   * - *eta21sc*
     - :math:`\eta_{21}\sin kx\cos kx`
   * - *eta12cs*
     - :math:`\eta_{12}\cos kx\sin kx`
   * - *eta22ss*
     - :math:`\eta_{22}\sin^2 kx`
   * - *alp11_x*
     - :math:`\alpha_{11}x`
   * - *alp21_x*
     - :math:`\alpha_{21}x`
   * - *alp12_x*
     - :math:`\alpha_{12}x`
   * - *alp22_x*
     - :math:`\alpha_{22}x`
   * - *eta11_x*
     - :math:`\eta_{11}kx`
   * - *eta21_x*
     - :math:`\eta_{21}kx`
   * - *eta12_x*
     - :math:`\eta_{12}kx`
   * - *eta22_x*
     - :math:`\eta_{22}kx`
   * - *alp11_x2*
     - :math:`\alpha_{11}x^2`
   * - *alp21_x2*
     - :math:`\alpha_{21}x^2`
   * - *alp12_x2*
     - :math:`\alpha_{12}x^2`
   * - *alp22_x2*
     - :math:`\alpha_{22}x^2`
   * - *eta11_x2*
     - :math:`\eta_{11}kx^2`
   * - *eta21_x2*
     - :math:`\eta_{21}kx^2`
   * - *eta12_x2*
     - :math:`\eta_{12}kx^2`
   * - *eta22_x2*
     - :math:`\eta_{22}kx^2`
   * - *b11rms*
     - :math:`\left<b_{11}^2\right>^{1/2}`
   * - *b21rms*
     - :math:`\left<b_{21}^2\right>^{1/2}`
   * - *b12rms*
     - :math:`\left<b_{12}^2\right>^{1/2}`
   * - *b22rms*
     - :math:`\left<b_{22}^2\right>^{1/2}`
   * - *b0rms*
     - :math:`\left<b_{0}^2\right>^{1/2}`
   * - *E11rms*
     - :math:`\left<{\cal E}_{11}^2\right>^{1/2}`
   * - *E21rms*
     - :math:`\left<{\cal E}_{21}^2\right>^{1/2}`
   * - *E12rms*
     - :math:`\left<{\cal E}_{12}^2\right>^{1/2}`
   * - *E22rms*
     - :math:`\left<{\cal E}_{22}^2\right>^{1/2}`
   * - *E0rms*
     - :math:`\left<{\cal E}_{0}^2\right>^{1/2}`
   * - *E111z*
     - :math:`{\cal E}_1^{11}`
   * - *E211z*
     - :math:`{\cal E}_2^{11}`
   * - *E311z*
     - :math:`{\cal E}_3^{11}`
   * - *E121z*
     - :math:`{\cal E}_1^{21}`
   * - *E221z*
     - :math:`{\cal E}_2^{21}`
   * - *E321z*
     - :math:`{\cal E}_3^{21}`
   * - *E112z*
     - :math:`{\cal E}_1^{12}`
   * - *E212z*
     - :math:`{\cal E}_2^{12}`
   * - *E312z*
     - :math:`{\cal E}_3^{12}`
   * - *E122z*
     - :math:`{\cal E}_1^{22}`
   * - *E222z*
     - :math:`{\cal E}_2^{22}`
   * - *E322z*
     - :math:`{\cal E}_3^{22}`
   * - *E10z*
     - :math:`{\cal E}_1^{0}`
   * - *E20z*
     - :math:`{\cal E}_2^{0}`
   * - *E30z*
     - :math:`{\cal E}_3^{0}`
   * - *EBpq*
     - :math:`{\cal E}\cdot\Bv^{pq}`
   * - *bx0mz*
     - :math:`\left<b_{x}\right>_{xy}`
   * - *by0mz*
     - :math:`\left<b_{y}\right>_{xy}`
   * - *bz0mz*
     - :math:`\left<b_{z}\right>_{xy}`
   * - *alp11x*
     - :math:`\alpha_{11}(x,t)`
   * - *alp21x*
     - :math:`\alpha_{21}(x,t)`
   * - *alp12x*
     - :math:`\alpha_{12}(x,t)`
   * - *alp22x*
     - :math:`\alpha_{22}(x,t)`
   * - *eta11x*
     - :math:`\eta_{11}(x,t)`
   * - *eta21x*
     - :math:`\eta_{21}(x,t)`
   * - *eta12x*
     - :math:`\eta_{12}(x,t)`
   * - *eta22x*
     - :math:`\eta_{22}(x,t)`

Module *test_methods/testfield_z.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *alp11*
     - :math:`\alpha_{11}`
   * - *alp21*
     - :math:`\alpha_{21}`
   * - *alp31*
     - :math:`\alpha_{31}`
   * - *alp12*
     - :math:`\alpha_{12}`
   * - *alp22*
     - :math:`\alpha_{22}`
   * - *alp32*
     - :math:`\alpha_{32}`
   * - *alp13*
     - :math:`\alpha_{13}`
   * - *alp23*
     - :math:`\alpha_{23}`
   * - *eta11*
     - :math:`\eta_{113}k`  or :math:`\eta_{11}k`  if leta_rank2=T
   * - *eta21*
     - :math:`\eta_{213}k`  or :math:`\eta_{21}k`  if leta_rank2=T
   * - *eta31*
     - :math:`\eta_{313}k`
   * - *eta12*
     - :math:`\eta_{123}k`  or :math:`\eta_{12}k`  if leta_rank2=T
   * - *eta22*
     - :math:`\eta_{223}k`  or :math:`\eta_{22}k`  if leta_rank2=T
   * - *eta32*
     - :math:`\eta_{323}k`
   * - *alp11cc*
     - :math:`\alpha_{11}\cos^2 kz`
   * - *alp21sc*
     - :math:`\alpha_{21}\sin kz\cos kz`
   * - *alp12cs*
     - :math:`\alpha_{12}\cos kz\sin kz`
   * - *alp22ss*
     - :math:`\alpha_{22}\sin^2 kz`
   * - *eta11cc*
     - :math:`\eta_{11}\cos^2 kz`
   * - *eta21sc*
     - :math:`\eta_{21}\sin kz\cos kz`
   * - *eta12cs*
     - :math:`\eta_{12}\cos kz\sin kz`
   * - *eta22ss*
     - :math:`\eta_{22}\sin^2 kz`
   * - *s2kzDFm*
     - :math:`\left<\sin2kz\nabla\cdot F\right>`
   * - *M11*
     - :math:`{\cal M}_{11}`
   * - *M22*
     - :math:`{\cal M}_{22}`
   * - *M33*
     - :math:`{\cal M}_{33}`
   * - *M11cc*
     - :math:`{\cal M}_{11}\cos^2 kz`
   * - *M11ss*
     - :math:`{\cal M}_{11}\sin^2 kz`
   * - *M22cc*
     - :math:`{\cal M}_{22}\cos^2 kz`
   * - *M22ss*
     - :math:`{\cal M}_{22}\sin^2 kz`
   * - *M12cs*
     - :math:`{\cal M}_{12}\cos kz\sin kz`
   * - *bx11pt*
     - :math:`b_x^{11}`
   * - *bx21pt*
     - :math:`b_x^{21}`
   * - *bx12pt*
     - :math:`b_x^{12}`
   * - *bx22pt*
     - :math:`b_x^{22}`
   * - *bx0pt*
     - :math:`b_x^{0}`
   * - *by11pt*
     - :math:`b_y^{11}`
   * - *by21pt*
     - :math:`b_y^{21}`
   * - *by12pt*
     - :math:`b_y^{12}`
   * - *by22pt*
     - :math:`b_y^{22}`
   * - *by0pt*
     - :math:`b_y^{0}`
   * - *b11rms*
     - :math:`\left<b_{11}^2\right>^{1/2}`
   * - *b21rms*
     - :math:`\left<b_{21}^2\right>^{1/2}`
   * - *b12rms*
     - :math:`\left<b_{12}^2\right>^{1/2}`
   * - *b22rms*
     - :math:`\left<b_{22}^2\right>^{1/2}`
   * - *b0rms*
     - :math:`\left<b_{0}^2\right>^{1/2}`
   * - *jb0m*
     - :math:`\left<jb_{0}\right>`
   * - *E11rms*
     - :math:`\left<{\cal E}_{11}^2\right>^{1/2}`
   * - *E21rms*
     - :math:`\left<{\cal E}_{21}^2\right>^{1/2}`
   * - *E12rms*
     - :math:`\left<{\cal E}_{12}^2\right>^{1/2}`
   * - *E22rms*
     - :math:`\left<{\cal E}_{22}^2\right>^{1/2}`
   * - *E0rms*
     - :math:`\left<{\cal E}_{0}^2\right>^{1/2}`
   * - *Ex11pt*
     - :math:`{\cal E}_x^{11}`
   * - *Ex21pt*
     - :math:`{\cal E}_x^{21}`
   * - *Ex12pt*
     - :math:`{\cal E}_x^{12}`
   * - *Ex22pt*
     - :math:`{\cal E}_x^{22}`
   * - *Ex0pt*
     - :math:`{\cal E}_x^{0}`
   * - *Ey11pt*
     - :math:`{\cal E}_y^{11}`
   * - *Ey21pt*
     - :math:`{\cal E}_y^{21}`
   * - *Ey12pt*
     - :math:`{\cal E}_y^{12}`
   * - *Ey22pt*
     - :math:`{\cal E}_y^{22}`
   * - *Ey0pt*
     - :math:`{\cal E}_y^{0}`
   * - *bamp*
     - bamp
   * - *alp11z*
     - :math:`\alpha_{11}(z,t)`
   * - *alp21z*
     - :math:`\alpha_{21}(z,t)`
   * - *alp12z*
     - :math:`\alpha_{12}(z,t)`
   * - *alp22z*
     - :math:`\alpha_{22}(z,t)`
   * - *alp13z*
     - :math:`\alpha_{13}(z,t)`
   * - *alp23z*
     - :math:`\alpha_{23}(z,t)`
   * - *eta11z*
     - :math:`\eta_{11}(z,t)`
   * - *eta21z*
     - :math:`\eta_{21}(z,t)`
   * - *eta12z*
     - :math:`\eta_{12}(z,t)`
   * - *eta22z*
     - :math:`\eta_{22}(z,t)`
   * - *uzjx1z*
     - :math:`u_z j^{11}_x`
   * - *uzjy1z*
     - :math:`u_z j^{11}_y`
   * - *uzjz1z*
     - :math:`u_z j^{11}_z`
   * - *uzjx2z*
     - :math:`u_z j^{21}_x`
   * - *uzjy2z*
     - :math:`u_z j^{21}_y`
   * - *uzjz2z*
     - :math:`u_z j^{21}_z`
   * - *uzjx3z*
     - :math:`u_z j^{12}_x`
   * - *uzjy3z*
     - :math:`u_z j^{12}_y`
   * - *uzjz3z*
     - :math:`u_z j^{12}_z`
   * - *uzjx4z*
     - :math:`u_z j^{22}_x`
   * - *uzjy4z*
     - :math:`u_z j^{22}_y`
   * - *uzjz4z*
     - :math:`u_z j^{22}_z`
   * - *E111z*
     - :math:`{\cal E}_1^{11}`
   * - *E211z*
     - :math:`{\cal E}_2^{11}`
   * - *E311z*
     - :math:`{\cal E}_3^{11}`
   * - *E121z*
     - :math:`{\cal E}_1^{21}`
   * - *E221z*
     - :math:`{\cal E}_2^{21}`
   * - *E321z*
     - :math:`{\cal E}_3^{21}`
   * - *E112z*
     - :math:`{\cal E}_1^{12}`
   * - *E212z*
     - :math:`{\cal E}_2^{12}`
   * - *E312z*
     - :math:`{\cal E}_3^{12}`
   * - *E122z*
     - :math:`{\cal E}_1^{22}`
   * - *E222z*
     - :math:`{\cal E}_2^{22}`
   * - *E322z*
     - :math:`{\cal E}_3^{22}`
   * - *E10z*
     - :math:`{\cal E}_1^{0}`
   * - *E20z*
     - :math:`{\cal E}_2^{0}`
   * - *E30z*
     - :math:`{\cal E}_3^{0}`
   * - *EBpq*
     - :math:`{\cal E}\cdot\Bv^{pq}`
   * - *E0Um*
     - :math:`{\cal E}^0\cdot\Uv`
   * - *E0Wm*
     - :math:`{\cal E}^0\cdot\Wv`
   * - *bx0mz*
     - :math:`\left<b_{x}\right>_{xy}`
   * - *by0mz*
     - :math:`\left<b_{y}\right>_{xy}`
   * - *bz0mz*
     - :math:`\left<b_{z}\right>_{xy}`
   * - *M11z*
     - :math:`\left<{\cal M}_{11}\right>_{xy}`
   * - *M22z*
     - :math:`\left<{\cal M}_{22}\right>_{xy}`
   * - *M33z*
     - :math:`\left<{\cal M}_{33}\right>_{xy}`

Module *test_methods/testscalar.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *gam11*
     - :math:`\gamma_{1}^{(1)}`
   * - *gam12*
     - :math:`\gamma_{2}^{(1)}`
   * - *gam13*
     - :math:`\gamma_{3}^{(1)}`
   * - *gam21*
     - :math:`\gamma_{1}^{(2)}`
   * - *gam22*
     - :math:`\gamma_{2}^{(2)}`
   * - *gam23*
     - :math:`\gamma_{3}^{(2)}`
   * - *gam31*
     - :math:`\gamma_{1}^{(3)}`
   * - *gam32*
     - :math:`\gamma_{2}^{(3)}`
   * - *gam33*
     - :math:`\gamma_{3}^{(3)}`
   * - *kap11*
     - :math:`\kappa_{11}`
   * - *kap21*
     - :math:`\kappa_{21}`
   * - *kap31*
     - :math:`\kappa_{31}`
   * - *kap12*
     - :math:`\kappa_{12}`
   * - *kap22*
     - :math:`\kappa_{22}`
   * - *kap32*
     - :math:`\kappa_{32}`
   * - *kap13*
     - :math:`\kappa_{13}`
   * - *kap23*
     - :math:`\kappa_{23}`
   * - *kap33*
     - :math:`\kappa_{33}`
   * - *gam11z*
     - :math:`\gamma_{1}^{(1)}(z,t)`
   * - *gam12z*
     - :math:`\gamma_{2}^{(1)}(z,t)`
   * - *gam13z*
     - :math:`\gamma_{3}^{(1)}(z,t)`
   * - *gam21z*
     - :math:`\gamma_{1}^{(2)}(z,t)`
   * - *gam22z*
     - :math:`\gamma_{2}^{(2)}(z,t)`
   * - *gam23z*
     - :math:`\gamma_{3}^{(2)}(z,t)`
   * - *gam31z*
     - :math:`\gamma_{1}^{(3)}(z,t)`
   * - *gam32z*
     - :math:`\gamma_{2}^{(3)}(z,t)`
   * - *gam33z*
     - :math:`\gamma_{3}^{(3)}(z,t)`
   * - *kap11z*
     - :math:`\kappa_{11}(z,t)`
   * - *kap21z*
     - :math:`\kappa_{21}(z,t)`
   * - *kap31z*
     - :math:`\kappa_{31}(z,t)`
   * - *kap12z*
     - :math:`\kappa_{12}(z,t)`
   * - *kap22z*
     - :math:`\kappa_{22}(z,t)`
   * - *kap32z*
     - :math:`\kappa_{32}(z,t)`
   * - *kap13z*
     - :math:`\kappa_{13}(z,t)`
   * - *kap23z*
     - :math:`\kappa_{23}(z,t)`
   * - *kap33z*
     - :math:`\kappa_{33}(z,t)`
   * - *mgam33*
     - :math:`\tilde\gamma_{33}`
   * - *mkap33*
     - :math:`\tilde\kappa_{33}`
   * - *ngam33*
     - :math:`\hat\gamma_{33}`
   * - *nkap33*
     - :math:`\hat\kappa_{33}`
   * - *c1rms*
     - :math:`\left<c_{1}^2\right>^{1/2}`
   * - *c2rms*
     - :math:`\left<c_{2}^2\right>^{1/2}`
   * - *c3rms*
     - :math:`\left<c_{3}^2\right>^{1/2}`
   * - *c4rms*
     - :math:`\left<c_{4}^2\right>^{1/2}`
   * - *c5rms*
     - :math:`\left<c_{5}^2\right>^{1/2}`
   * - *c6rms*
     - :math:`\left<c_{6}^2\right>^{1/2}`
   * - *c1pt*
     - :math:`c^{1}`
   * - *c2pt*
     - :math:`c^{2}`
   * - *c3pt*
     - :math:`c^{3}`
   * - *c4pt*
     - :math:`c^{4}`
   * - *c5pt*
     - :math:`c^{5}`
   * - *c6pt*
     - :math:`c^{6}`
   * - *F11z*
     - :math:`{\cal F}_1^{1}`
   * - *F21z*
     - :math:`{\cal F}_2^{1}`
   * - *F31z*
     - :math:`{\cal F}_3^{1}`
   * - *F12z*
     - :math:`{\cal F}_1^{2}`
   * - *F22z*
     - :math:`{\cal F}_2^{2}`
   * - *F32z*
     - :math:`{\cal F}_3^{2}`

Module *test_methods/testscalar_axisym.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *muc1*
     - :math:`\mu^{(c1)}`
   * - *muc2*
     - :math:`\mu^{(c2)}`
   * - *gamc*
     - :math:`\gamma^{(c)}`
   * - *kapcPERP1*
     - :math:`\kappa_\perp^{(1)}`
   * - *kapcPERP2*
     - :math:`\kappa_\perp^{(2)}`
   * - *kapcPARA*
     - :math:`\kappa_\parallel`
   * - *mucz*
     - :math:`\mu^{(c)}(z,t)`
   * - *gamcz*
     - :math:`\gamma^{(c)}(z,t)`
   * - *kapcPERPz*
     - :math:`\kappa_\perp(z,t)`
   * - *kapcPARAz*
     - :math:`\kappa_\parallel(z,t)`
   * - *gam11*
     - :math:`\gamma_{1}^{(1)}`
   * - *gam12*
     - :math:`\gamma_{2}^{(1)}`
   * - *gam13*
     - :math:`\gamma_{3}^{(1)}`
   * - *gam21*
     - :math:`\gamma_{1}^{(2)}`
   * - *gam22*
     - :math:`\gamma_{2}^{(2)}`
   * - *gam23*
     - :math:`\gamma_{3}^{(2)}`
   * - *gam31*
     - :math:`\gamma_{1}^{(3)}`
   * - *gam32*
     - :math:`\gamma_{2}^{(3)}`
   * - *gam33*
     - :math:`\gamma_{3}^{(3)}`
   * - *kap11*
     - :math:`\kappa_{11}`
   * - *kap21*
     - :math:`\kappa_{21}`
   * - *kap31*
     - :math:`\kappa_{31}`
   * - *kap12*
     - :math:`\kappa_{12}`
   * - *kap22*
     - :math:`\kappa_{22}`
   * - *kap32*
     - :math:`\kappa_{32}`
   * - *kap13*
     - :math:`\kappa_{13}`
   * - *kap23*
     - :math:`\kappa_{23}`
   * - *kap33*
     - :math:`\kappa_{33}`
   * - *gam11z*
     - :math:`\gamma_{1}^{(1)}(z,t)`
   * - *gam12z*
     - :math:`\gamma_{2}^{(1)}(z,t)`
   * - *gam13z*
     - :math:`\gamma_{3}^{(1)}(z,t)`
   * - *gam21z*
     - :math:`\gamma_{1}^{(2)}(z,t)`
   * - *gam22z*
     - :math:`\gamma_{2}^{(2)}(z,t)`
   * - *gam23z*
     - :math:`\gamma_{3}^{(2)}(z,t)`
   * - *gam31z*
     - :math:`\gamma_{1}^{(3)}(z,t)`
   * - *gam32z*
     - :math:`\gamma_{2}^{(3)}(z,t)`
   * - *gam33z*
     - :math:`\gamma_{3}^{(3)}(z,t)`
   * - *gam3z*
     - :math:`\gamma^{(c)}(z,t)`
   * - *kap11z*
     - :math:`\kappa_{11}(z,t)`
   * - *kap21z*
     - :math:`\kappa_{21}(z,t)`
   * - *kap31z*
     - :math:`\kappa_{31}(z,t)`
   * - *kap12z*
     - :math:`\kappa_{12}(z,t)`
   * - *kap22z*
     - :math:`\kappa_{22}(z,t)`
   * - *kap32z*
     - :math:`\kappa_{32}(z,t)`
   * - *kap13z*
     - :math:`\kappa_{13}(z,t)`
   * - *kap23z*
     - :math:`\kappa_{23}(z,t)`
   * - *kap33z*
     - :math:`\kappa_{33}(z,t)`
   * - *mgam33*
     - :math:`\tilde\gamma_{33}`
   * - *mkap33*
     - :math:`\tilde\kappa_{33}`
   * - *ngam33*
     - :math:`\hat\gamma_{33}`
   * - *nkap33*
     - :math:`\hat\kappa_{33}`
   * - *c1rms*
     - :math:`\left<c_{1}^2\right>^{1/2}`
   * - *c2rms*
     - :math:`\left<c_{2}^2\right>^{1/2}`
   * - *c3rms*
     - :math:`\left<c_{3}^2\right>^{1/2}`
   * - *c4rms*
     - :math:`\left<c_{4}^2\right>^{1/2}`
   * - *c5rms*
     - :math:`\left<c_{5}^2\right>^{1/2}`
   * - *c6rms*
     - :math:`\left<c_{6}^2\right>^{1/2}`
   * - *c1pt*
     - :math:`c^{1}`
   * - *c2pt*
     - :math:`c^{2}`
   * - *c3pt*
     - :math:`c^{3}`
   * - *c4pt*
     - :math:`c^{4}`
   * - *c5pt*
     - :math:`c^{5}`
   * - *c6pt*
     - :math:`c^{6}`
   * - *F11z*
     - :math:`{\cal F}_1^{1}`
   * - *F21z*
     - :math:`{\cal F}_2^{1}`
   * - *F31z*
     - :math:`{\cal F}_3^{1}`
   * - *F12z*
     - :math:`{\cal F}_1^{2}`
   * - *F22z*
     - :math:`{\cal F}_2^{2}`
   * - *F32z*
     - :math:`{\cal F}_3^{2}`

Module *test_methods/testscalar_simple.f90*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.. list-table::
   :header-rows: 1
   :widths: 25 75

   * - Variable
     - Meaning
   * - *gam11*
     - :math:`\gamma_{1}^{(1)}`
   * - *gam12*
     - :math:`\gamma_{2}^{(1)}`
   * - *gam13*
     - :math:`\gamma_{3}^{(1)}`
   * - *gam21*
     - :math:`\gamma_{1}^{(2)}`
   * - *gam22*
     - :math:`\gamma_{2}^{(2)}`
   * - *gam23*
     - :math:`\gamma_{3}^{(2)}`
   * - *gam31*
     - :math:`\gamma_{1}^{(3)}`
   * - *gam32*
     - :math:`\gamma_{2}^{(3)}`
   * - *gam33*
     - :math:`\gamma_{3}^{(3)}`
   * - *kap11*
     - :math:`\kappa_{11}`
   * - *kap21*
     - :math:`\kappa_{21}`
   * - *kap31*
     - :math:`\kappa_{31}`
   * - *kap12*
     - :math:`\kappa_{12}`
   * - *kap22*
     - :math:`\kappa_{22}`
   * - *kap32*
     - :math:`\kappa_{32}`
   * - *kap13*
     - :math:`\kappa_{13}`
   * - *kap23*
     - :math:`\kappa_{23}`
   * - *kap33*
     - :math:`\kappa_{33}`
   * - *gam11z*
     - :math:`\gamma_{1}^{(1)}(z,t)`
   * - *gam12z*
     - :math:`\gamma_{2}^{(1)}(z,t)`
   * - *gam13z*
     - :math:`\gamma_{3}^{(1)}(z,t)`
   * - *gam21z*
     - :math:`\gamma_{1}^{(2)}(z,t)`
   * - *gam22z*
     - :math:`\gamma_{2}^{(2)}(z,t)`
   * - *gam23z*
     - :math:`\gamma_{3}^{(2)}(z,t)`
   * - *gam31z*
     - :math:`\gamma_{1}^{(3)}(z,t)`
   * - *gam32z*
     - :math:`\gamma_{2}^{(3)}(z,t)`
   * - *gam33z*
     - :math:`\gamma_{3}^{(3)}(z,t)`
   * - *kap11z*
     - :math:`\kappa_{11}(z,t)`
   * - *kap21z*
     - :math:`\kappa_{21}(z,t)`
   * - *kap31z*
     - :math:`\kappa_{31}(z,t)`
   * - *kap12z*
     - :math:`\kappa_{12}(z,t)`
   * - *kap22z*
     - :math:`\kappa_{22}(z,t)`
   * - *kap32z*
     - :math:`\kappa_{32}(z,t)`
   * - *kap13z*
     - :math:`\kappa_{13}(z,t)`
   * - *kap23z*
     - :math:`\kappa_{23}(z,t)`
   * - *kap33z*
     - :math:`\kappa_{33}(z,t)`
   * - *mgam33*
     - :math:`\tilde\gamma_{33}`
   * - *mkap33*
     - :math:`\tilde\kappa_{33}`
   * - *ngam33*
     - :math:`\hat\gamma_{33}`
   * - *nkap33*
     - :math:`\hat\kappa_{33}`
   * - *c1rms*
     - :math:`\left<c_{1}^2\right>^{1/2}`
   * - *c2rms*
     - :math:`\left<c_{2}^2\right>^{1/2}`
   * - *c3rms*
     - :math:`\left<c_{3}^2\right>^{1/2}`
   * - *c4rms*
     - :math:`\left<c_{4}^2\right>^{1/2}`
   * - *c5rms*
     - :math:`\left<c_{5}^2\right>^{1/2}`
   * - *c6rms*
     - :math:`\left<c_{6}^2\right>^{1/2}`
   * - *c1pt*
     - :math:`c^{1}`
   * - *c2pt*
     - :math:`c^{2}`
   * - *c3pt*
     - :math:`c^{3}`
   * - *c4pt*
     - :math:`c^{4}`
   * - *c5pt*
     - :math:`c^{5}`
   * - *c6pt*
     - :math:`c^{6}`
   * - *F11z*
     - :math:`{\cal F}_1^{1}`
   * - *F21z*
     - :math:`{\cal F}_2^{1}`
   * - *F31z*
     - :math:`{\cal F}_3^{1}`
   * - *F12z*
     - :math:`{\cal F}_1^{2}`
   * - *F22z*
     - :math:`{\cal F}_2^{2}`
   * - *F32z*
     - :math:`{\cal F}_3^{2}`