Fixed negativity in sqrt in costs.py by np.clip.

Author: OceanNuclear

process/costs.py

@@ -3026,7 +3026,9 @@ def coelc(self):
         #  Annual cost of operation and maintenance
 
         annoam = cost_variables.ucoam[cost_variables.lsa - 1] * np.sqrt(
-            heat_transport_variables.p_plant_electric_net_mw / 1200.0e0
+            np.clip(
+                heat_transport_variables.p_plant_electric_net_mw, 0, np.inf
+            ) / 1200.0e0
         )
 
         #  Additional cost due to pulsed reactor thermal storage
@@ -3097,7 +3099,9 @@ def coelc(self):
         #  Annual cost of waste disposal
 
         annwst = cost_variables.ucwst[cost_variables.lsa - 1] * np.sqrt(
-            heat_transport_variables.p_plant_electric_net_mw / 1200.0e0
+            np.clip(
+                heat_transport_variables.p_plant_electric_net_mw, 0, np.inf
+            ) / 1200.0e0
         )
 
         #  Cost of electricity due to waste disposal