Index: /issm/trunk/src/m/solutions/control_core.m
===================================================================
--- /issm/trunk/src/m/solutions/control_core.m	(revision 4554)
+++ /issm/trunk/src/m/solutions/control_core.m	(revision 4555)
@@ -90,7 +90,7 @@
 %	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parametersn,AdjointxEnum);
 %	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parametersn,AdjointyEnum);
-	if (dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
-	femmodel.results.JEnum=J;
-	femmodel.results.ControlTypeEnum=EnumAsString(control_type);
+if (dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);end
+	femmodel.results=AddExternalResult(femmodel.results,JEnum,J);
+	femmodel.results=AddExternalResult(femmodel.results,ControlTypeEnum,control_type);
 
 end %end function
Index: /issm/trunk/src/m/solutions/stokescontrolinit.m
===================================================================
--- /issm/trunk/src/m/solutions/stokescontrolinit.m	(revision 4554)
+++ /issm/trunk/src/m/solutions/stokescontrolinit.m	(revision 4555)
@@ -19,8 +19,8 @@
 	%Run a complete diagnostic to update the Stokes spcs:
 	femmodel=SetCurrentConfiguration(femmodel,DiagnosticHorizAnalysisEnum);
-	femmodel=solver_diagnostic_nonlinear(NULL,NULL,NULL,femmodel,conserve_loads);
+	femmodel=solver_diagnostic_nonlinear([],[],[],femmodel,conserve_loads);
 
 	%vertical velocity
-	femmodel=solver_linear(NULL,femmodel);
+	femmodel=solver_linear([],femmodel);
 
 	%recondition pressure computed previously:
@@ -33,3 +33,3 @@
 	displaystring(verbose,'\n%s',['      computing stokes velocity and pressure...']);
 	femmodel=SetCurrentConfiguration(femmodel,DiagnosticStokesAnalysisEnum);
-	femmode=solver_diagnostic_nonlinear(NULL,NULL,NULL,femmodel,conserve_loads);
+	femmode=solver_diagnostic_nonlinear([],[],[],femmodel,conserve_loads);