diff --git a/trunk/constitutive_dislo.f90 b/trunk/constitutive_dislo.f90
index 8c5cfd527..8ad80812c 100644
--- a/trunk/constitutive_dislo.f90
+++ b/trunk/constitutive_dislo.f90
@@ -720,7 +720,7 @@ do i=1,material_maxN
    do k=1,constitutive_maxNstatevars
 !* Hardening type *
       do l=1,constitutive_maxNstatevars
-	     if (crystal_SlipIntType(j,k,l)==l) then
+	     if (crystal_SlipIntType(j,k,i)==l) then
 		    K_inter=material_SlipIntCoeff(l,i)
 		 else
 		    K_inter=0.0_pReal
diff --git a/trunk/constitutive_pheno.f90 b/trunk/constitutive_pheno.f90
index 103b40512..f41997d11 100644
--- a/trunk/constitutive_pheno.f90
+++ b/trunk/constitutive_pheno.f90
@@ -691,11 +691,13 @@ do i=1,material_maxN
    do k=1,constitutive_maxNstatevars
 !* Hardening type *
       do l=1,constitutive_maxNstatevars
-	     if (crystal_SlipIntType(j,k,l)==l) then
-		    K_inter=material_SlipIntCoeff(l,i)
-		 else
+	     if (crystal_SlipIntType(j,k,i)==0) then
 		    K_inter=0.0_pReal
-		 endif
+		 elseif (crystal_SlipIntType(j,k,i)==1) then
+		    K_inter=material_SlipIntCoeff(1,i)
+		 else
+		    K_inter=material_SlipIntCoeff(2,i)
+	     endif
       enddo
       constitutive_HardeningMatrix(j,k,i)=K_inter
    enddo
diff --git a/trunk/mattex.mpie b/trunk/mattex.mpie
index 49c4dc2a6..96d6d39da 100644
--- a/trunk/mattex.mpie
+++ b/trunk/mattex.mpie
@@ -5,15 +5,19 @@ Nslip				12
 C11				106.75e3
 C12				60.41e3
 C44				28.34e3
-# Parameters for phenomenological modeling
+
+## Parameters for phenomenological modeling
+# Unit in [MPa]
 s0_slip				31.0
 gdot0_slip			0.001
 n_slip				20
 h0				75
 s_sat				63
 w0				2.25
+# Self and latent hardening coefficients
 hardening_coefficients        	1.0 1.4
-# Parameters for dislocation-based modeling
+
+## Parameters for dislocation-based modeling
 # Initial dislocation density [m]²
 rho0				              1.5e11 
 # Burgers vector [m]