Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 630)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 631)
@@ -265,4 +265,10 @@
 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)
Index: /issm/trunk/src/m/classes/public/mechanicalproperties.m
===================================================================
--- /issm/trunk/src/m/classes/public/mechanicalproperties.m	(revision 630)
+++ /issm/trunk/src/m/classes/public/mechanicalproperties.m	(revision 631)
@@ -68,5 +68,4 @@
 nu(location)=B_bar(location);
 clear B_bar location second_inv power
-mean(nu)
 
 %compute stress
@@ -104,4 +103,6 @@
 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,4 +114,3 @@
 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;