From 9bc6d978119a29f1f707008fa92b4a85817a85d4 Mon Sep 17 00:00:00 2001
From: mn3981 <christopher.ashe@ukaea.uk>
Date: Mon, 17 Feb 2025 11:31:25 +0000
Subject: [PATCH] =?UTF-8?q?=F0=9F=94=84=20Update=20f=5Fres=5Fplasma=5Fneo?=
 =?UTF-8?q?=20calculation=20to=20include=20valid=20aspect=20ratio=20range?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 process/physics.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/process/physics.py b/process/physics.py
index 5c05b2710..9e6c6643a 100644
--- 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