This page was generated from
docs/notebooks/plasma/1_equilibrium.ipynb.
Interactive online version:
.
Download notebook.
Interactive online version:
Plasma Equilibrium¶
This notebook demonstrates how to load and visualize a plasma equilibrium using the cherab.imas interface.
The example test data was calculated by JINTRAC for an ITER 15 MA H-mode scenario.
[1]:
from cherab.imas.datasets import iter_jintrac
from cherab.imas.plasma import load_equilibrium
from cherab.tools.equilibrium import plot_equilibrium
Retrieve ITER JINTRAC sample data¶
[2]:
path = iter_jintrac()
Downloading file 'iter_scenario_53298_seq1_DD4.nc' from 'doi:10.5281/zenodo.17062699/iter_scenario_53298_seq1_DD4.nc' to '/home/runner/.cache/cherab/imas'.
17:01:39 INFO Parsing data dictionary version 4.1.1 @dd_zip.py:89
=== Apply patch to fix ===
17:01:39 INFO Parsing data dictionary version 4.0.0 @dd_zip.py:89
/home/runner/work/imas/imas/src/cherab/imas/datasets/_patch.py:46: RuntimeWarning: The 'get_slice' method is not implemented for the URI '/home/runner/.cache/cherab/imas/iter_scenario_53298_seq1_DD4.nc'. Falling back to 'get' method because the returned IDS contains a single time slice.
ids_core = get_ids_time_slice(entry, "core_profiles", 0.0)
Load the equilibrium data¶
The equilibrium information is loaded from the IMAS equilibrium IDS.
[3]:
equilibrium = load_equilibrium(path, "r")
/home/runner/work/imas/imas/src/cherab/imas/plasma/equilibrium.py:97: RuntimeWarning: The 'get_slice' method is not implemented for the URI '/home/runner/.cache/cherab/imas/iter_scenario_53298_seq1_DD4_mod.nc'. Falling back to 'get' method because the returned IDS contains a single time slice.
equilibrium_ids = get_ids_time_slice(
Plot several quantities of the equilibrium¶
Using cherab’s built-in plotting utilities, we can visualize the magnetic equilibrium.
[4]:
plot_equilibrium(equilibrium, detail=True)
Sampling psi...
Plotting summary...
Sampling B-field...
Sampling LCFS interior...
Calculating B-field magnitude...
Sampling q...
Plotting details...