From 0d0d37ae57fb2716fa3a481f411c1fd0baa0a79c Mon Sep 17 00:00:00 2001
From: Brandon Barker <bbarker.py@protonmail.ch>
Date: Mon, 15 Dec 2025 08:07:15 -0700
Subject: [PATCH] feat: can select kumamoto opacities. Not implemented so this
 prints an error and exits

---
 core/input.c            | 4 ++++
 prob/torus_cbc/build.py | 4 +++-
 script/config.py        | 4 ++++
 3 files changed, 11 insertions(+), 1 deletion(-)

diff --git a/core/input.c b/core/input.c
index 9a121b7..019d5b6 100644
--- a/core/input.c
+++ b/core/input.c
@@ -158,6 +158,10 @@ void set_core_params() {
 #if HDF5_OPACITIES
   set_param("opac_file", &opac_file);
 #endif // HDF5_OPACITIES
+#if KUMAMOTO_OPACITIES
+  printf("Trying to use Kumamoto opacities! These are not currently implemented, crashing out...\n");
+  exit(-1);
+#endif // KUMAMOTO_OPACITIES
 #endif // RADTYPE_NEUTRINOS
 #endif // RADIATION
 
diff --git a/prob/torus_cbc/build.py b/prob/torus_cbc/build.py
index b39e2d5..ae28725 100644
--- a/prob/torus_cbc/build.py
+++ b/prob/torus_cbc/build.py
@@ -40,6 +40,7 @@
 TRACERTEST = '-tracertest' in sys.argv
 RESTARTTEST = '-restarttest' in sys.argv
 HDF = '-hdf' in sys.argv
+KUMAMOTO = '-kumamoto' in sys.argv
 OSCILLATIONS = "-oscillations" in sys.argv
 N1N2N3CPU_FROM_CLI = '-n1n2n3cpu' in sys.argv
 N1N2N3TOT_FROM_CLI = '-n1n2n3tot' in sys.argv
@@ -71,7 +72,7 @@
 EMISS = not NOEMISS
 SCATT = not (NOSCATT or KILL)
 ABS = not (NOABS or KILL)
-FORTRAN = NEUTRINOS and not HDF
+FORTRAN = NEUTRINOS and not HDF and not KUMAMOTO
 TRACERS = not NOTRACE
 
 USE_TABLE = GAMTABLE or RELTABLE
@@ -411,6 +412,7 @@
     bhl.config.set_cparm('KILL_ALL_PACKETS', True)
 bhl.config.set_cparm('BURROWS_OPACITIES', FORTRAN)
 bhl.config.set_cparm('HDF5_OPACITIES', HDF)
+bhl.config.set_cparm('KUMAMOTO_OPACITIES',KUMAMOTO)
 bhl.config.set_cparm('NU_BINS', 61)
 bhl.config.set_cparm('ESTIMATE_THETAE', False)
 bhl.config.set_cparm('GRAYABSORPTION',  False)
diff --git a/script/config.py b/script/config.py
index 892091d..1781ac3 100644
--- a/script/config.py
+++ b/script/config.py
@@ -337,6 +337,10 @@ def build(PROBLEM, PATHS):
       print_config("HDF5_OPACITIES", CPARMS["HDF5_OPACITIES"])
     else:
       set_cparm("HDF5_OPACITIES", 0)
+    if util.parm_is_active(CPARMS, "KUMAMOTO_OPACITIES"):
+      print_config("KUMAMOTO_OPACITIES", CPARMS["KUMAMOTO_OPACITIES"])
+    else:
+      set_cparm("KUMAMOTO_OPACITIES", 0)
     if util.parm_is_active(CPARMS, "RAD_NUM_TYPES"):
       print_config("RAD_NUM_TYPES", CPARMS["RAD_NUM_TYPES"])
     else: