sed-model 0.3.0__tar.gz → 0.3.2__tar.gz

{sed_model-0.3.0 → sed_model-0.3.2}/.github/workflows/build.yml

@@ -31,16 +31,16 @@ jobs:
   # ---------------------------------------------------------------------------
   sdist:
     name: Source distribution
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     outputs:
       version: ${{ steps.version.outputs.version }}
       is_new: ${{ steps.pypi.outputs.is_new }}
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
         with:
           fetch-depth: 0   # meson dist builds from git
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: "3.12"
 
@@ -75,7 +75,7 @@ jobs:
           python -c "import sed_model, sed_model.cc_api"
           python -m pytest tests/test_kernels.py -p no:cacheprovider
 
-      - uses: actions/upload-artifact@v4
+      - uses: actions/upload-artifact@v7
         with:
           name: sdist
           path: dist/*.tar.gz
@@ -86,13 +86,13 @@ jobs:
   # ---------------------------------------------------------------------------
   pure_wheel:
     name: Pure-Python wheel
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
         with:
           fetch-depth: 0
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: "3.12"
 
@@ -112,7 +112,7 @@ jobs:
           print("sed_model", sed_model.__version__, "- pure wheel OK")
           PY
 
-      - uses: actions/upload-artifact@v4
+      - uses: actions/upload-artifact@v7
         with:
           name: pure-wheel
           path: dist/*.whl
@@ -129,7 +129,7 @@ jobs:
       fail-fast: false
       matrix:
         include:
-          - { os: ubuntu-latest,    label: linux-x86_64 }
+          - { os: ubuntu-24.04,     label: linux-x86_64 }
           - { os: ubuntu-24.04-arm, label: linux-aarch64 }
           # The Fortran runtime is linked statically on macOS (see meson.build),
           # so the wheel's minimum macOS is this target, not whatever the
@@ -138,11 +138,11 @@ jobs:
           - { os: macos-14,         label: macos-arm64,  macos_target: "14.0" }
           - { os: macos-15-intel,   label: macos-x86_64, macos_target: "15.0" }
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
         with:
           fetch-depth: 0
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: "3.12"
 
@@ -157,7 +157,7 @@ jobs:
         env:
           CIBW_ENVIRONMENT_MACOS: MACOSX_DEPLOYMENT_TARGET=${{ matrix.macos_target }}
 
-      - uses: actions/upload-artifact@v4
+      - uses: actions/upload-artifact@v7
         with:
           name: wheels-${{ matrix.label }}
           path: wheelhouse/*.whl
@@ -173,12 +173,12 @@ jobs:
       fail-fast: false
       matrix:
         include:
-          - { os: ubuntu-latest,    label: linux-x86_64,  image: quay.io/pypa/manylinux_2_28_x86_64 }
+          - { os: ubuntu-24.04,     label: linux-x86_64,  image: quay.io/pypa/manylinux_2_28_x86_64 }
           - { os: ubuntu-24.04-arm, label: linux-aarch64, image: quay.io/pypa/manylinux_2_28_aarch64 }
           - { os: macos-14,         label: macos-arm64,   macos_target: "14.0" }
           - { os: macos-15-intel,   label: macos-x86_64,  macos_target: "15.0" }
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
         with:
           fetch-depth: 0
 
@@ -201,7 +201,7 @@ jobs:
           version=$(git describe --tags --always)
           tar -C stage -czf "sed_model-fortran-${version}-${{ matrix.label }}.tar.gz" sed_model-fortran
 
-      - uses: actions/upload-artifact@v4
+      - uses: actions/upload-artifact@v7
         with:
           name: fortran-${{ matrix.label }}
           path: sed_model-fortran-*.tar.gz
@@ -226,7 +226,7 @@ jobs:
       && (github.ref == 'refs/heads/main' || startsWith(github.ref, 'refs/tags/v'))
       && needs.sdist.outputs.is_new == 'true'
     needs: [sdist, pure_wheel, wheels, fortran_lib]
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     environment: pypi
     concurrency:
       group: publish
@@ -234,11 +234,11 @@ jobs:
     permissions:
       id-token: write   # trusted publishing
     steps:
-      - uses: actions/download-artifact@v4
+      - uses: actions/download-artifact@v8
         with:
           path: artifacts
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: "3.12"
 
@@ -292,12 +292,12 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        os: [ubuntu-latest, ubuntu-24.04-arm, macos-14, macos-15-intel]
+        os: [ubuntu-24.04, ubuntu-24.04-arm, macos-14, macos-15-intel]
         python: ["3.10", "3.14"]
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: ${{ matrix.python }}
 
@@ -328,11 +328,11 @@ jobs:
   release:
     name: GitHub release
     needs: [sdist, verify_pypi]
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     permissions:
       contents: write
     steps:
-      - uses: actions/download-artifact@v4
+      - uses: actions/download-artifact@v8
         with:
           path: artifacts
 
@@ -345,11 +345,14 @@ jobs:
           fi
           # An existing tag must already point at the commit that was just built;
           # otherwise the release would attach these artifacts to different code.
-          tag_sha=$(gh api "repos/$GITHUB_REPOSITORY/commits/$tag" --jq .sha 2>/dev/null || true)
-          if [ -n "$tag_sha" ] && [ "$tag_sha" != "$GITHUB_SHA" ]; then
-            echo "::error::Tag $tag already exists at $tag_sha, but this build is $GITHUB_SHA." \
-                 "Move or delete the tag (or bump the version) and re-run."
-            exit 1
+          # gh prints the API error body to stdout on failure, so test the exit
+          # status rather than the output: a missing tag means nothing to compare.
+          if tag_sha=$(gh api "repos/$GITHUB_REPOSITORY/commits/$tag" --jq .sha 2>/dev/null); then
+            if [ "$tag_sha" != "$GITHUB_SHA" ]; then
+              echo "::error::Tag $tag already exists at $tag_sha, but this build is $GITHUB_SHA." \
+                   "Move or delete the tag (or bump the version) and re-run."
+              exit 1
+            fi
           fi
           gh release create "$tag" --repo "$GITHUB_REPOSITORY" --target "$GITHUB_SHA" \
             --title "$tag" --generate-notes \

{sed_model-0.3.0 → sed_model-0.3.2}/.github/workflows/docs.yml

@@ -16,11 +16,11 @@ concurrency:
 
 jobs:
   build:
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v7
 
-      - uses: actions/setup-python@v5
+      - uses: actions/setup-python@v7
         with:
           python-version: "3.11"
 
@@ -34,16 +34,16 @@ jobs:
       - name: Build site
         run: mkdocs build --strict
 
-      - uses: actions/upload-pages-artifact@v3
+      - uses: actions/upload-pages-artifact@v5
         with:
           path: site
 
   deploy:
     needs: build
-    runs-on: ubuntu-latest
+    runs-on: ubuntu-24.04
     environment:
       name: github-pages
       url: ${{ steps.deployment.outputs.page_url }}
     steps:
       - id: deployment
-        uses: actions/deploy-pages@v4
+        uses: actions/deploy-pages@v5

{sed_model-0.3.0 → sed_model-0.3.2}/CHANGELOG.md

@@ -1,5 +1,27 @@
 # Changelog
 
+## [0.3.2] — 2026-10-25
+
+- **Filter files from SED_Tools 0.4.0 load again.** SED_Tools 0.4.0 puts its
+  metadata on the column row (`Wavelength,Transmission,FilterID=...;DetectorType=...`);
+  the loader now reads only the first two columns instead of failing with
+  `expected two-column file`. Required for SED_Tools >= 0.4.0 filter files.
+  The metadata itself is not read yet.
+- Tests find their data in `~/.sed_tools/data` as well as the older locations
+  and `$MESA_DIR/data/colors_data`.
+
+## [0.3.1] — 2026-10-23
+
+Packaging and CI only; no numerical changes.
+
+- README logos use absolute URLs, so they render on the PyPI project page.
+- Square `docs/favicon.png` for the documentation site.
+- Release job: the tag guard no longer misreads the API error body as a commit
+  when the tag does not yet exist.
+- CI: actions moved to current majors (Node 24), and runners pinned to
+  `ubuntu-24.04` instead of `ubuntu-latest`.
+- Removed the legacy, uncompiled `fortran/private/cc_api.f90`.
+
 ## [0.3.0] — 2026-10-21
 
 Fortran kernels synchronised with the current MESA `colors` module, and

{sed_model-0.3.0 → sed_model-0.3.2}/CITATION.cff

@@ -7,7 +7,7 @@ authors:
   - family-names: "Miller"
     given-names: "Niall"
     email: "niall.j.miller@gmail.com"
-version: "0.3.0"
+version: "0.3.2"
 repository-code: "https://github.com/nialljmiller/SED_Model"
 url: "https://nialljmiller.github.io/SED_Model/"
 license: MIT

{sed_model-0.3.0 → sed_model-0.3.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: sed-model
-Version: 0.3.0
+Version: 0.3.2
 Summary: Forward and inverse stellar photometry from SED grids.
 Keywords: astronomy,astrophysics,stellar astronomy,synthetic photometry,spectral energy distributions,sed,stellar parameters,bayesian inference,mcmc
 Author-Email: Niall Miller <niall.j.miller@gmail.com>
@@ -27,7 +27,7 @@ Requires-Dist: matplotlib; extra == "demo"
 Description-Content-Type: text/markdown
 
 <p align="center">
-  <img src="docs/SED_Model_logo.jpeg" alt="SED_Model logo" width="600"/>
+  <img src="https://raw.githubusercontent.com/nialljmiller/SED_Model/main/docs/SED_Model_logo.jpeg" alt="SED_Model logo" width="600"/>
 </p>
 
 <h1 align="center">SED_Model</h1>
@@ -36,7 +36,7 @@ Description-Content-Type: text/markdown
 
 SED_Model is a Python package for computing observer-ready synthetic spectral energy distributions (SEDs) and broadband magnitudes from stellar atmosphere grids, and for recovering stellar parameters from observed photometry via Bayesian MCMC inference. It is designed to work directly downstream of `SED_Tools` and as a standalone drop-in for workflows that use MESA-style atmosphere grids and filter transmission curves.
 
-Version: `0.3.0`
+Version: `0.3.2`
 
 **Documentation:** https://nialljmiller.github.io/SED_Model/
 
@@ -461,5 +461,5 @@ SED_Model uses a mixed-license structure. The Python source code, documentation,
 ---
 
 <p align="center">
-  <img src="docs/SED_Model_Logo_full.png" alt="SED_Model: The Pythonic Inside Story of MCMC-Based Stellar Parameter Inference from Photometric Grids, by Niall Miller" width="360"/>
+  <img src="https://raw.githubusercontent.com/nialljmiller/SED_Model/main/docs/SED_Model_Logo_full.png" alt="SED_Model: The Pythonic Inside Story of MCMC-Based Stellar Parameter Inference from Photometric Grids, by Niall Miller" width="360"/>
 </p>

{sed_model-0.3.0 → sed_model-0.3.2}/README.md

@@ -1,5 +1,5 @@
 <p align="center">
-  <img src="docs/SED_Model_logo.jpeg" alt="SED_Model logo" width="600"/>
+  <img src="https://raw.githubusercontent.com/nialljmiller/SED_Model/main/docs/SED_Model_logo.jpeg" alt="SED_Model logo" width="600"/>
 </p>
 
 <h1 align="center">SED_Model</h1>
@@ -8,7 +8,7 @@
 
 SED_Model is a Python package for computing observer-ready synthetic spectral energy distributions (SEDs) and broadband magnitudes from stellar atmosphere grids, and for recovering stellar parameters from observed photometry via Bayesian MCMC inference. It is designed to work directly downstream of `SED_Tools` and as a standalone drop-in for workflows that use MESA-style atmosphere grids and filter transmission curves.
 
-Version: `0.3.0`
+Version: `0.3.2`
 
 **Documentation:** https://nialljmiller.github.io/SED_Model/
 
@@ -433,5 +433,5 @@ SED_Model uses a mixed-license structure. The Python source code, documentation,
 ---
 
 <p align="center">
-  <img src="docs/SED_Model_Logo_full.png" alt="SED_Model: The Pythonic Inside Story of MCMC-Based Stellar Parameter Inference from Photometric Grids, by Niall Miller" width="360"/>
+  <img src="https://raw.githubusercontent.com/nialljmiller/SED_Model/main/docs/SED_Model_Logo_full.png" alt="SED_Model: The Pythonic Inside Story of MCMC-Based Stellar Parameter Inference from Photometric Grids, by Niall Miller" width="360"/>
 </p>

{sed_model-0.3.0 → sed_model-0.3.2}/docs/index.md

@@ -8,7 +8,7 @@
 
 SED_Model is a Python package for computing observer-ready synthetic spectral energy distributions (SEDs) and broadband magnitudes from stellar atmosphere grids, and for recovering stellar parameters from observed photometry via Bayesian MCMC inference. It works directly downstream of [SED_Tools](https://github.com/nialljmiller/SED_Tools) and as a standalone drop-in for workflows that use MESA-style atmosphere grids and filter transmission curves.
 
-Current version: **0.3.0** ([PyPI](https://pypi.org/project/sed-model/))
+Current version: **0.3.2** ([PyPI](https://pypi.org/project/sed-model/))
 
 ## Two paired models, one parameter language
 

{sed_model-0.3.0 → sed_model-0.3.2}/docs/testing.md

@@ -9,13 +9,13 @@ Both suites skip gracefully — with an explanatory reason — when the Fortran
 
 ## Data locations
 
-The tests look for SED_Tools data products at `~/SED_Tools/data/...` by default, overridable via environment variables:
+The tests look for SED_Tools data products in `~/.sed_tools/data`, then `~/SED_Tools/data`, then `$MESA_DIR/data/colors_data`, using the first location that exists. Each path can be overridden with an environment variable:
 
-| Variable | Default |
+| Variable | Path under the data root |
 |---|---|
-| `CC_KURUCZ_DIR` | `~/SED_Tools/data/stellar_models/Kurucz2003all` |
-| `CC_JOHNSON_DIR` | `~/SED_Tools/data/filters/Generic/Johnson` |
-| `CC_VEGA_SED` | `~/SED_Tools/data/stellar_models/vega_flam.csv` |
+| `CC_KURUCZ_DIR` | `stellar_models/Kurucz2003all` |
+| `CC_JOHNSON_DIR` | `filters/Generic/Johnson` |
+| `CC_VEGA_SED` | `stellar_models/vega_flam.csv` |
 
 ## Forward model: validated against MESA
 

{sed_model-0.3.0 → sed_model-0.3.2}/meson.build

@@ -2,7 +2,7 @@
 # and for the others the checks below give a readable error if one is missing.
 project(
   'sed_model',
-  version: '0.3.0',
+  version: '0.3.2',
   meson_version: '>=1.1.0',
 )
 

{sed_model-0.3.0 → sed_model-0.3.2}/mkdocs.yml

@@ -8,7 +8,7 @@ edit_uri: edit/main/docs/
 theme:
   name: material
   logo: SED_Model_logo.jpeg
-  favicon: SED_Model_logo.jpeg
+  favicon: favicon.png
   palette:
     - media: "(prefers-color-scheme: light)"
       scheme: default

{sed_model-0.3.0 → sed_model-0.3.2}/sed_model/filters.py

@@ -227,8 +227,11 @@ def _load_filter_dat(path: Path) -> tuple[np.ndarray, np.ndarray]:
         raise FileNotFoundError(f"Filter file not found: {path}")
 
     delimiter = _sniff_delimiter(path)
+    # usecols: SED_Tools >= 0.4.0 appends a metadata field to the column row
+    # (``Wavelength,Transmission,FilterID=...;DetectorType=...``), which would
+    # otherwise make genfromtxt expect three columns.
     data = np.genfromtxt(path, comments="#", delimiter=delimiter,
-                         invalid_raise=False)
+                         usecols=(0, 1), invalid_raise=False)
 
     if data.ndim == 2:
         mask = np.isfinite(data[:, 0]) & np.isfinite(data[:, 1])

sed_model-0.3.2/tests/_data_paths.py

@@ -0,0 +1,29 @@
+"""Locations of the SED_Tools data products the tests use.
+
+Each path can be set with an environment variable; otherwise the first
+existing candidate is used (SED_Tools' own cache, the older ``~/SED_Tools``
+layout, then the data shipped with a MESA checkout).
+"""
+
+from __future__ import annotations
+
+import os
+
+_ROOTS = [
+    os.path.expanduser("~/.sed_tools/data"),
+    os.path.expanduser("~/SED_Tools/data"),
+]
+if os.environ.get("MESA_DIR"):
+    _ROOTS.append(os.path.join(os.environ["MESA_DIR"], "data", "colors_data"))
+
+
+def _find(env: str, *parts: str) -> str:
+    if os.environ.get(env):
+        return os.environ[env]
+    candidates = [os.path.join(root, *parts) for root in _ROOTS]
+    return next((c for c in candidates if os.path.exists(c)), candidates[0])
+
+
+KURUCZ_DIR  = _find("CC_KURUCZ_DIR", "stellar_models", "Kurucz2003all")
+JOHNSON_DIR = _find("CC_JOHNSON_DIR", "filters", "Generic", "Johnson")
+VEGA_SED    = _find("CC_VEGA_SED", "stellar_models", "vega_flam.csv")

{sed_model-0.3.0 → sed_model-0.3.2}/tests/test_forward.py

@@ -33,18 +33,7 @@ import numpy as np
 # ---------------------------------------------------------------------------
 # Path configuration
 # ---------------------------------------------------------------------------
-KURUCZ_DIR  = os.environ.get(
-    "CC_KURUCZ_DIR",
-    os.path.expanduser("~/SED_Tools/data/stellar_models/Kurucz2003all"),
-)
-JOHNSON_DIR = os.environ.get(
-    "CC_JOHNSON_DIR",
-    os.path.expanduser("~/SED_Tools/data/filters/Generic/Johnson"),
-)
-VEGA_SED    = os.environ.get(
-    "CC_VEGA_SED",
-    os.path.expanduser("~/SED_Tools/data/stellar_models/vega_flam.csv"),
-)
+from _data_paths import KURUCZ_DIR, JOHNSON_DIR, VEGA_SED
 
 R_SUN  = 6.957e10
 D_10PC = 3.0857e19

{sed_model-0.3.0 → sed_model-0.3.2}/tests/test_inverse.py

@@ -34,18 +34,7 @@ import pytest
 # ---------------------------------------------------------------------------
 # Paths
 # ---------------------------------------------------------------------------
-KURUCZ_DIR  = os.environ.get(
-    "CC_KURUCZ_DIR",
-    os.path.expanduser("~/SED_Tools/data/stellar_models/Kurucz2003all"),
-)
-JOHNSON_DIR = os.environ.get(
-    "CC_JOHNSON_DIR",
-    os.path.expanduser("~/SED_Tools/data/filters/Generic/Johnson"),
-)
-VEGA_SED    = os.environ.get(
-    "CC_VEGA_SED",
-    os.path.expanduser("~/SED_Tools/data/stellar_models/vega_flam.csv"),
-)
+from _data_paths import KURUCZ_DIR, JOHNSON_DIR, VEGA_SED
 
 R_SUN  = 6.957e10
 D_10PC = 3.0857e19

{sed_model-0.3.0 → sed_model-0.3.2}/tests/test_kernels.py

@@ -142,3 +142,26 @@ def test_linear_interp_handles_nan_cube_node():
     flux, ierr = cc.interp_sed_linear(0.4, 0.0, 0.0, x, one, one, cube)
     assert np.all(np.isfinite(flux))
     assert ierr >= 1
+
+
+def test_filter_header_row_with_metadata_field(tmp_path):
+    """SED_Tools >= 0.4.0 appends ``;``-separated metadata as a third comma
+    field of the column row.  The loader reads only the two data columns, for
+    that layout, the bare row, and the older leading ``# Key: value`` lines."""
+    from sed_model.filters import load_filters
+
+    rows = "3200.0,2.4e-08\n3300.0,0.5\n3400.0,0.9\n3500.0,0.0\n"
+    layouts = {
+        "bare": "Wavelength,Transmission\n" + rows,
+        "meta": ("Wavelength,Transmission,FilterID=GAIA/GAIA3.G;DetectorType=photon;"
+                 "DetectorTypeSource=svo-curve-match;WavelengthUnit=Angstrom;"
+                 "Components=filter+qe+instrument\n" + rows),
+        "comments": "# FilterID: GAIA/GAIA3.G\n# DetectorType: photon\nWavelength,Transmission\n" + rows,
+    }
+    waves = {}
+    for name, text in layouts.items():
+        (tmp_path / f"{name}.dat").write_text(text)
+        (f,) = load_filters([tmp_path / f"{name}.dat"])
+        waves[name] = (f.wavelengths.tolist(), f.transmission.tolist())
+    assert waves["meta"] == waves["bare"] == waves["comments"]
+    assert waves["meta"][0] == [3200.0, 3300.0, 3400.0, 3500.0]

sed_model-0.3.0/fortran/private/cc_api.f90

@@ -1,271 +0,0 @@
-! ***********************************************************************
-! fortran/private/cc_api.f90  —  NOT COMPILED
-!
-! This file is a legacy reference copy that routes through colors_lib.f90
-! (the public aggregation module).  It is NOT included in the Makefile or
-! meson.build build sources and is therefore never compiled.
-!
-! The file that IS compiled is fortran/cc_api.f90 (one level up), which
-! uses the private kernels directly.  That workaround exists because the
-! pinned f2py cannot wrap a module that only re-exports procedures from
-! other modules (see fortran/cc_api.f90 header and MIGRATION.md).
-!
-! This copy is kept for reference only.  All changes to the compiled API
-! must go into fortran/cc_api.f90, NOT this file.
-! ***********************************************************************
-
-module cc_api
-   use colors_lib, only: dp, &
-      hermite_interp_vector, &
-      trilinear_interp_vector, &
-      lib_dilute_flux    => dilute_flux, &
-      lib_trapz          => trapezoidal_integration, &
-      lib_simpson        => simpson_integration, &
-      lib_interp_filter  => interp_filter_onto_sed, &
-      lib_synthetic_flux => calculate_synthetic_flux, &
-      lib_magnitude      => magnitude, &
-      lib_bol_phot       => calculate_bolometric_phot, &
-      lib_vega_zp        => compute_vega_zero_point, &
-      lib_ab_zp          => compute_ab_zero_point, &
-      lib_st_zp          => compute_st_zero_point
-
-   implicit none
-   private
-
-   public :: interp_sed_hermite
-   public :: interp_sed_linear
-   public :: dilute_flux
-   public :: synthetic_magnitude
-   public :: bolometric
-   public :: vega_zero_point
-   public :: ab_zero_point
-   public :: st_zero_point
-   public :: trapz
-   public :: simpson
-
-contains
-
-   ! -------------------------------------------------------------------------
-   ! interp_sed_hermite
-   !
-   ! Interpolate a full SED at (teff, logg, meta) using cubic Hermite
-   ! tensor interpolation on the preloaded flux cube. Forwards straight
-   ! to hermite_interp_vector in fortran/private/hermite_interp.f90.
-   ! -------------------------------------------------------------------------
-   subroutine interp_sed_hermite(teff, logg, meta, &
-                                  teff_grid, nt, &
-                                  logg_grid, nl, &
-                                  meta_grid, nm, &
-                                  flux_cube, nw, &
-                                  result_flux, ierr)
-      !f2py intent(in)  :: teff, logg, meta
-      !f2py intent(in)  :: teff_grid, logg_grid, meta_grid
-      !f2py intent(in)  :: nt, nl, nm, nw
-      !f2py intent(in)  :: flux_cube
-      !f2py intent(out) :: result_flux
-      !f2py intent(out) :: ierr
-      real(dp), intent(in)  :: teff, logg, meta
-      integer,  intent(in)  :: nt, nl, nm, nw
-      real(dp), intent(in)  :: teff_grid(nt), logg_grid(nl), meta_grid(nm)
-      real(dp), intent(in)  :: flux_cube(nt, nl, nm, nw)
-      real(dp), intent(out) :: result_flux(nw)
-      integer,  intent(out) :: ierr
-
-      call hermite_interp_vector(teff, logg, meta, &
-                                 teff_grid, logg_grid, meta_grid, &
-                                 flux_cube, nw, result_flux, ierr)
-   end subroutine interp_sed_hermite
-
-
-   ! -------------------------------------------------------------------------
-   ! interp_sed_linear
-   !
-   ! Same interface as interp_sed_hermite but using trilinear interpolation.
-   ! -------------------------------------------------------------------------
-   subroutine interp_sed_linear(teff, logg, meta, &
-                                 teff_grid, nt, &
-                                 logg_grid, nl, &
-                                 meta_grid, nm, &
-                                 flux_cube, nw, &
-                                 result_flux, ierr)
-      !f2py intent(in)  :: teff, logg, meta
-      !f2py intent(in)  :: teff_grid, logg_grid, meta_grid
-      !f2py intent(in)  :: nt, nl, nm, nw
-      !f2py intent(in)  :: flux_cube
-      !f2py intent(out) :: result_flux
-      !f2py intent(out) :: ierr
-      real(dp), intent(in)  :: teff, logg, meta
-      integer,  intent(in)  :: nt, nl, nm, nw
-      real(dp), intent(in)  :: teff_grid(nt), logg_grid(nl), meta_grid(nm)
-      real(dp), intent(in)  :: flux_cube(nt, nl, nm, nw)
-      real(dp), intent(out) :: result_flux(nw)
-      integer,  intent(out) :: ierr
-
-      call trilinear_interp_vector(teff, logg, meta, &
-                                   teff_grid, logg_grid, meta_grid, &
-                                   flux_cube, nw, result_flux, ierr)
-   end subroutine interp_sed_linear
-
-
-   ! -------------------------------------------------------------------------
-   ! dilute_flux
-   !
-   ! Apply (R/d)^2 dilution to convert surface flux to observed flux.
-   ! -------------------------------------------------------------------------
-   subroutine dilute_flux(surface_flux, nw, R, d, observed_flux)
-      !f2py intent(in)  :: surface_flux, nw, R, d
-      !f2py intent(out) :: observed_flux
-      integer,  intent(in)  :: nw
-      real(dp), intent(in)  :: surface_flux(nw), R, d
-      real(dp), intent(out) :: observed_flux(nw)
-      call lib_dilute_flux(surface_flux, R, d, observed_flux)
-   end subroutine dilute_flux
-
-
-   ! -------------------------------------------------------------------------
-   ! synthetic_magnitude
-   !
-   ! Compute a synthetic magnitude in a single filter from a diluted SED.
-   !
-   ! Steps performed internally:
-   !   1. Interpolate filter transmission onto the SED wavelength grid
-   !   2. Photon-counting in-band flux integration
-   !   3. m = -2.5 log10(F_band / F_zp)
-   !
-   ! Parameters
-   ! ----------
-   ! sed_wave(nw)        : SED wavelength grid (A)
-   ! obs_flux(nw)        : diluted (observer-frame) SED flux (erg/s/cm^2/A)
-   ! filt_wave(nf)       : filter wavelength grid (A)
-   ! filt_trans(nf)      : filter transmission [0,1]
-   ! zero_point          : precomputed photometric zero-point
-   ! mag                 : output magnitude
-   ! band_flux           : output in-band flux (before zero-point)
-   ! ierr                : 0=ok, 1=integration failure, 2=non-positive flux
-   ! -------------------------------------------------------------------------
-   subroutine synthetic_magnitude(sed_wave, obs_flux, nw, &
-                                   filt_wave, filt_trans, nf, &
-                                   zero_point, &
-                                   mag, band_flux, ierr)
-      !f2py intent(in)  :: sed_wave, obs_flux, nw
-      !f2py intent(in)  :: filt_wave, filt_trans, nf
-      !f2py intent(in)  :: zero_point
-      !f2py intent(out) :: mag, band_flux, ierr
-      integer,  intent(in)  :: nw, nf
-      real(dp), intent(in)  :: sed_wave(nw), obs_flux(nw)
-      real(dp), intent(in)  :: filt_wave(nf), filt_trans(nf)
-      real(dp), intent(in)  :: zero_point
-      real(dp), intent(out) :: mag, band_flux
-      integer,  intent(out) :: ierr
-
-      real(dp) :: filt_on_sed(nw)
-      integer  :: ierr2
-
-      call lib_interp_filter(filt_wave, filt_trans, sed_wave, filt_on_sed, ierr)
-      if (ierr /= 0) then
-         mag = -99.9_dp; band_flux = -1.0_dp; return
-      end if
-
-      call lib_synthetic_flux(sed_wave, obs_flux * filt_on_sed, filt_on_sed, band_flux)
-      if (band_flux <= 0.0_dp) then
-         mag = -99.9_dp; ierr = 1; return
-      end if
-
-      call lib_magnitude(band_flux, zero_point, mag, ierr2)
-      if (ierr2 /= 0) ierr = ierr2
-   end subroutine synthetic_magnitude
-
-
-   ! -------------------------------------------------------------------------
-   ! bolometric
-   !
-   ! Compute bolometric flux and magnitude from a diluted SED.
-   ! -------------------------------------------------------------------------
-   subroutine bolometric(sed_wave, obs_flux, nw, bol_flux, bol_mag, ierr)
-      !f2py intent(in)  :: sed_wave, obs_flux, nw
-      !f2py intent(out) :: bol_flux, bol_mag, ierr
-      integer,  intent(in)  :: nw
-      real(dp), intent(in)  :: sed_wave(nw), obs_flux(nw)
-      real(dp), intent(out) :: bol_flux, bol_mag
-      integer,  intent(out) :: ierr
-
-      call lib_bol_phot(sed_wave, obs_flux, bol_mag, bol_flux, ierr)
-   end subroutine bolometric
-
-
-   ! -------------------------------------------------------------------------
-   ! Zero-point subroutines (called once at init from Python)
-   !
-   ! colors_lib's compute_*_zero_point functions return the MESA-style
-   ! function value with -1.0 as the failure sentinel; ierr here is
-   ! derived from that sentinel to keep the existing ierr-based Python
-   ! contract unchanged.
-   ! -------------------------------------------------------------------------
-
-   subroutine vega_zero_point(vega_wave, vega_flux, nv, &
-                               filt_wave, filt_trans, nf, &
-                               zp, ierr)
-      !f2py intent(in)  :: vega_wave, vega_flux, nv
-      !f2py intent(in)  :: filt_wave, filt_trans, nf
-      !f2py intent(out) :: zp, ierr
-      integer,  intent(in)  :: nv, nf
-      real(dp), intent(in)  :: vega_wave(nv), vega_flux(nv)
-      real(dp), intent(in)  :: filt_wave(nf), filt_trans(nf)
-      real(dp), intent(out) :: zp
-      integer,  intent(out) :: ierr
-
-      zp = lib_vega_zp(vega_wave, vega_flux, filt_wave, filt_trans)
-      ierr = merge(0, 1, zp > 0.0_dp)
-   end subroutine vega_zero_point
-
-
-   subroutine ab_zero_point(filt_wave, filt_trans, nf, zp, ierr)
-      !f2py intent(in)  :: filt_wave, filt_trans, nf
-      !f2py intent(out) :: zp, ierr
-      integer,  intent(in)  :: nf
-      real(dp), intent(in)  :: filt_wave(nf), filt_trans(nf)
-      real(dp), intent(out) :: zp
-      integer,  intent(out) :: ierr
-
-      zp = lib_ab_zp(filt_wave, filt_trans)
-      ierr = merge(0, 1, zp > 0.0_dp)
-   end subroutine ab_zero_point
-
-
-   subroutine st_zero_point(filt_wave, filt_trans, nf, zp, ierr)
-      !f2py intent(in)  :: filt_wave, filt_trans, nf
-      !f2py intent(out) :: zp, ierr
-      integer,  intent(in)  :: nf
-      real(dp), intent(in)  :: filt_wave(nf), filt_trans(nf)
-      real(dp), intent(out) :: zp
-      integer,  intent(out) :: ierr
-
-      zp = lib_st_zp(filt_wave, filt_trans)
-      ierr = merge(0, 1, zp > 0.0_dp)
-   end subroutine st_zero_point
-
-
-   ! -------------------------------------------------------------------------
-   ! Standalone integration (exposed for testing / Python use)
-   ! -------------------------------------------------------------------------
-
-   subroutine trapz(x, y, n, result)
-      !f2py intent(in)  :: x, y, n
-      !f2py intent(out) :: result
-      integer,  intent(in)  :: n
-      real(dp), intent(in)  :: x(n), y(n)
-      real(dp), intent(out) :: result
-      call lib_trapz(x, y, result)
-   end subroutine trapz
-
-
-   subroutine simpson(x, y, n, result)
-      !f2py intent(in)  :: x, y, n
-      !f2py intent(out) :: result
-      integer,  intent(in)  :: n
-      real(dp), intent(in)  :: x(n), y(n)
-      real(dp), intent(out) :: result
-      call lib_simpson(x, y, result)
-   end subroutine simpson
-
-end module cc_api