Changeset 631

Timestamp:

05/28/09 15:05:08 (16 years ago)

Author:

Mathieu Morlighem

Message:

minor

Location:

issm/trunk/src/m/classes/public

Files:

• ismodelselfconsistent.m (modified) (1 diff)
• mechanicalproperties.m (modified) (3 diffs)

issm/trunk/src/m/classes/public/ismodelselfconsistent.m

@@ -265 +265 @@
 if strcmp(md.analysis_type,'parameters')
 
+        %PACKAGE
+        if ~strcmpi(packagem,'ice'),
+                disp('parameter solution only supported by package ice yet');
+                bool=0;return;
+        end
+
         %OUTPUT
         if ~iscell(md.parameteroutput)

issm/trunk/src/m/classes/public/mechanicalproperties.m

@@ -68 +68 @@
 nu(location)=B_bar(location);
 clear B_bar location second_inv power
-mean(nu)
 
 %compute stress
@@ -104 +103 @@
 stress.principalaxis1=directionsstress(:,3:4);
 stress.effectivevalue=1/sqrt(2)*sqrt(stress.xx.^2+stress.yy.^2+2*stress.xy.^2);
+md.stress=stress;
+
 strainrate=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]);
 strainrate.xx=ux;
@@ -113 +114 @@
 strainrate.principalaxis1=directionsstrain(:,3:4);
 strainrate.effectivevalue=1/sqrt(2)*sqrt(strainrate.xx.^2+strainrate.yy.^2+2*strainrate.xy.^2);
-md.stress=stress;
 md.strainrate=strainrate;