diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90
index 4b41f9ad2..bcb87a25d 100644
--- a/src/plastic_phenopowerlaw.f90
+++ b/src/plastic_phenopowerlaw.f90
@@ -603,7 +603,7 @@ subroutine plastic_phenopowerlaw_init(fileUnit)
      state   (instance)%s_twin=>plasticState(phase)%state   (startIndex:endIndex,:)
      dotState(instance)%s_twin=>plasticState(phase)%dotState(startIndex:endIndex,:)
      plasticState(phase)%state0(startIndex:endIndex,:) = &
-       spread(param(instance)%tau0_twin(1:totalNtwin(instance)),2,NipcMyPhase)
+       spread(math_expand(param(instance)%tau0_twin, param(instance)%Ntwin), 2, NipcMyPhase)
      plasticState(phase)%aTolState(startIndex:endIndex) = param(instance)%aTolResistance
 
      startIndex = endIndex + 1_pInt