Index: ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 19024)
+++ ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 19025)
@@ -74,6 +74,7 @@
 }
 /*}}}*/
 void LevelsetAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
+	parameters->AddObject(iomodel->CopyConstantObject(LevelsetStabilizationEnum));
 	return;
 }
 /*}}}*/
@@ -86,10 +87,10 @@
 	#endif
 
 	/*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: */
 	femmodel->SetCurrentConfiguration(LevelsetAnalysisEnum);
@@ -122,8 +123,7 @@
 	Element* basalelement = element->SpawnBasalElement();
 
 	/*Intermediaries */
-	int  stabilization=2;
-	int  dim, domaintype, calvinglaw;
+	int  stabilization,dim, domaintype, calvinglaw;
 	bool iscalving;
 	int i, row, col;
 	IssmDouble kappa;
@@ -137,6 +137,7 @@
 	basalelement->FindParam(&iscalving,TransientIscalvingEnum);
 	basalelement->FindParam(&domaintype,DomainTypeEnum);
 	basalelement->FindParam(&calvinglaw,CalvingLawEnum);
+	basalelement->FindParam(&stabilization,LevelsetStabilizationEnum);
 	switch(domaintype){
 		case Domain2DverticalEnum:   dim = 1; break;
 		case Domain2DhorizontalEnum: dim = 2; break;
Index: ../trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- ../trunk-jpl/src/c/cores/transient_core.cpp	(revision 19024)
+++ ../trunk-jpl/src/c/cores/transient_core.cpp	(revision 19025)
@@ -111,6 +111,7 @@
 			if(iscalving && calvinglaw==CalvingPiEnum){
 				if(VerboseSolution()) _printf0_("   computing calving rate\n");
 				femmodel->StrainRateparallelx();
+				femmodel->DeviatoricStressx();
 				femmodel->CalvingRatePix();
 			}
 			if(VerboseSolution()) _printf0_("   computing movement of ice boundaries\n");
