Regression - MINOR - Drop the gal_flag provenance pin now that it is retired

The Galerkin cases pinned both integrator and gal_flag because the Galerkin
Delta-prime used to be the additive path -- gal_flag layered onto a riccati
sweep -- so the pair was what named the formalism. The standalone galerkin
integrator has since landed and gal_flag is gone from the control struct, so the
decks now declare integrator = "galerkin" and that pin alone identifies the
formalism. The gal_flag pin would only ever read "<unset>".

Also reattaches two section comments that the original insertion orphaned onto
the provenance block.

Co-Authored-By: Claude Opus 5 <noreply@anthropic.com>

Author: d-burg

regression-harness/cases/gal_resistive_diiid.toml

@@ -9,13 +9,10 @@ name = "gal_resistive_diiid"
 description = "DIII-D-like, n=1, RDCON outer-region Galerkin Δ′ with rpec coil columns (delta_coil)"
 example_dir = "examples/DIIID-like_gal_resistive_example"
 
-# Number of resonant surfaces in the gal domain
-# Which formalism produced the Galerkin Δ′ below. Both this deck's `integrator` and its
-# `gal_flag` are pinned because the Galerkin Δ′ is currently the additive path (gal_flag
-# layered on a riccati sweep) rather than a standalone integrator, so the pair is what names
-# the formalism. Once the standalone `integrator = "galerkin"` lands and gal_flag is retired,
-# these pins fail loudly — which is the intent: a deliberate formalism change must regenerate
-# the gold rather than silently re-pointing it at different physics.
+# Which formalism produced the Galerkin Δ′ below. The Δ′ matrix is written to one canonical
+# SingularSurfaces/ path whichever integrator computed it, so the numbers alone cannot say whether
+# a gold pinned from the Galerkin solve is being compared against a BVP run. Pinning the deck's
+# declared choice makes a formalism swap a loud failure instead of a silent re-gold.
 [quantities.integrator]
 h5path = "Input/gpec_toml_raw"
 type = "token"
@@ -24,14 +21,7 @@ label = "integrator (Δ′ provenance)"
 noise_threshold = 0
 order = 1
 
-[quantities.gal_flag]
-h5path = "Input/gpec_toml_raw"
-type = "token"
-extract = "toml_key:ForceFreeStates.gal_flag"
-label = "gal_flag (Δ′ provenance)"
-noise_threshold = 0
-order = 2
-
+# Number of resonant surfaces in the gal domain
 [quantities.gal_msing]
 h5path = "ForceFreeStates/Solutions/GalerkinIntegration/rational_count"
 type = "int_scalar"

regression-harness/cases/gal_resistive_pe.toml

@@ -13,13 +13,10 @@ example_dir = "examples/DIIID-like_gal_resistive_pe_example"
 # FINITE, per-rational-surface singular-coupling scalars (the physics-meaningful driven quantities) — not
 # any whole-profile ξ/ξ′ norm (those diverge at the rationals and are ill-posed as benchmarks).
 
-# Surface locations / count (sanity)
-# Which formalism produced the Galerkin Δ′ below. Both this deck's `integrator` and its
-# `gal_flag` are pinned because the Galerkin Δ′ is currently the additive path (gal_flag
-# layered on a riccati sweep) rather than a standalone integrator, so the pair is what names
-# the formalism. Once the standalone `integrator = "galerkin"` lands and gal_flag is retired,
-# these pins fail loudly — which is the intent: a deliberate formalism change must regenerate
-# the gold rather than silently re-pointing it at different physics.
+# Which formalism produced the Galerkin Δ′ below. The Δ′ matrix is written to one canonical
+# SingularSurfaces/ path whichever integrator computed it, so the numbers alone cannot say whether
+# a gold pinned from the Galerkin solve is being compared against a BVP run. Pinning the deck's
+# declared choice makes a formalism swap a loud failure instead of a silent re-gold.
 [quantities.integrator]
 h5path = "Input/gpec_toml_raw"
 type = "token"
@@ -28,14 +25,7 @@ label = "integrator (Δ′ provenance)"
 noise_threshold = 0
 order = 1
 
-[quantities.gal_flag]
-h5path = "Input/gpec_toml_raw"
-type = "token"
-extract = "toml_key:ForceFreeStates.gal_flag"
-label = "gal_flag (Δ′ provenance)"
-noise_threshold = 0
-order = 2
-
+# Surface locations / count (sanity)
 [quantities.pe_rational_q]
 h5path = "PerturbedEquilibrium/SingularCoupling/rational_q"
 type = "real_vector"