Index: /issm/trunk-jpl/src/m/mech/mechanicalproperties.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/mechanicalproperties.m	(revision 22753)
+++ /issm/trunk-jpl/src/m/mech/mechanicalproperties.m	(revision 22754)
@@ -111,16 +111,4 @@
 %plug onto the model
 %NB: Matlab sorts the eigen value in increasing order, we want the reverse
-%NOTE: all stresses below are deviatoric
-stress=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]);
-stress.xx=tau_xx;
-stress.yy=tau_yy;
-stress.xy=tau_xy;
-stress.principalvalue1=valuesstress(:,1);
-stress.principalaxis1=directionsstress(:,1:2);
-stress.principalvalue2=valuesstress(:,2);
-stress.principalaxis2=directionsstress(:,3:4);
-stress.effectivevalue=1/sqrt(2)*sqrt(stress.xx.^2+stress.yy.^2+2*stress.xy.^2); % effective shear stress (sqrt(J2))
-md.results.stress=stress;
-
 strainrate=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]);
 strainrate.xx=ux*md.constants.yts; %strain rate in 1/a instead of 1/s
Index: /issm/trunk-jpl/src/m/mech/mechanicalproperties.py
===================================================================
--- /issm/trunk-jpl/src/m/mech/mechanicalproperties.py	(revision 22753)
+++ /issm/trunk-jpl/src/m/mech/mechanicalproperties.py	(revision 22754)
@@ -125,14 +125,4 @@
 	##plug onto the model
 	##NB: Matlab sorts the eigen value in increasing order, we want the reverse
-	stress=results()
-	stress.xx=tau_xx
-	stress.yy=tau_yy
-	stress.xy=tau_xy
-	stress.principalvalue1=valuesstress[:,0]
-	stress.principalaxis1=directionsstress[:,0:2]
-	stress.principalvalue2=valuesstress[:,1]
-	stress.principalaxis2=directionsstress[:,2:4]
-	stress.effectivevalue=1./np.sqrt(2.)*np.sqrt(stress.xx**2+stress.yy**2+2.*stress.xy**2)
-	md.results.stress=stress
 	
 	strainrate=results()