diff --git a/code/constitutive_phenopowerlaw.f90 b/code/constitutive_phenopowerlaw.f90
index 0e0934c8f..3599a479a 100644
--- a/code/constitutive_phenopowerlaw.f90
+++ b/code/constitutive_phenopowerlaw.f90
@@ -605,7 +605,7 @@ pure function constitutive_phenopowerlaw_homogenizedC(state,ipc,ip,el)
    state                                                                                            ! state variables
  
  matID = phase_plasticityInstance(material_phase(ipc,ip,el))
- constitutive_phenopowerlaw_homogenizedC = constitutive_phenopowerlaw_Cslip_66(6,6,matID)
+ constitutive_phenopowerlaw_homogenizedC = constitutive_phenopowerlaw_Cslip_66(1:6,1:6,matID)
 
 end function constitutive_phenopowerlaw_homogenizedC