Heatflux post processing

Documenation of the EMC3-EIRENE is part of the EMC3-EIRENE documentation

• map particle and density deposition ENERGY_DEPO and PARTICLE_DEPO to the target structures

• Target heatflux contains contribution from potential energy

• Neutrals are not included

• plasma radiation not included

# Ensure the example data is present
from get_data import load_example_data

_ = load_example_data()

# Input file
! cat ../../example-data/emc3_example/fort.3

# Input file
! cat ../../example-data/emc3_example/fort.79

Run with e.g.

cd /raven/u/mp004/runs/casegroup/case_docs/
srun emc3_eirene

# xemc3 can be used for plotting:
! xemc3-divertor ../../example-data/emc3_example/ -gls -t 'Example data'

Analysing the data

import xemc3

plates = xemc3.load.plates("../../example-data/emc3_example/")

plates

The data contains padding such that they can be combined into a numpy-data-block

plates["f_E"].isel(plate_ind=13).plot()

Using the ds.emc3[] the padding can be removed, if it not not needed for the current data, i.e. if only a single plate is included:

nosym = list(plates.emc3.iter_plates())
nosym[13].emc3["f_E"].plot()
nosym[13]

It is possible to get all plates, rather then just the simulated half-module using the ds.emc3.iter_plates function:

sym = list(plates.emc3.iter_plates(symmetry=True, segments=5))
len(sym)