Changeset 12918

Timestamp:

08/06/12 17:35:49 (13 years ago)

Author:

cborstad

Message:

CHG: return all strain rates in units of 1/a for consistency

File:

• issm/branches/trunk-jpl-damage/src/m/model/mechanicalproperties.m (modified) (2 diffs)

issm/branches/trunk-jpl-damage/src/m/model/mechanicalproperties.m ¶

@@ -3 +3 @@
 %
 %   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 +99 @@
 
 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);