Fix CCE 19.0.0 + OpenACC: patch pyrometheus thermochem GPU_ROUTINE macro

Pyrometheus generates !DIR$ INLINEALWAYS for Cray+OpenACC but omits
!acc routine seq, so thermochem routines are not registered as OpenACC
device routines. Calling them from GPU kernels causes runtime memory
access faults at time step 1.

Post-process the generated m_thermochem.f90 to replace the broken
Cray-specific #ifdef block with a plain #define GPU_ROUTINE(name)
!acc routine seq, matching the fix already applied to parallel_macros.fpp
in commit 66ad9b2.

Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>

Author: Spencer Bryngelson

toolchain/mfc/run/input.py

@@ -90,13 +90,25 @@ def generate_fpp(self, target) -> None:
             directive_str = None
 
         # Write the generated Fortran code to the m_thermochem.f90 file with the chosen precision
+        thermochem_code = pyro.FortranCodeGenerator().generate(
+            "m_thermochem",
+            self.get_cantera_solution(),
+            pyro.CodeGenerationOptions(scalar_type = real_type, directive_offload = directive_str)
+        )
+
+        # CCE 19.0.0 workaround: pyrometheus generates !DIR$ INLINEALWAYS for Cray+ACC
+        # but omits !$acc routine seq, so thermochem routines are not registered as
+        # OpenACC device routines. Replace with plain !$acc routine seq (no INLINEALWAYS).
+        if directive_str == 'acc':
+            thermochem_code = thermochem_code.replace(
+                "#ifdef _CRAYFTN\n#define GPU_ROUTINE(name) !DIR$ INLINEALWAYS name\n"
+                "#else\n#define GPU_ROUTINE(name) !$acc routine seq\n#endif",
+                "#define GPU_ROUTINE(name) !$acc routine seq"
+            )
+
         common.file_write(
             os.path.join(modules_dir, "m_thermochem.f90"),
-            pyro.FortranCodeGenerator().generate(
-                "m_thermochem",
-                self.get_cantera_solution(),
-                pyro.CodeGenerationOptions(scalar_type = real_type, directive_offload = directive_str)
-            ),
+            thermochem_code,
             True
         )