Index: /issm/trunk/src/c/solutions/control_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/control_core.cpp	(revision 4522)
+++ /issm/trunk/src/c/solutions/control_core.cpp	(revision 4523)
@@ -105,5 +105,5 @@
 		/*Temporary saving every 5 control steps: */
 		if (((n+1)%5)==0){
-			_printf_("%s","      saving temporary results...");
+			_printf_("%s\n","      saving temporary results...");
 			controlrestart(femmodel,J);
 		}
@@ -112,5 +112,5 @@
 
 	convergence_point:
-	_printf_("%s","      preparing final velocity solution");
+	_printf_("%s\n","      preparing final velocity solution");
 	if (control_steady) steadystate_core(femmodel);
 	else diagnostic_core(femmodel);
@@ -121,7 +121,7 @@
 	InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VyEnum); 
 	InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VelEnum); 
-	InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,GradientEnum); 
-	InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointxEnum); 
-	InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointyEnum); 
+	//InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,GradientEnum); 
+	//InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointxEnum); 
+	//InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,AdjointyEnum); 
 	if(dim==3)InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VzEnum); 
 	femmodel->results->AddObject(new DoubleVecExternalResult(femmodel->results->Size()+1,JEnum,J,nsteps,1,0));
Index: /issm/trunk/src/m/utils/Nightly/testsgetfields.m
===================================================================
--- /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4522)
+++ /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4523)
@@ -111,8 +111,12 @@
 		vel=zeros(md.numberofgrids,1);
 		vel(md.results.SteadystateSolution.Vel.index)=md.results.SteadystateSolution.Vel.value;
-		j=zeros(md.numberofgrids,1);
-		j(md.results.SteadystateSolution.J.index)=md.results.SteadystateSolution.J.value;
-		parameter=zeros(md.numberofgrids,1);
-		parameter(md.results.SteadystateSolution.parameter.index)=md.results.SteadystateSolution.parameter.value;
+		j=md.results.DiagnosticSolution.J;
+		if strcmpi(md.control_type,'rheology_B'),
+			parameter(md.results.DiagnosticSolution.RheologyB.index)=md.results.DiagnosticSolution.RheologyB.value;
+		elseif strcmpi(md.control_type,'drag_coefficient'),
+			parameter(md.results.DiagnosticSolution.DragCoefficient.index)=md.results.DiagnosticSolution.DragCoefficient.value;
+		else
+			error('not implemented yet');
+		end
 		field_values={vel,j,parameter};
 	else