Index: /issm/branches/trunk-jpl-damage/src/m/model/mechanicalproperties.m
===================================================================
--- /issm/branches/trunk-jpl-damage/src/m/model/mechanicalproperties.m	(revision 12917)
+++ /issm/branches/trunk-jpl-damage/src/m/model/mechanicalproperties.m	(revision 12918)
@@ -3,5 +3,5 @@
 %
 %   this routine computes the components of the stress tensor
-%   strain rate tensor and their respective principal directions.
+%   strain rate tensor (1/a) and their respective principal directions.
 %   the results are in the model md: md.results
 %
@@ -99,7 +99,7 @@
 
 strainrate=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]);
-strainrate.xx=ux;
-strainrate.yy=vy;
-strainrate.xy=uyvx;
+strainrate.xx=ux*(365.25*24*3600); %strain rate in 1/a instead of 1/s
+strainrate.yy=vy*(365.25*24*3600); %strain rate in 1/a instead of 1/s
+strainrate.xy=uyvx*(365.25*24*3600); %strain rate in 1/a instead of 1/s
 strainrate.principalvalue2=valuesstrain(:,1)*(365.25*24*3600); %strain rate in 1/a instead of 1/s
 strainrate.principalaxis2=directionsstrain(:,1:2);