Index: /issm/trunk-jpl/test/NightlyRun/test236.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test236.py	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test236.py	(revision 24296)
@@ -69,6 +69,4 @@
 
 
-#
-md.transient.requested_outputs = ['default', 'SmbMonthlytemperatures']
 md = setflowequation(md, 'SSA', 'all')
 md.cluster = generic('name', gethostname(), 'np', 3)
@@ -76,7 +74,7 @@
 
 #Fields and tolerances to track changes
-field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMonthlytemperatures1', 'SmbMassBalance1',
-               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMonthlytemperatures2', 'SmbMassBalance2',
-               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMonthlytemperatures3', 'SmbMassBalance3']
+field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMassBalance1',
+               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMassBalance2',
+               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMassBalance3']
 field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
                     1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
@@ -89,5 +87,4 @@
                 md.results.TransientSolution[0].Surface,
                 md.results.TransientSolution[0].Thickness,
-                md.results.TransientSolution[0].SmbMonthlytemperatures,
                 md.results.TransientSolution[0].SmbMassBalance,
                 md.results.TransientSolution[1].Vx,
@@ -98,5 +95,4 @@
                 md.results.TransientSolution[1].Surface,
                 md.results.TransientSolution[1].Thickness,
-                md.results.TransientSolution[1].SmbMonthlytemperatures,
                 md.results.TransientSolution[1].SmbMassBalance,
                 md.results.TransientSolution[2].Vx,
@@ -107,4 +103,3 @@
                 md.results.TransientSolution[2].Surface,
                 md.results.TransientSolution[2].Thickness,
-                md.results.TransientSolution[2].SmbMonthlytemperatures,
                 md.results.TransientSolution[2].SmbMassBalance]
Index: /issm/trunk-jpl/test/NightlyRun/test237.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.m	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test237.m	(revision 24296)
@@ -55,6 +55,4 @@
 end
 
-%
-md.transient.requested_outputs={'default','SmbMonthlytemperatures'};
 md=extrude(md,3,1);
 
@@ -64,7 +62,7 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','SmbMonthlytemperatures1','SmbMassBalance1',...
-						'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','SmbMonthlytemperatures2','SmbMassBalance2',...
-						'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3','SmbMonthlytemperatures3','SmbMassBalance3'};
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','SmbMassBalance1',...
+						'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','SmbMassBalance2',...
+						'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3','SmbMassBalance3'};
 field_tolerances={...
 	1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,1e-8,1e-8,1e-13,1e-8,1e-13,1e-13,...
@@ -82,5 +80,4 @@
 	(md.results.TransientSolution(1).Temperature),...
 	(md.results.TransientSolution(1).BasalforcingsGroundediceMeltingRate),...
-	(md.results.TransientSolution(1).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(1).SmbMassBalance),...
 	(md.results.TransientSolution(2).Vx),...
@@ -94,5 +91,4 @@
 	(md.results.TransientSolution(2).Temperature),...
 	(md.results.TransientSolution(2).BasalforcingsGroundediceMeltingRate),...
-	(md.results.TransientSolution(2).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(2).SmbMassBalance),...
 	(md.results.TransientSolution(3).Vx),...
@@ -106,5 +102,4 @@
 	(md.results.TransientSolution(3).Temperature),...
 	(md.results.TransientSolution(3).BasalforcingsGroundediceMeltingRate),...
-	(md.results.TransientSolution(3).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(3).SmbMassBalance),...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test237.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 24296)
@@ -69,6 +69,4 @@
     md.smb.sealev[1, iint] = (float(iint)) * 20
 
-#
-md.transient.requested_outputs = ['default', 'SmbMonthlytemperatures']
 md.extrude(3, 1.)
 md = setflowequation(md, 'SSA', 'all')
@@ -77,7 +75,7 @@
 
 #Fields and tolerances to track changes
-field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'SmbMonthlytemperatures1', 'SmbMassBalance1',
-               'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'BasalforcingsGroundediceMeltingRate2', 'SmbMonthlytemperatures2', 'SmbMassBalance2',
-               'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'BasalforcingsGroundediceMeltingRate3', 'SmbMonthlytemperatures3', 'SmbMassBalance3']
+field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'SmbMassBalance1',
+               'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'BasalforcingsGroundediceMeltingRate2', 'SmbMassBalance2',
+               'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'BasalforcingsGroundediceMeltingRate3', 'SmbMassBalance3']
 field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-8, 1e-8, 1e-8, 1e-13, 1e-8, 1e-13, 1e-13,
                     1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-8, 1e-8, 1e-8, 7e-13, 2e-7, 1e-13, 1e-13,
@@ -93,5 +91,4 @@
                 md.results.TransientSolution[0].Temperature,
                 md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,
-                md.results.TransientSolution[0].SmbMonthlytemperatures,
                 md.results.TransientSolution[0].SmbMassBalance,
                 md.results.TransientSolution[1].Vx,
@@ -105,5 +102,4 @@
                 md.results.TransientSolution[1].Temperature,
                 md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,
-                md.results.TransientSolution[1].SmbMonthlytemperatures,
                 md.results.TransientSolution[1].SmbMassBalance,
                 md.results.TransientSolution[2].Vx,
@@ -117,4 +113,3 @@
                 md.results.TransientSolution[2].Temperature,
                 md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,
-                md.results.TransientSolution[2].SmbMonthlytemperatures,
                 md.results.TransientSolution[2].SmbMassBalance]
Index: /issm/trunk-jpl/test/NightlyRun/test238.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test238.m	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test238.m	(revision 24296)
@@ -44,6 +44,4 @@
 md.timestepping.final_time=60;
 
-% 
-md.transient.requested_outputs={'default','SmbMonthlytemperatures'};
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',1); % 3 for the cluster
@@ -51,7 +49,7 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',...
-	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',...
-	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3'};
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1',...
+	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2',...
+	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3'};
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
 	   1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
@@ -65,5 +63,4 @@
 	(md.results.TransientSolution(1).Surface),...
 	(md.results.TransientSolution(1).Thickness),...
-	(md.results.TransientSolution(1).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(1).SmbMassBalance),...
 	(md.results.TransientSolution(2).Vx),...
@@ -74,5 +71,4 @@
 	(md.results.TransientSolution(2).Surface),...
 	(md.results.TransientSolution(2).Thickness),...
-	(md.results.TransientSolution(2).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(2).SmbMassBalance),...
 	(md.results.TransientSolution(3).Vx),...
@@ -83,5 +79,4 @@
 	(md.results.TransientSolution(3).Surface),...
 	(md.results.TransientSolution(3).Thickness),...
-	(md.results.TransientSolution(3).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(3).SmbMassBalance),...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test238.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test238.py	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test238.py	(revision 24296)
@@ -49,6 +49,4 @@
 md.timestepping.final_time = 60
 
-#
-md.transient.requested_outputs = ['default', 'SmbMonthlytemperatures']
 md = setflowequation(md, 'SSA', 'all')
 md.cluster = generic('name', gethostname(), 'np', 3)
@@ -56,7 +54,7 @@
 
 #Fields and tolerances to track changes
-field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMonthlytemperatures1', 'SmbMassBalance1',
-               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMonthlytemperatures2', 'SmbMassBalance2',
-               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMonthlytemperatures3', 'SmbMassBalance3']
+field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMassBalance1',
+               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMassBalance2',
+               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMassBalance3']
 field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
                     1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
@@ -69,5 +67,4 @@
                 md.results.TransientSolution[0].Surface,
                 md.results.TransientSolution[0].Thickness,
-                md.results.TransientSolution[0].SmbMonthlytemperatures,
                 md.results.TransientSolution[0].SmbMassBalance,
                 md.results.TransientSolution[1].Vx,
@@ -78,5 +75,4 @@
                 md.results.TransientSolution[1].Surface,
                 md.results.TransientSolution[1].Thickness,
-                md.results.TransientSolution[1].SmbMonthlytemperatures,
                 md.results.TransientSolution[1].SmbMassBalance,
                 md.results.TransientSolution[2].Vx,
@@ -87,4 +83,3 @@
                 md.results.TransientSolution[2].Surface,
                 md.results.TransientSolution[2].Thickness,
-                md.results.TransientSolution[2].SmbMonthlytemperatures,
                 md.results.TransientSolution[2].SmbMassBalance]
Index: /issm/trunk-jpl/test/NightlyRun/test239.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test239.m	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test239.m	(revision 24296)
@@ -45,6 +45,4 @@
 md.timestepping.final_time=2;
 
-% 
-md.transient.requested_outputs={'default','SmbMonthlytemperatures'};
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',1); % 3 for the cluster
@@ -52,7 +50,7 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',...
-	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',...
-	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3'};
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1',...
+	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2',...
+	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3'};
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
 	   1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
@@ -66,5 +64,4 @@
 	(md.results.TransientSolution(1).Surface),...
 	(md.results.TransientSolution(1).Thickness),...
-	(md.results.TransientSolution(1).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(1).SmbMassBalance),...
 	(md.results.TransientSolution(2).Vx),...
@@ -75,5 +72,4 @@
 	(md.results.TransientSolution(2).Surface),...
 	(md.results.TransientSolution(2).Thickness),...
-	(md.results.TransientSolution(2).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(2).SmbMassBalance),...
 	(md.results.TransientSolution(3).Vx),...
@@ -84,5 +80,4 @@
 	(md.results.TransientSolution(3).Surface),...
 	(md.results.TransientSolution(3).Thickness),...
-	(md.results.TransientSolution(3).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(3).SmbMassBalance),...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test239.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 24296)
@@ -49,6 +49,4 @@
 md.timestepping.final_time = 2
 
-#
-md.transient.requested_outputs = ['default', 'SmbMonthlytemperatures']
 md = setflowequation(md, 'SSA', 'all')
 md.cluster = generic('name', gethostname(), 'np', 3)
@@ -56,7 +54,7 @@
 
 #Fields and tolerances to track changes
-field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMonthlytemperatures1', 'SmbMassBalance1',
-               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMonthlytemperatures2', 'SmbMassBalance2',
-               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMonthlytemperatures3', 'SmbMassBalance3']
+field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMassBalance1',
+               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMassBalance2',
+               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMassBalance3']
 field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
                     1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
@@ -69,5 +67,4 @@
                 md.results.TransientSolution[0].Surface,
                 md.results.TransientSolution[0].Thickness,
-                md.results.TransientSolution[0].SmbMonthlytemperatures,
                 md.results.TransientSolution[0].SmbMassBalance,
                 md.results.TransientSolution[1].Vx,
@@ -78,5 +75,4 @@
                 md.results.TransientSolution[1].Surface,
                 md.results.TransientSolution[1].Thickness,
-                md.results.TransientSolution[1].SmbMonthlytemperatures,
                 md.results.TransientSolution[1].SmbMassBalance,
                 md.results.TransientSolution[2].Vx,
@@ -87,4 +83,3 @@
                 md.results.TransientSolution[2].Surface,
                 md.results.TransientSolution[2].Thickness,
-                md.results.TransientSolution[2].SmbMonthlytemperatures,
                 md.results.TransientSolution[2].SmbMassBalance]
Index: /issm/trunk-jpl/test/NightlyRun/test240.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test240.m	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test240.m	(revision 24296)
@@ -45,6 +45,4 @@
 md.timestepping.interp_forcings=0;
 
-% 
-md.transient.requested_outputs={'default','SmbMonthlytemperatures'};
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',1); % 3 for the cluster
@@ -52,7 +50,7 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',...
-	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',...
-	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3'};
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1',...
+	      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2',...
+	      'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3'};
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
 	   1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
@@ -66,5 +64,4 @@
 	(md.results.TransientSolution(1).Surface),...
 	(md.results.TransientSolution(1).Thickness),...
-	(md.results.TransientSolution(1).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(1).SmbMassBalance),...
 	(md.results.TransientSolution(2).Vx),...
@@ -75,5 +72,4 @@
 	(md.results.TransientSolution(2).Surface),...
 	(md.results.TransientSolution(2).Thickness),...
-	(md.results.TransientSolution(2).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(2).SmbMassBalance),...
 	(md.results.TransientSolution(3).Vx),...
@@ -84,5 +80,4 @@
 	(md.results.TransientSolution(3).Surface),...
 	(md.results.TransientSolution(3).Thickness),...
-	(md.results.TransientSolution(3).SmbMonthlytemperatures),...
 	(md.results.TransientSolution(3).SmbMassBalance),...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test240.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test240.py	(revision 24295)
+++ /issm/trunk-jpl/test/NightlyRun/test240.py	(revision 24296)
@@ -50,6 +50,4 @@
 md.timestepping.interp_forcings = 0
 
-#
-md.transient.requested_outputs = ['default', 'SmbMonthlytemperatures']
 md = setflowequation(md, 'SSA', 'all')
 md.cluster = generic('name', gethostname(), 'np', 3)
@@ -57,7 +55,7 @@
 
 #Fields and tolerances to track changes
-field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMonthlytemperatures1', 'SmbMassBalance1',
-               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMonthlytemperatures2', 'SmbMassBalance2',
-               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMonthlytemperatures3', 'SmbMassBalance3']
+field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'SmbMassBalance1',
+               'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'SmbMassBalance2',
+               'Vx3', 'Vy3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'SmbMassBalance3']
 field_tolerances = [1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
                     1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13,
@@ -70,5 +68,4 @@
                 md.results.TransientSolution[0].Surface,
                 md.results.TransientSolution[0].Thickness,
-                md.results.TransientSolution[0].SmbMonthlytemperatures,
                 md.results.TransientSolution[0].SmbMassBalance,
                 md.results.TransientSolution[1].Vx,
@@ -79,5 +76,4 @@
                 md.results.TransientSolution[1].Surface,
                 md.results.TransientSolution[1].Thickness,
-                md.results.TransientSolution[1].SmbMonthlytemperatures,
                 md.results.TransientSolution[1].SmbMassBalance,
                 md.results.TransientSolution[2].Vx,
@@ -88,4 +84,3 @@
                 md.results.TransientSolution[2].Surface,
                 md.results.TransientSolution[2].Thickness,
-                md.results.TransientSolution[2].SmbMonthlytemperatures,
                 md.results.TransientSolution[2].SmbMassBalance]
