Use properly computed Ip_non_inductive in ActorQED

src/actors/current/qed_actor.jl

@@ -60,7 +60,6 @@ function _step(actor::ActorQED)
     cp1d = dd.core_profiles.profiles_1d[]
 
     # non_inductive contribution
-    Ip_non_inductive = trapz(cp1d.grid.area, cp1d.j_non_inductive)
     B0 = eqt.global_quantities.vacuum_toroidal_field.b0
     JBni = QED.FE(cp1d.grid.rho_tor_norm, cp1d.j_non_inductive .* B0)
 
@@ -93,6 +92,7 @@ function _step(actor::ActorQED)
         for time0 in range(t0 + δt / 2.0, t1 + δt / 2.0, No + 1)[1:end-1]
             if par.solve_for == :ip
                 Ip = IMAS.get_from(dd, Val{:ip}, par.ip_from; time0)
+                Ip_non_inductive = dd.core_profiles.global_quantities.current_non_inductive[time0]
                 Vedge = nothing
                 if par.allow_floating_plasma_current && abs(Ip) < abs(Ip_non_inductive)
                     Ip = nothing
@@ -109,6 +109,7 @@ function _step(actor::ActorQED)
         # steady state solution
         if par.solve_for == :ip
             Ip = IMAS.get_from(dd, Val{:ip}, par.ip_from)
+            Ip_non_inductive = dd.core_profiles.global_quantities.current_non_inductive[]
             Vedge = nothing
             if par.allow_floating_plasma_current && abs(Ip) < abs(Ip_non_inductive)
                 Ip = nothing