Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 19024)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 19025)
@@ -75,4 +75,5 @@
 /*}}}*/
 void LevelsetAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
+	parameters->AddObject(iomodel->CopyConstantObject(LevelsetStabilizationEnum));
 	return;
 }
@@ -87,8 +88,8 @@
 
 	/*parameters: */
+	int  stabilization;
 	bool save_results;
-	int  stabilization = 1;
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	//femmodel->parameters->FindParam(&stabilization,LevelsetStabilizationEnum);
+	femmodel->parameters->FindParam(&stabilization,LevelsetStabilizationEnum);
 
 	/*activate formulation: */
@@ -123,6 +124,5 @@
 
 	/*Intermediaries */
-	int  stabilization=2;
-	int  dim, domaintype, calvinglaw;
+	int  stabilization,dim, domaintype, calvinglaw;
 	bool iscalving;
 	int i, row, col;
@@ -138,4 +138,5 @@
 	basalelement->FindParam(&domaintype,DomainTypeEnum);
 	basalelement->FindParam(&calvinglaw,CalvingLawEnum);
+	basalelement->FindParam(&stabilization,LevelsetStabilizationEnum);
 	switch(domaintype){
 		case Domain2DverticalEnum:   dim = 1; break;
Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 19024)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 19025)
@@ -112,4 +112,5 @@
 				if(VerboseSolution()) _printf0_("   computing calving rate\n");
 				femmodel->StrainRateparallelx();
+				femmodel->DeviatoricStressx();
 				femmodel->CalvingRatePix();
 			}