diff --git a/trunk/constitutive.f90 b/trunk/constitutive.f90
index fbd862003..273c3e691 100644
--- a/trunk/constitutive.f90
+++ b/trunk/constitutive.f90
@@ -280,7 +280,7 @@ do l=1,3
       constitutive_Sslip_v(3,k,l)=constitutive_Sslip(3,3,k,l)
 	  !* be compatible with Mandel notation of Tstar
       constitutive_Sslip_v(4,k,l)=(constitutive_Sslip(1,2,k,l)+constitutive_Sslip(2,1,k,l))/dsqrt(2.0_pReal)
-      constitutive_Sslip_v(5,k,l)=(constitutive_Sslip(2,3,k,l)+constitutive_Sslip(3,3,k,l))/dsqrt(2.0_pReal)
+      constitutive_Sslip_v(5,k,l)=(constitutive_Sslip(2,3,k,l)+constitutive_Sslip(3,2,k,l))/dsqrt(2.0_pReal)
       constitutive_Sslip_v(6,k,l)=(constitutive_Sslip(1,3,k,l)+constitutive_Sslip(3,1,k,l))/dsqrt(2.0_pReal)
    enddo
 enddo
diff --git a/trunk/mpie_cpfem_marc.f90 b/trunk/mpie_cpfem_marc.f90
index 9cbd39b58..308eac601 100644
--- a/trunk/mpie_cpfem_marc.f90
+++ b/trunk/mpie_cpfem_marc.f90
@@ -178,7 +178,7 @@
  cp_en = mesh_FEasCP('elem', n(1))
  s(1:ngens)=invnrmMandel(1:ngens)*CPFEM_stress_all(1:ngens, nn, cp_en)
  d(1:ngens,1:ngens)=CPFEM_jaco_old(1:ngens,1:ngens, nn, cp_en)
- forall(i=1:ngens) d(i,1:ngens)=d(i,1:ngens)*invnrmMandel(1:ngens)
+ forall(i=1:ngens) d(1:ngens,i)=d(1:ngens,1)*invnrmMandel(1:ngens)
  return
  
  END SUBROUTINE