Index: /issm/trunk/src/m/solutions/jpl/balancedthickness.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedthickness.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/balancedthickness.m	(revision 4120)
@@ -21,5 +21,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		balancedthickness_core(femmodel);
+		femmodel=balancedthickness_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/jpl/balancedthickness2.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedthickness2.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/balancedthickness2.m	(revision 4120)
@@ -18,5 +18,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		balancedthickness2_core(femmodel);
+		femmodel=balancedthickness2_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/jpl/balancedvelocities.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedvelocities.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/balancedvelocities.m	(revision 4120)
@@ -21,5 +21,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		balancedvelocities_core(femmodel);
+		femmodel=balancedvelocities_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/jpl/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/prognostic.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/prognostic.m	(revision 4120)
@@ -22,5 +22,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		prognostic_core(femmodel);
+		femmodel=prognostic_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/jpl/prognostic2.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/prognostic2.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/prognostic2.m	(revision 4120)
@@ -21,5 +21,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		prognostic2_core(femmodel);
+		femmodel=prognostic2_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/jpl/steadystate.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/steadystate.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/steadystate.m	(revision 4120)
@@ -25,10 +25,10 @@
 			
 			displaystring(verbose,'%s',['call computational core']);
-			steadystate_core(femmodel);
+			femmodel=steadystate_core(femmodel);
 
 		else,
 			
 			displaystring(verbose,'%s',['call computational core']);
-			control_core(femmodel);
+			femmodel=control_core(femmodel);
 
 		end
Index: /issm/trunk/src/m/solutions/jpl/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/thermal.m	(revision 4119)
+++ /issm/trunk/src/m/solutions/jpl/thermal.m	(revision 4120)
@@ -21,5 +21,5 @@
 	if ~qmu_analysis,
 		displaystring(verbose,'%s',['call computational core']);
-		thermal_core(femmodel);
+		femmodel=thermal_core(femmodel);
 	else
 		%launch dakota driver for diagnostic core solution