Index: ../trunk-jpl/src/c/cores/balancevelocity_core.cpp
===================================================================
--- ../trunk-jpl/src/c/cores/balancevelocity_core.cpp	(revision 18606)
+++ ../trunk-jpl/src/c/cores/balancevelocity_core.cpp	(revision 18607)
@@ -12,13 +12,14 @@
 void balancevelocity_core(FemModel* femmodel){
 
 	/*parameters: */
-	bool save_results;
+	bool        save_results;
+	IssmPDouble l = 8.;
 
 	/*recover parameters: */
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
 
 	if(VerboseSolution()) _printf0_("computing smooth driving stress:\n");
-	femmodel->parameters->SetParam(IssmPDouble(8.),SmoothThicknessMultiplierEnum);
+	femmodel->parameters->SetParam(l,SmoothThicknessMultiplierEnum);
 	femmodel->SetCurrentConfiguration(SmoothAnalysisEnum);
 	femmodel->parameters->SetParam(DrivingStressXEnum,InputToSmoothEnum);
 	solutionsequence_linear(femmodel);
Index: ../trunk-jpl/src/c/cores/balancethickness2_core.cpp
===================================================================
--- ../trunk-jpl/src/c/cores/balancethickness2_core.cpp	(revision 18606)
+++ ../trunk-jpl/src/c/cores/balancethickness2_core.cpp	(revision 18607)
@@ -12,13 +12,14 @@
 void balancethickness2_core(FemModel* femmodel){
 
 	/*parameters: */
-	bool save_results;
+	bool        save_results;
+	IssmPDouble l = 8.;
 
 	/*recover parameters: */
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
 
 	if(VerboseSolution()) _printf0_("computing smooth surface slopes:\n");
-	femmodel->parameters->SetParam(IssmPDouble(8.),SmoothThicknessMultiplierEnum);
+	femmodel->parameters->SetParam(l,SmoothThicknessMultiplierEnum);
 	femmodel->SetCurrentConfiguration(SmoothAnalysisEnum);
 	femmodel->parameters->SetParam(SurfaceSlopeXEnum,InputToSmoothEnum);
 	solutionsequence_linear(femmodel);
