Support for time-averaged radial profiles

class magic.MagicRadial(datadir='.', field='eKin', iplot=True, tag=None, tags=None, normalize_radius=False, quiet=False)

This class can be used to read and display the time and horizontally averaged files:

• Kinetic energy: eKinR.TAG

• Magnetic energy: eMagR.TAG

• Anelastic reference state: anel.TAG

• Variable electrical conductivity: varCond.TAG

• Variable thermal diffusivity: varDiff.TAG

• Variable kinematic viscosity: varVisc.TAG

• Diagnostic parameters: parR.TAG

• Power budget: powerR.TAG

• Phase field: phiR.TAG

• Heat fluxes: fluxesR.TAG

• Mean entropy, temperature and pressure: heatR.TAG

• Radial profiles used for boundary layers: bLayersR.TAG

• Parallel/perpendicular decomposition: perpParR.TAG

>>> rad = MagicRadial(field='eKinR') # display the content of eKinR.tag
>>> print(rad.radius, rad.ekin_pol_axi) # print radius and poloidal energy

__init__(datadir='.', field='eKin', iplot=True, tag=None, tags=None, normalize_radius=False, quiet=False)

Parameters:

• datadir (str) – working directory

• field (str) – the field you want to plot

• iplot (bool) – to plot the output, default is True

• tag (str) – a specific tag, default is None

• tags (list) – a list that contains multiple tags: useful to sum several radial files

• quiet (bool) – when set to True, makes the output silent (default False)

plot()

Display the result when iplot=True