Index: /issm/trunk-jpl/src/c/solutions/diagnostic_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/diagnostic_core.cpp	(revision 14413)
+++ /issm/trunk-jpl/src/c/solutions/diagnostic_core.cpp	(revision 14414)
@@ -22,8 +22,8 @@
 	bool  isl1l2            = false;
 	bool  isstokes          = false;
-	bool  isnewton          = false;
 	bool  conserve_loads    = true;
 	bool  modify_loads      = true;
 	bool  save_results;
+	int   newton;
 	int   solution_type;
 	int   numoutputs        = 0;
@@ -36,5 +36,5 @@
 	femmodel->parameters->FindParam(&isl1l2,FlowequationIsl1l2Enum);
 	femmodel->parameters->FindParam(&isstokes,FlowequationIsstokesEnum);
-	femmodel->parameters->FindParam(&isnewton,DiagnosticIsnewtonEnum);
+	femmodel->parameters->FindParam(&newton,DiagnosticIsnewtonEnum);
 	femmodel->parameters->FindParam(&dakota_analysis,QmuIsdakotaEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
@@ -76,6 +76,6 @@
 		if(VerboseSolution()) _pprintLine_("   computing velocities");
 		femmodel->SetCurrentConfiguration(DiagnosticHorizAnalysisEnum);
-		if(isnewton)
-		 solver_newton(femmodel); 
+		if(newton>0)
+		 solver_newton(femmodel);
 		else
 		 solver_nonlinear(femmodel,conserve_loads); 
Index: /issm/trunk-jpl/src/c/solvers/solver_newton.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solvers/solver_newton.cpp	(revision 14413)
+++ /issm/trunk-jpl/src/c/solvers/solver_newton.cpp	(revision 14414)
@@ -15,5 +15,5 @@
 	/*intermediary: */
 	bool   converged;
-	int    count;
+	int    count,newton;
 	IssmDouble kmax;
 	Matrix<IssmDouble>* Kff = NULL;
@@ -37,4 +37,5 @@
 	femmodel->parameters->FindParam(&max_nonlinear_iterations,DiagnosticMaxiterEnum);
 	femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
+	femmodel->parameters->FindParam(&newton,DiagnosticIsnewtonEnum);
 	femmodel->UpdateConstraintsx();