Prohibit hyper_cleaning in 2D simulations (p == 0)

The psi field used by hyperbolic cleaning is initialized with a
3D loop that accesses z_cc, which is out-of-bounds when p == 0.
The 1D prohibition already existed; extend it to 2D.

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

Author: sbryngelson

toolchain/mfc/case_validator.py

@@ -1029,6 +1029,7 @@ def check_mhd_simulation(self):
         relativity = self.get('relativity', 'F') == 'T'
         hyper_cleaning = self.get('hyper_cleaning', 'F') == 'T'
         n = self.get('n', 0)
+        p = self.get('p', 0)
 
         self.prohibit(mhd and riemann_solver is not None and riemann_solver not in [1, 4],
                      "MHD simulations require riemann_solver = 1 (HLL) or riemann_solver = 4 (HLLD)")
@@ -1040,6 +1041,8 @@ def check_mhd_simulation(self):
                      "Hyperbolic cleaning requires mhd to be enabled")
         self.prohibit(hyper_cleaning and n is not None and n == 0,
                      "Hyperbolic cleaning is not supported for 1D simulations")
+        self.prohibit(hyper_cleaning and p is not None and p == 0,
+                     "Hyperbolic cleaning is not supported for 2D simulations")
 
 
     def check_igr_simulation(self):  # pylint: disable=too-many-locals