Index: /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 20965)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 20966)
@@ -644,7 +644,7 @@
 	}
 	else{
-		/* eps_eff^2 = 1/2 (exx^2 + 2*exy^2 )*/
+		/* eps_eff^2 = 1/2 (2*exx^2 + 2*exy^2 ) (since eps_zz = - eps_xx)*/
 		element->StrainRateHO2dvertical(&epsilon2d[0],xyz_list,gauss,vx_input,vy_input);
-		eps_eff = 1./sqrt(2.)*sqrt(epsilon2d[0]*epsilon2d[0] + 2.*epsilon2d[1]*epsilon2d[1]);
+		eps_eff = 1./sqrt(2.)*sqrt(2*epsilon2d[0]*epsilon2d[0] + 2*epsilon2d[1]*epsilon2d[1]);
 	}
 
@@ -732,7 +732,7 @@
 	}
 	else{
-		/* eps_eff^2 = 1/2 exx^2*/
+		/* eps_eff^2 = exx^2*/
 		element->StrainRateSSA1d(&epsilon1d,xyz_list,gauss,vx_input);
-		eps_eff = sqrt(epsilon1d*epsilon1d/2.);
+		eps_eff = fabs(epsilon1d);
 	}