From ab5ff594b59d7dfab51b86e328fa47ba2df9de8f Mon Sep 17 00:00:00 2001
From: Jonathan Maddock <78556175+jonmaddock@users.noreply.github.com>
Date: Fri, 7 Mar 2025 14:17:25 +0000
Subject: [PATCH] Warn about fcohbop > 1 rather than preventing it

---
 process/pfcoil.py | 12 +++++-------
 1 file changed, 5 insertions(+), 7 deletions(-)

diff --git a/process/pfcoil.py b/process/pfcoil.py
index 4a91b843a9..8e650e978f 100644
--- a/process/pfcoil.py
+++ b/process/pfcoil.py
@@ -488,13 +488,11 @@ def pfcoil(self):
             pfv.f_j_cs_start_end_flat_top = (
                 (-ioheof * pfv.f_j_cs_start_pulse_end_flat_top) + dics
             ) / ioheof
-            pfv.f_j_cs_start_end_flat_top = min(
-                pfv.f_j_cs_start_end_flat_top, 1.0e0
-            )  # constrains abs(f_j_cs_start_end_flat_top) <= 1.0;
-            pfv.f_j_cs_start_end_flat_top = max(
-                pfv.f_j_cs_start_end_flat_top, -1.0e0
-            )  # probably un-necessary
-
+            if np.abs(pfv.f_j_cs_start_end_flat_top) > 1.0:
+                logger.warning(
+                    "Ratio of central solenoid overall current density at "
+                    "beginning of flat-top / end of flat-top > 1 (|f_j_cs_start_end_flat_top| > 1)"
+                )
         else:
             dics = 0.0e0
             pfv.f_j_cs_start_end_flat_top = 1.0e0