Location 1 PF coils don't work if there is no CS supersedes #3364

process/pfcoil.py

@@ -166,7 +166,7 @@ def pfcoil(self):
         for j in range(pfv.ngrp):
             if pfv.ipfloc[j] == 1:
                 # PF coil is stacked on top of the Central Solenoid
-                for k in pfv.ncls[j]:
+                for k in range(pfv.ncls[j]):
                     pf.rcls[j, k] = pfv.rohc + pfv.rpf1
 
                     # Z coordinate of coil enforced so as not
@@ -327,6 +327,9 @@ def pfcoil(self):
                 nocoil = 0
                 for i in range(pfv.ngrp):
                     if pfv.ipfloc[i] == 1:
+                        # Do not allow if no central solenoid
+                        if bv.iohcl == 0:
+                            eh.report_error(288)
                         # PF coil is stacked on top of the Central Solenoid
                         # This coil is to balance Central Solenoid flux and should not be involved
                         # in equilibrium calculation -- RK 07/12

process/utilities/errorlist.json

@@ -8,7 +8,7 @@
   "comment2": [
     "Increment n_errortypes if an error is added to this list"
   ],
-  "n_errortypes": 287,
+  "n_errortypes": 288,
   "errors": [
     {
       "no": 1,
@@ -1444,6 +1444,11 @@
       "no": 287,
       "level": 2,
       "message": "[pfcoil][cntrpost]: fjohc0 shouldn't be above 0.7 for engineering reliability."
+    },
+    {
+      "no": 288,
+      "level": 3,
+      "message": "PFCOIL: ipfloc(i) should not be 1 if iohcl=0"
     }
   ]
 }