TokaMaker: Fix bug in time-dependent solves introduced in 6b3439b

hansec · Mar 4, 2024 · 56e47e3 · 56e47e3
1 parent 6b3439b
commit 56e47e3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/physics/grad_shaf_td.F90 b/src/physics/grad_shaf_td.F90
@@ -370,11 +370,12 @@ subroutine step_gs_td(time,dt,nl_its,lin_its,nretry)
     CALL adv_pre%update(.TRUE.)
 END IF
 ! Update coil currents (end of time step)
+tMaker_td_obj%curr_reg=0.d0
 DO i=1,tMaker_td_obj%gs_eq%ncoils
     DO k=1,tMaker_td_obj%gs_eq%ncoil_regs
         j=tMaker_td_obj%gs_eq%coil_regions(k)%id
         tMaker_td_obj%curr_reg(j)=tMaker_td_obj%curr_reg(j) &
-        + tMaker_td_obj%gs_eq%coil_currs(i)*tMaker_td_obj%gs_eq%coil_nturns(j,i)
+            + tMaker_td_obj%gs_eq%coil_currs(i)*tMaker_td_obj%gs_eq%coil_nturns(j,i)
     END DO
 END DO
 ! Point to profiles in case they changed