Index: /issm/trunk/src/m/utils/Nightly/testsgetfields.m
===================================================================
--- /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4322)
+++ /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4323)
@@ -25,9 +25,9 @@
 		field_names={'Vel','J','parameter'};
 		vel=zeros(md.numberofgrids,1);
-		vel(md.results.DiagnosticAnalysis.Vel.index)=md.results.DiagnosticAnalysis.Vel.value;
+		vel(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 		j=zeros(md.numberofgrids,1);
-		j(md.results.DiagnosticAnalysis.J.index)=md.results.DiagnosticAnalysis.J.value;
+		j(md.results.DiagnosticSolution.J.index)=md.results.DiagnosticSolution.J.value;
 		parameter=zeros(md.numberofgrids,1);
-		parameter(md.results.DiagnosticAnalysis.parameter.index)=md.results.DiagnosticAnalysis.parameter.value;
+		parameter(md.results.DiagnosticSolution.parameter.index)=md.results.DiagnosticSolution.parameter.value;
 		field_values={vel,j,parameter};
 	elseif qmu,
@@ -37,12 +37,12 @@
 			field_names={'Vy','Vz'};
 			vy=zeros(md.numberofgrids,1);
-			vy(md.results.DiagnosticAnalysis.Vy.index)=md.results.DiagnosticAnalysis.Vy.value;
+			vy(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value;
 			vz=zeros(md.numberofgrids,1);
-			vz(md.results.DiagnosticAnalysis.Vz.index)=md.results.DiagnosticAnalysis.Vz.value;
+			vz(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value;
 			field_values={vy,vz};
 		else
 			field_names={'Vel'};
 			vel=zeros(md.numberofgrids,1);
-			vel(md.results.DiagnosticAnalysis.Vel.index)=md.results.DiagnosticAnalysis.Vel.value;
+			vel(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 			field_values={vel};
 		end
@@ -52,7 +52,7 @@
 	field_names={'Temperature','Melting'};
 	temperature=zeros(md.numberofgrids,1);
-	temperature(md.results.ThermalAnalysis(end).Temperature.index)=md.results.ThermalAnalysis(end).Temperature.value;
+	temperature(md.results.ThermalSolution(end).Temperature.index)=md.results.ThermalSolution(end).Temperature.value;
 	melting=zeros(md.numberofgrids,1);
-	melting(md.results.ThermalAnalysis(end).Melting.index)=md.results.ThermalAnalysis(end).Melting.value;
+	melting(md.results.ThermalSolution(end).Melting.index)=md.results.ThermalSolution(end).Melting.value;
 	field_values={temperature,melting};
 
@@ -60,5 +60,5 @@
 	field_names={'Thickness'};
 	thickness=zeros(md.numberofgrids,1);
-	thickness(md.results.PrognosticAnalysis.Thickness.index)=md.results.PrognosticAnalysis.Thickness.value;
+	thickness(md.results.PrognosticSolution.Thickness.index)=md.results.PrognosticSolution.Thickness.value;
 	field_values={thickness};
 
@@ -66,5 +66,5 @@
 	field_names={'Thickness'};
 	thickness=zeros(md.numberofgrids,1);
-	thickness(md.results.Prognostic2Analysis.Thickness.index)=md.results.Prognostic2Analysis.Thickness.value;
+	thickness(md.results.Prognostic2Solution.Thickness.index)=md.results.Prognostic2Solution.Thickness.value;
 	field_values={thickness};
 
@@ -72,11 +72,11 @@
 	field_names={'Vel','Thickness','Surface','Bed'};
 	vel=zeros(md.numberofgrids,1);
-	vel(md.results.Transient2DAnalysis(end).Vel.index)=md.results.Transient2DAnalysis(end).Vel.value;
+	vel(md.results.Transient2DSolution(end).Vel.index)=md.results.Transient2DSolution(end).Vel.value;
 	thickness=zeros(md.numberofgrids,1);
-	thickness(md.results.Transient2DAnalysis(end).Thickness.index)=md.results.Transient2DAnalysis(end).Thickness.value;
+	thickness(md.results.Transient2DSolution(end).Thickness.index)=md.results.Transient2DSolution(end).Thickness.value;
 	surface=zeros(md.numberofgrids,1);
-	surface(md.results.Transient2DAnalysis(end).Surface.index)=md.results.Transient2DAnalysis(end).Surface.value;
+	surface(md.results.Transient2DSolution(end).Surface.index)=md.results.Transient2DSolution(end).Surface.value;
 	bed=zeros(md.numberofgrids,1);
-	bed(md.results.Transient2DAnalysis(end).Bed.index)=md.results.Transient2DAnalysis(end).Bed.value;
+	bed(md.results.Transient2DSolution(end).Bed.index)=md.results.Transient2DSolution(end).Bed.value;
 	field_values={vel,thickness,surface,bed};
 
@@ -84,17 +84,17 @@
 	field_names={'Vel','Pressure','Temperature','Melting','Thickness','Surface','Bed'};
 	vel=zeros(md.numberofgrids,1);
-	vel(md.results.Transient3DAnalysis(end).Vel.index)=md.results.Transient3DAnalysis(end).Vel.value;
+	vel(md.results.Transient3DSolution(end).Vel.index)=md.results.Transient3DSolution(end).Vel.value;
 	pressure=zeros(md.numberofgrids,1);
-	pressure(md.results.Transient3DAnalysis(end).Pressure.index)=md.results.Transient3DAnalysis(end).Pressure.value;
+	pressure(md.results.Transient3DSolution(end).Pressure.index)=md.results.Transient3DSolution(end).Pressure.value;
 	temperature=zeros(md.numberofgrids,1);
-	temperature(md.results.Transient3DAnalysis(end).Temperature.index)=md.results.Transient3DAnalysis(end).Temperature.value;
+	temperature(md.results.Transient3DSolution(end).Temperature.index)=md.results.Transient3DSolution(end).Temperature.value;
 	melting=zeros(md.numberofgrids,1);
-	melting(md.results.Transient3DAnalysis(end).Melting.index)=md.results.Transient3DAnalysis(end).Melting.value;
+	melting(md.results.Transient3DSolution(end).Melting.index)=md.results.Transient3DSolution(end).Melting.value;
 	thickness=zeros(md.numberofgrids,1);
-	thickness(md.results.Transient3DAnalysis(end).Thickness.index)=md.results.Transient3DAnalysis(end).Thickness.value;
+	thickness(md.results.Transient3DSolution(end).Thickness.index)=md.results.Transient3DSolution(end).Thickness.value;
 	surface=zeros(md.numberofgrids,1);
-	surface(md.results.Transient3DAnalysis(end).Surface.index)=md.results.Transient3DAnalysis(end).Surface.value;
+	surface(md.results.Transient3DSolution(end).Surface.index)=md.results.Transient3DSolution(end).Surface.value;
 	bed=zeros(md.numberofgrids,1);
-	bed(md.results.Transient3DAnalysis(end).Bed.index)=md.results.Transient3DAnalysis(end).Bed.value;
+	bed(md.results.Transient3DSolution(end).Bed.index)=md.results.Transient3DSolution(end).Bed.value;
 	field_values={vel,pressure,temperature,melting,thickness,surface,bed};
 
@@ -103,18 +103,18 @@
 		field_names={'Vel','J','parameter'};
 		vel=zeros(md.numberofgrids,1);
-		vel(md.results.SteadystateAnalysis.Vel.index)=md.results.SteadystateAnalysis.Vel.value;
+		vel(md.results.SteadystateSolution.Vel.index)=md.results.SteadystateSolution.Vel.value;
 		j=zeros(md.numberofgrids,1);
-		j(md.results.SteadystateAnalysis.J.index)=md.results.SteadystateAnalysis.J.value;
+		j(md.results.SteadystateSolution.J.index)=md.results.SteadystateSolution.J.value;
 		parameter=zeros(md.numberofgrids,1);
-		parameter(md.results.SteadystateAnalysis.parameter.index)=md.results.SteadystateAnalysis.parameter.value;
+		parameter(md.results.SteadystateSolution.parameter.index)=md.results.SteadystateSolution.parameter.value;
 		field_values={vel,j,parameter};
 	else
 		field_names={'Vel','Pressure','Temperature'};
 		vel=zeros(md.numberofgrids,1);
-		vel(md.results.SteadystateAnalysis.Vel.index)=md.results.SteadystateAnalysis.Vel.value;
+		vel(md.results.SteadystateSolution.Vel.index)=md.results.SteadystateSolution.Vel.value;
 		presure=zeros(md.numberofgrids,1);
-		presure(md.results.SteadystateAnalysis.Pressure.index)=md.results.SteadystateAnalysis.Pressure.value;
+		presure(md.results.SteadystateSolution.Pressure.index)=md.results.SteadystateSolution.Pressure.value;
 		temperature=zeros(md.numberofgrids,1);
-		temperature(md.results.SteadystateAnalysis.Temperature.index)=md.results.SteadystateAnalysis.Temperature.value;
+		temperature(md.results.SteadystateSolution.Temperature.index)=md.results.SteadystateSolution.Temperature.value;
 		field_values={vel,pressure,temperature};
 	end
@@ -123,5 +123,5 @@
 	field_names={'Thickness'};
 	thickness=zeros(md.numberofgrids,1);
-	thickness(md.results.BalancedthicknessAnalysis.Thickness.index)=md.results.BalancedthicknessAnalysis.Thickness.value;
+	thickness(md.results.BalancedthicknessSolution.Thickness.index)=md.results.BalancedthicknessSolution.Thickness.value;
 	field_values={thickness};
 
@@ -129,7 +129,7 @@
 	field_names={'Slopex','Slopey'};
 	slopex=zeros(md.numberofgrids,1);
-	slopex(md.results.DiagnosticAnalysis.Slopex.index)=md.results.DiagnosticAnalysis.Slopex.value;
+	slopex(md.results.DiagnosticSolution.Slopex.index)=md.results.DiagnosticSolution.Slopex.value;
 	slopey=zeros(md.numberofgrids,1);
-	slopey(md.results.DiagnosticAnalysis.Slopey.index)=md.results.DiagnosticAnalysis.Slopey.value;
+	slopey(md.results.DiagnosticSolution.Slopey.index)=md.results.DiagnosticSolution.Slopey.value;
 	field_values={slopex,slopey};
 
@@ -137,7 +137,7 @@
 	field_names={'Slopex','Slopey'};
 	slopex=zeros(md.numberofgrids,1);
-	slopex(md.results.DiagnosticAnalysis.Slopex.index)=md.results.DiagnosticAnalysis.Slopex.value;
+	slopex(md.results.DiagnosticSolution.Slopex.index)=md.results.DiagnosticSolution.Slopex.value;
 	slopey=zeros(md.numberofgrids,1);
-	slopey(md.results.DiagnosticAnalysis.Slopey.index)=md.results.DiagnosticAnalysis.Slopey.value;
+	slopey(md.results.DiagnosticSolution.Slopey.index)=md.results.DiagnosticSolution.Slopey.value;
 	field_values={slopex,slopey};