⚡ Enhance core functionality and documentation improvements

CHANGELOG.md

@@ -5,6 +5,22 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.5.1] - 2026-07-28
+
+### Added
+
+- Add synthetic JINTRAC radiation values dataset support and related regression tests
+- Extend radiation emitter loading to support dual emissivity sources with improved validation
+
+### Changed
+
+- Improve radiation emitter loading checks for duplicate emissivity values and core-profile grid data
+
+### Fixed
+
+- Fix grid data loading checks for radiation core profiles
+- Improve error handling in radiation emitter loading workflows
+
 ## [0.5.0] - 2026-06-24
 
 ### Added

README.md

@@ -49,6 +49,9 @@ pip install cherab-imas
 mamba install -c conda-forge cherab-imas
 ```
 
+> [!NOTE]
+> Some IMAS features that rely on the memory backend require a working `imas_core` runtime. On platforms where that backend is unavailable (for example, some macOS Intel environments), a subset of tests that exercise those features may be skipped automatically, while the core package functionality remains available.
+
 ## 📝 Documentation
 
 See the [official documentation](https://cherab.github.io/imas/) to learn more.

docs/notebooks/plasma/4_emission.ipynb

@@ -721,7 +721,7 @@
     "    title=\"Ray-traced emission spectrum integrated over wavelength bands\",\n",
     "    ylabel=\"Radiance [W/m²/sr]\",\n",
     "    xlocator=range(len(wavelength_bands)),\n",
-    "    xticklabels=[f\"{band[0]} – {band[1]}\" for band in wavelength_bands],\n",
+    "    xticklabels=[f\"{band[0]:g} – {band[1]:g}\" for band in wavelength_bands],\n",
     ")\n",
     "\n",
     "axs.format(\n",

docs/notebooks/radiation/radiation_3d.ipynb

@@ -21,12 +21,8 @@
    "source": [
     "import numpy as np\n",
     "import ultraplot as uplt\n",
-    "from imas import DBEntry\n",
     "from raysect.optical import World\n",
-    "from rich import print as rprint\n",
-    "from rich.table import Table\n",
     "\n",
-    "from cherab.core.math import sample3d_grid\n",
     "from cherab.imas.datasets import iter_jorek\n",
     "from cherab.imas.emitter import load_radiation_emitter\n",
     "\n",

docs/source/api.md

@@ -1,4 +1,6 @@
-# API Reference
+(api-reference)=
+
+# 📚 API Reference
 
 This page contains auto-generated API reference documentation.
 

docs/source/contributing.md

@@ -1,4 +1,6 @@
-# Contributing
+(contributing)=
+
+# 🤝 Contributing
 
 We welcome contributions to `cherab-imas`! Whether you're fixing bugs, adding new features, or improving documentation, your help is appreciated.
 
@@ -17,7 +19,7 @@ If you encounter any issues or have suggestions for improvements, please open an
 
 ### Prerequisites
 
-Development tasks are managed using `pixi`. If you don't have `pixi` installed, please refer to the https://pixi.sh documentation for installation instructions.
+Development tasks are managed using `pixi`. If you don't have `pixi` installed, please refer to the <https://pixi.sh> documentation for installation instructions.
 Other tools like `git` can be installed globally via `pixi`:
 
 ```bash
@@ -74,7 +76,7 @@ If you want to host the documentation locally, you can do so with:
 pixi run doc-serve
 ```
 
-The documentation will be served at http://localhost:8000.
+The documentation will be served at <http://localhost:8000>.
 
 :::
 :::{md-tab-item} lint/format

docs/source/index.md

@@ -66,7 +66,7 @@ examples.md
 ```
 
 ```{toctree}
-:caption: Reference 📖
+:caption: Reference
 :hidden:
 :maxdepth: 1
 

src/cherab/imas/datasets/__init__.py

@@ -46,14 +46,15 @@
 the internet connectivity.
 """
 
-from ._builtin import bolometer_moc
+from ._builtin import bolometer_moc, iter_jintrac_radiation_values
 from ._fetchers import iter_jintrac, iter_jorek, iter_solps
 from ._utils import clear_cache
 
 __all__ = [
     "iter_jintrac",
     "iter_solps",
     "iter_jorek",
+    "iter_jintrac_radiation_values",
     "bolometer_moc",
     "clear_cache",
 ]

src/cherab/imas/datasets/_builtin.py

@@ -1,13 +1,17 @@
 """Provide functionality to create builtin IMAS sample datasets."""
 
 import datetime
+from pathlib import Path
 
 import numpy as np
 import pooch
 from raysect.core.math import Point3D, Vector3D, rotate_z, to_cylindrical
 
 from imas import DBEntry, IDSFactory
 from imas.ids_defs import IDS_TIME_MODE_HOMOGENEOUS
+from imas.ids_toplevel import IDSToplevel
+
+from ..ids.common.ggd import load_grid
 
 N_CH = 5  # Number of channels per camera
 N_APERTURE = 3  # Number of apertures per channel (for collimator cameras)
@@ -27,6 +31,79 @@
 Y_AXIS = Vector3D(0, 1, 0)
 
 
+def _iter_jintrac_radiation_values_data() -> tuple[IDSToplevel, IDSToplevel]:
+    """Create synthetic radiation/equilibrium IDS objects for values-based emitter tests.
+
+    Returns
+    -------
+    tuple[IDSToplevel, IDSToplevel]
+        A pair containing ``(equilibrium_ids, radiation_ids)``.
+
+    Raises
+    ------
+    RuntimeError
+        If the source JINTRAC dataset does not contain required equilibrium/core/edge data.
+    """
+    from ._fetchers import iter_jintrac
+
+    with DBEntry(iter_jintrac(), "r") as entry:
+        equilibrium = entry.get("equilibrium", autoconvert=False)
+        core_profiles = entry.get("core_profiles", autoconvert=False)
+        edge_profiles = entry.get("edge_profiles", autoconvert=False)
+
+    if not len(equilibrium.time):
+        raise RuntimeError("The source equilibrium IDS has no time slices.")
+    if not len(core_profiles.profiles_1d):
+        raise RuntimeError("The source core_profiles IDS has no 1D profile data.")
+    if not len(edge_profiles.grid_ggd):
+        raise RuntimeError("The source edge_profiles IDS has no GGD grid.")
+
+    src_grid = core_profiles.profiles_1d[0].grid
+    rho_tor_norm = np.asarray(src_grid.rho_tor_norm, dtype=np.float64)
+    psi = np.asarray(src_grid.psi, dtype=np.float64)
+
+    if rho_tor_norm.size == 0 or psi.size == 0:
+        raise RuntimeError("The source core_profiles grid does not contain psi/rho_tor_norm data.")
+
+    radiation = IDSFactory(equilibrium._version).new("radiation")
+    radiation.ids_properties.homogeneous_time = equilibrium.ids_properties.homogeneous_time
+    radiation.ids_properties.comment = "Synthetic radiation IDS for CHERAB-IMAS values tests"
+    radiation.ids_properties.creation_date = datetime.date.today().isoformat()
+    radiation.time = np.asarray(equilibrium.time, dtype=np.float64)
+
+    radiation.grid_ggd.resize(1)
+    radiation.grid_ggd[0] = edge_profiles.grid_ggd[0]
+
+    radiation.process.resize(1)
+    process = radiation.process[0]
+    process.identifier.index = 901  # custom emission, referenced: https://imas-data-dictionary.readthedocs.io/en/latest/generated/identifier/radiation_identifier.html
+    process.identifier.name = "total"
+    process.profiles_1d.resize(1)
+    process.profiles_1d[0].grid.rho_tor_norm = rho_tor_norm
+    process.profiles_1d[0].grid.psi = psi
+
+    # Smooth, strictly positive core emissivity profile.
+    core_values = 2.0e5 * (1.0 - 0.7 * np.clip(rho_tor_norm, 0.0, 1.0) ** 1.5) + 2.0e4
+    process.profiles_1d[0].electrons.emissivity = core_values
+
+    _, subsets, subset_id = load_grid(radiation.grid_ggd[0], with_subsets=True)
+    subset_name = next((name for name, index in subset_id.items() if index == 5), None)
+    if subset_name is None:
+        raise RuntimeError("Unable to find GGD subset id=5 (cells) in source grid_ggd.")
+
+    num_cells = len(subsets[subset_name])
+    if num_cells == 0:
+        raise RuntimeError("The selected GGD subset (id=5) contains no cells.")
+
+    edge_values = np.linspace(4.0e4, 1.0e5, num_cells, dtype=np.float64)
+    process.ggd.resize(1)
+    process.ggd[0].electrons.emissivity.resize(1)
+    process.ggd[0].electrons.emissivity[0].grid_subset_index = 5
+    process.ggd[0].electrons.emissivity[0].values = edge_values
+
+    return equilibrium, radiation
+
+
 def _bolo_data():
     """
     Create a mock bolometer IDS dataset.
@@ -283,3 +360,40 @@ def bolometer_moc() -> str:
             entry.put(ids)
 
     return str(path)
+
+
+def iter_jintrac_radiation_values() -> str:
+    """Return a synthetic radiation dataset for values-based emitter examples/tests.
+
+    This builtin dataset is generated from the fetched ``iter_jintrac`` sample and contains:
+    - one ``equilibrium`` IDS (copied from ``iter_jintrac``), and
+    - one ``radiation`` IDS with both core-profile emissivity and GGD emissivity values.
+
+    .. note::
+
+        This dataset is intended for testing and demonstration purposes only, and does not
+        represent any real physical scenario.
+
+    Returns
+    -------
+    str
+        Path to the synthetic dataset file.
+
+    Examples
+    --------
+    >>> from cherab.imas import datasets
+    >>> data_path = datasets.iter_jintrac_radiation_values()
+    >>> data_path
+    '.../cherab/imas/iter_jintrac_radiation_values.nc'
+    """
+    path = Path(pooch.os_cache("cherab/imas")) / "iter_jintrac_radiation_values.nc"
+
+    path.parent.mkdir(parents=True, exist_ok=True)
+
+    if not path.exists():
+        equilibrium, radiation = _iter_jintrac_radiation_values_data()
+        with DBEntry(str(path), "w", dd_version=equilibrium._version) as entry:
+            entry.put(equilibrium)
+            entry.put(radiation)
+
+    return str(path)

src/cherab/imas/datasets/_registry.py

@@ -12,5 +12,6 @@
     "iter_jintrac": ["iter_scenario_53298_seq1_DD4.nc", "iter_scenario_53298_seq1_DD4_mod.nc"],
     "iter_solps": ["iter_scenario_123364_1.nc"],
     "iter_jorek": ["iter_disruption_113112_1.nc"],
+    "iter_jintrac_radiation_values": ["iter_jintrac_radiation_values.nc"],
     "bolometer_moc": ["bolometer_moc.nc"],
 }