Changeset 14296 for issm/trunk-jpl/src/m/mech/mechanicalproperties.m

Timestamp:

01/31/13 19:43:33 (12 years ago)

Author:

cborstad

Message:

NEW: function to calculate error in strain rate components from velocity error

File:

• issm/trunk-jpl/src/m/mech/mechanicalproperties.m (modified) (1 diff)

issm/trunk-jpl/src/m/mech/mechanicalproperties.m ¶

@@ -119 +119 @@
 
 strainrate=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]);
-strainrate.xx=ux;
-strainrate.yy=vy;
-strainrate.xy=uyvx;
-strainrate.principalvalue1=valuesstrain(:,1)*(365.25*24*3600); %strain rate in 1/a instead of 1/s
+strainrate.xx=ux*md.constants.yts; %strain rate in 1/a instead of 1/s
+strainrate.yy=vy*md.constants.yts; 
+strainrate.xy=uyvx*md.constants.yts; 
+strainrate.principalvalue1=valuesstrain(:,1)*md.constants.yts; 
 strainrate.principalaxis1=directionsstrain(:,1:2);
-strainrate.principalvalue2=valuesstrain(:,2)*(365.25*24*3600); %strain rate in 1/a instead of 1/s
+strainrate.principalvalue2=valuesstrain(:,2)*md.constants.yts; 
 strainrate.principalaxis2=directionsstrain(:,3:4);
 strainrate.effectivevalue=1/sqrt(2)*sqrt(strainrate.xx.^2+strainrate.yy.^2+2*strainrate.xy.^2);