🔄 Update f_res_plasma_neo calculation to include valid aspect ratio r…

…ange
ukaea · Feb 17, 2025 · 9bc6d97 · 9bc6d97
1 parent 94650de
commit 9bc6d97
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/process/physics.py b/process/physics.py
@@ -3147,7 +3147,9 @@ def plasma_ohmic_heating(
         # Taken from ITER Physics Design Guidelines: 1989
         # The expression is valid for aspect ratios in the range 2.5 to 4.0
 
-        f_res_plasma_neo = 1.0 if rmajor / rminor < 2.5 else 4.3 - 0.6 * rmajor / rminor
+        f_res_plasma_neo = (
+            1.0 if 2.5 >= rmajor / rminor <= 4.0 else 4.3 - 0.6 * rmajor / rminor
+        )
 
         res_plasma = res_plasma * f_res_plasma_neo