diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index c60fb150a339d9c4a4db8353245d24975fa1c348..277a12626e3abac6aa5b3d56110b72620713e231 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -6192,7 +6192,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
             for(int k=0; k<3; k++)
               for(int l=0; l<3; l++)
                 for(int p=0; p<3; p++){
-                  dUdF(m,n,k,l) += I4_first_type(m,p,k,l)*R_def(n,p);
+                  dUdF(m,n,k,l) += R_def(p,m)*I4_first_type(p,n,k,l); // R_def_T is used here
                   for(int q=0; q<3; q++)
                     dUrot_dU(m,n,k,l) += R(m,p)*I4_symm(p,q,k,l)*RT(q,n);
                 }
@@ -6223,7 +6223,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
         // dPdF
         static STensor43 dSdF, dSrot_dF, L_rot;
         static STensor33 dqdF_rot;
-        static STensor3 dqdgradT_rot, dwdf_rot, dmechSourcedf_rot;
+        static STensor3 dwdf_rot, dmechSourcedf_rot;
         
         STensorOperation::zero(L_rot);
         STensorOperation::zero(dqdF_rot);
@@ -6271,7 +6271,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
         STensorOperation::multSTensor3STensor43(dmechSourcedf,dFrot_dF,dmechSourcedf_rot);
           
         // Reassign
-        L = L_rot; dqdF = dqdF_rot; dqdgradT = dqdgradT_rot; dwdf = dwdf_rot; dmechSourcedf = dmechSourcedf_rot; 
+        L = L_rot; dqdF = dqdF_rot; dwdf = dwdf_rot; dmechSourcedf = dmechSourcedf_rot; 
       }
     
       STensorOperation::multSTensor3(F_unrot,S,P);