diff --git a/code/plastic_phenopowerlaw.f90 b/code/plastic_phenopowerlaw.f90
index 2f9051c71..28a8946c8 100644
--- a/code/plastic_phenopowerlaw.f90
+++ b/code/plastic_phenopowerlaw.f90
@@ -1013,7 +1013,7 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el)
      enddo nonSchmidSystems
      gdot_slip(j) = plastic_phenopowerlaw_gdot0_slip(instance)*0.5_pReal* &
                   ((abs(tau_slip_pos)/(plasticState(ph)%state(j,of)))**plastic_phenopowerlaw_n_slip(instance) &
-                  *sign(1.0_pReal,tau_slip_pos)
+                  *sign(1.0_pReal,tau_slip_pos) &
                   +(abs(tau_slip_neg)/(plasticState(ph)%state(j,of)))**plastic_phenopowerlaw_n_slip(instance) &
                   *sign(1.0_pReal,tau_slip_neg))
    enddo slipSystems1