Index: /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp
===================================================================
--- /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp	(revision 8807)
+++ /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp	(revision 8808)
@@ -27,4 +27,7 @@
 	if(pf && global_m*global_n){
 
+		/*Some checks in debugging mode*/
+		_assert_(y_s);
+
 		/*pf = pf - Kfs * y_s;*/
 		MatGetLocalSize(Kfs,&Kfsm,&Kfsn);
Index: /issm/trunk/src/c/objects/Node.cpp
===================================================================
--- /issm/trunk/src/c/objects/Node.cpp	(revision 8807)
+++ /issm/trunk/src/c/objects/Node.cpp	(revision 8808)
@@ -119,5 +119,5 @@
 		if (iomodel->dim==3){
 			/*On a 3d mesh, we may have collapsed elements, hence dead nodes. Freeze them out: */
-			if (!iomodel->nodeonbed) _error_("iomodel->nodeonbed is NULL");
+			_assert_(iomodel->nodeonbed);
 			if (!iomodel->nodeonbed[io_index]){
 				for(k=1;k<=gsize;k++){
Index: /issm/trunk/src/c/solvers/solver_linear.cpp
===================================================================
--- /issm/trunk/src/c/solvers/solver_linear.cpp	(revision 8807)
+++ /issm/trunk/src/c/solvers/solver_linear.cpp	(revision 8808)
@@ -17,11 +17,11 @@
 	Vec  df  = NULL;
 	Vec  ys  = NULL;
-	int  analysis_type;
+	int  configuration_type;
 
 	/*Recover parameters: */
-	femmodel->parameters->FindParam(&analysis_type,AnalysisTypeEnum);
+	femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
 
 	SystemMatricesx(&Kff, &Kfs, &pf, &df, NULL,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
-	CreateNodalConstraintsx(&ys,femmodel->nodes,analysis_type);
+	CreateNodalConstraintsx(&ys,femmodel->nodes,configuration_type);
 	Reduceloadx(pf, Kfs, ys); MatFree(&Kfs);
 	Solverx(&uf, Kff, pf, NULL, df, femmodel->parameters); MatFree(&Kff); VecFree(&pf); VecFree(&df);