Index: /issm/trunk/src/c/parallel/OutputControl.cpp
===================================================================
--- /issm/trunk/src/c/parallel/OutputControl.cpp	(revision 871)
+++ /issm/trunk/src/c/parallel/OutputControl.cpp	(revision 872)
@@ -22,7 +22,9 @@
 	Result*  result=NULL;
 	int numberofnodes;
+	int gsize;
 
 	/*recover numberofnodes: */
 	fem->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
+	gsize=fem->nodesets->GetGSize();
 
 	/*Initialize results*/
@@ -33,5 +35,5 @@
 	results->AddObject(result);
 
-	result=new Result(results->Size()+1,0,1,"param_g",param_g,numberofnodes);
+	result=new Result(results->Size()+1,0,1,"param_g",param_g,gsize);
 	results->AddObject(result);
 
Index: /issm/trunk/src/c/parallel/control.cpp
===================================================================
--- /issm/trunk/src/c/parallel/control.cpp	(revision 871)
+++ /issm/trunk/src/c/parallel/control.cpp	(revision 872)
@@ -151,5 +151,5 @@
 		ControlConstrainx( param_g,gsize,mincontrolconstraint,maxcontrolconstraint,control_type);
 		_printf_("%s\n","      done.");
-	
+
 		//some temporary saving 
 		if (((n+1)%5)==0){
@@ -157,4 +157,5 @@
 			inputs->Add(control_type,param_g,2,numberofnodes);
 			inputs->Add("fit",fit[n]);
+			UpdateFromInputsx(femmodel.elements,femmodel.nodes,femmodel.loads, femmodel.materials,inputs);
 			diagnostic_core_nonlinear(&u_g,NULL,NULL,&femmodel,inputs,analysis_type,sub_analysis_type);
 			OutputControl(u_g,param_g,J,nsteps,&femmodel,outputfilename);