f_psi_sign has been added for matfile in torpex.py for signed Bp

boutproject · Sep 18, 2024 · 7653245 · 7653245
1 parent f7e7882
commit 7653245
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/hypnotoad/cases/torpex.py b/hypnotoad/cases/torpex.py
@@ -297,6 +297,9 @@ def __init__(self, equilibOptions, meshOptions):
                     f"{eqfile['Z'][0, 0][ZindMid, RindMid]}"
                 )
             self.Bt_axis = Bt[ZindMid, RindMid]
+
+            self.f_psi_sign = numpy.sign(0.0 - psi[ZindMid, RindMid])  # H.seto
+
         else:
             raise ValueError("Failed to initialise psi function from inputs")