Transport properties of the reference state

These files define the radial transport properties of the reference state. These arrays are calculated in the subroutines radial and transportProperties. The output files are written in the subroutine preCalc.

anel.TAG

Note

This output is only calculated when an anelastic model is run, that is when l_anel=.true. or l_anelastic_liquid=.true..

This file contains the radial profiles of the reference state (density, temperature, gravity, etc.).

No. of column

Contents

1

radial level: \(r\)

2

temperature: \(\tilde{T}(r)\)

3

density: \(\tilde{\rho}(r)\)

4

radial derivative of the log of the density: \(\beta={\rm d} \ln\tilde{\rho}/{\rm d} r\)

5

radial derivative of \(\beta\): \({\rm d} \beta/{\rm d} r\)

6

gravity: \(g(r)\)

7

entropy gradient: \({\rm d} s_0/{\rm d} r\)

8

thermal diffusion operator: \(\nabla \cdot (K(r)\tilde{T}(r)\nabla s_0)\)

9

inverse of the Gruneisen parameter :math`1/Gamma`: \((\partial\ln\tilde{\rho}/\partial\ln\tilde{T})_S\)

10

radial derivative of the log of temperature: \(\beta={\rm d} \ln\tilde{T}/{\rm d} r\)

This file can be read using MagicRadial with the following options:

>>> rad = MagicRadial(field='anel')
>>> # print radius and density
>>> print(rad.radius, rad.rho0)

varCond.TAG

Note

This output is only calculated when the electrical conductivity varies with radius, i.e. when nVarCond /= 0

This file contains the radial profiles of the electrical conductivity, the electrical diffusivity and its radial derivative.

No. of column

Contents

1

radial level: \(r\)

2

electrical conductivity: \(\sigma(r)\)

3

electrical diffusivity: \(\lambda(r)=1/\sigma(r)\)

4

radial derivative of the electrical diffusivity: \({\rm d} \ln\lambda/{\rm d} r\)

This file can be read using MagicRadial with the following options:

>>> rad = MagicRadial(field='varCond')
>>> print(rad.conduc) # Electrical conductivity

varDiff.TAG

Note

This output is only calculated when the thermal diffusivity varies with radius, i.e. when nVarDiff /= 0

This file contains the radial profiles of the thermal conductivity, the thermal diffusivity and its radial derivative.

No. of column

Contents

1

radial level: \(r\)

2

thermal conductivity: \(K(r)\)

3

thermal diffusivity: \(\kappa(r)=K(r)/\tilde{\rho}(r)\)

4

radial derivative of the electrical diffusivity: \({\rm d} \ln\kappa/{\rm d} r\)

5

Prandtl number: \(Pr(r)=\nu(r)/\kappa(r)\)

This file can be read using MagicRadial with the following options:

>>> rad = MagicRadial(field='varDiff')
>>> print(rad.kappa) # Thermal diffusivity

varVisc.TAG

Note

This output is only calculated when the kinematic viscosity varies with radius, i.e. when nVarVisc /= 0

This file contains the radial profiles of the dynamic viscosity, the kinematic viscosity and its radial derivative.

No. of column

Contents

1

radial level: \(r\)

2

dynamic viscosity: \(\mu(r)\)

3

kinetmatic viscosity: \(\nu(r)=\mu(r)/\tilde{\rho}(r)\)

4

radial derivative of the kinematic viscosity: \({\rm d} \ln\nu/{\rm d} r\)

5

Prandtl number: \(Pr(r)=\nu(r)/\kappa(r)\)

6

magnetic Prandtl number \(Pm(r)=\nu(r)/\lambda(r)\)

This file can be read using MagicRadial with the following options:

>>> rad = MagicRadial(field='varVisc')
>>> # print kinematic viscosity and Ekman
>>> print(rad.kinVisc, rad.ekman)

Nonlinear mapping of the Chebyshev grid

rNM.TAG

Note

This file is only written when l_newmap=.true..

This file contains the profile of the radial mapping and its derivatives:

No. of column

Contents

1

Grid point index

2

Radius of a grid point

3

First derivative of the mapping at a grid point

4

Second derivative of the mapping at a grid point

5

Third derivative of the mapping at a grid point