diff --git a/src/misc.f90 b/src/misc.f90
index 1d546dd30..ac8f37839 100644
--- a/src/misc.f90
+++ b/src/misc.f90
@@ -143,7 +143,7 @@ end function misc_prefixOptions
 !--------------------------------------------------------------------------------------------------
 pure function misc_zeros(N)
 
-  integer, intent(in) :: N                                                                          !< number of ones
+  integer, intent(in) :: N                                                                          !< number of zeros
   real(pREAL), dimension(N) :: misc_zeros
 
 
diff --git a/src/phase_mechanical_plastic_phenopowerlaw.f90 b/src/phase_mechanical_plastic_phenopowerlaw.f90
index 2368e7e07..095f72b74 100644
--- a/src/phase_mechanical_plastic_phenopowerlaw.f90
+++ b/src/phase_mechanical_plastic_phenopowerlaw.f90
@@ -147,7 +147,7 @@ module function plastic_phenopowerlaw_init() result(myPlasticity)
       xi_0_sl            = math_expand(pl%get_as1dReal('xi_0_sl',       requiredSize=size(N_sl)), N_sl)
       prm%xi_inf_sl      = math_expand(pl%get_as1dReal('xi_inf_sl',     requiredSize=size(N_sl)), N_sl)
       prm%c_1            = math_expand(pl%get_as1dReal('c_1',           requiredSize=size(N_sl), &
-                                                                        defaultVal=misc_ones(size(N_sl))), N_sl)
+                                                                        defaultVal=misc_zeros(size(N_sl))), N_sl)
       prm%c_2            = math_expand(pl%get_as1dReal('c_2',           requiredSize=size(N_sl), &
                                                                         defaultVal=misc_ones(size(N_sl))), N_sl)
       prm%f_sat_sl_tw    = math_expand(pl%get_as1dReal('f_sat_sl-tw',   requiredSize=size(N_sl), &