Index: /issm/trunk-jpl/test/NightlyRun/test232.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test232.py	(revision 21303)
+++ /issm/trunk-jpl/test/NightlyRun/test232.py	(revision 21304)
@@ -15,5 +15,5 @@
 md=setflowequation(md,'SSA','all')
 md.cluster=generic('name',oshostname(),'np',3)
-md.thermal.spctemperature=numpy.vstack((numpy.hstack((md.thermal.spctemperature, md.thermal.spctemperature+5., md.thermal.spctemperature+10., md.thermal.spctemperature+15.)), [1.5,2.5,3.5,4.]))
+md.thermal.spctemperature=numpy.vstack((numpy.vstack((md.thermal.spctemperature, md.thermal.spctemperature+5., md.thermal.spctemperature+10., md.thermal.spctemperature+15.)).T, [1.5,2.5,3.5,4.]))
 md.timestepping.time_step=1.
 md.timestepping.final_time=4.
@@ -26,14 +26,15 @@
 
 #Fields and tolerances to track changes
-field_names     =['Temperature1','BasalforcingsGroundediceMeltingRate1','Temperature2','BasalforcingsGroundediceMeltingRate2','Temperature3','BasalforcingsGroundediceMeltingRate3','Temperature4','BasalforcingsGroundediceMeltingRate4']
+field_names     =['Temperature1','BasalforcingsGroundediceMeltingRate1',
+									'Temperature2','BasalforcingsGroundediceMeltingRate2',
+									'Temperature3','BasalforcingsGroundediceMeltingRate3',
+									'Temperature4','BasalforcingsGroundediceMeltingRate4']
 field_tolerances=[1e-13,1e-6,1e-13,1e-6,1e-13,1e-6,1e-13,1e-6]
-field_values=[\
-	md.results.TransientSolution[0].Temperature,\
-	md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[1].Temperature,\
-	md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[2].Temperature,\
-	md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[3].Temperature,\
-	md.results.TransientSolution[3].BasalforcingsGroundediceMeltingRate,\
-	]
+field_values=[md.results.TransientSolution[0].Temperature,
+							md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[1].Temperature,
+							md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[2].Temperature,
+							md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[3].Temperature,
+							md.results.TransientSolution[3].BasalforcingsGroundediceMeltingRate]
Index: /issm/trunk-jpl/test/NightlyRun/test233.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test233.py	(revision 21303)
+++ /issm/trunk-jpl/test/NightlyRun/test233.py	(revision 21304)
@@ -15,5 +15,5 @@
 md=setflowequation(md,'HO','all')
 md.cluster=generic('name',oshostname(),'np',3)
-md.thermal.spctemperature=numpy.vstack((numpy.hstack((md.thermal.spctemperature, md.thermal.spctemperature+5.)), [1.,2.]))
+md.thermal.spctemperature=numpy.vstack((numpy.vstack((md.thermal.spctemperature, md.thermal.spctemperature+5.)).T, [1.,2.]))
 md.timestepping.time_step=0.5
 md.timestepping.final_time=2.
@@ -21,52 +21,50 @@
 
 #Fields and tolerances to track changes
-field_names     =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', \
-	'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2', \
-	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3', \
-	'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsGroundediceMeltingRate4']
-field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, \
-	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, \
-	1e-08,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, \
-	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06]
-field_values=[\
-	md.results.TransientSolution[0].Vx,\
-	md.results.TransientSolution[0].Vy,\
-	md.results.TransientSolution[0].Vz,\
-	md.results.TransientSolution[0].Vel,\
-	md.results.TransientSolution[0].Pressure,\
-	md.results.TransientSolution[0].Base,\
-	md.results.TransientSolution[0].Surface,\
-	md.results.TransientSolution[0].Thickness,\
-	md.results.TransientSolution[0].Temperature,\
-	md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[1].Vx,\
-	md.results.TransientSolution[1].Vy,\
-	md.results.TransientSolution[1].Vz,\
-	md.results.TransientSolution[1].Vel,\
-	md.results.TransientSolution[1].Pressure,\
-	md.results.TransientSolution[1].Base,\
-	md.results.TransientSolution[1].Surface,\
-	md.results.TransientSolution[1].Thickness,\
-	md.results.TransientSolution[1].Temperature,\
-	md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[2].Vx,\
-	md.results.TransientSolution[2].Vy,\
-	md.results.TransientSolution[2].Vz,\
-	md.results.TransientSolution[2].Vel,\
-	md.results.TransientSolution[2].Pressure,\
-	md.results.TransientSolution[2].Base,\
-	md.results.TransientSolution[2].Surface,\
-	md.results.TransientSolution[2].Thickness,\
-	md.results.TransientSolution[2].Temperature,\
-	md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,\
-	md.results.TransientSolution[3].Vx,\
-	md.results.TransientSolution[3].Vy,\
-	md.results.TransientSolution[3].Vz,\
-	md.results.TransientSolution[3].Vel,\
-	md.results.TransientSolution[3].Pressure,\
-	md.results.TransientSolution[3].Base,\
-	md.results.TransientSolution[3].Surface,\
-	md.results.TransientSolution[3].Thickness,\
-	md.results.TransientSolution[3].Temperature,\
-	md.results.TransientSolution[3].BasalforcingsGroundediceMeltingRate,\
-	]
+field_names     =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1',
+									'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2',
+									'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3',
+									'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsGroundediceMeltingRate4']
+field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,
+									1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06,
+									1e-08,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06,
+									1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06]
+field_values=[md.results.TransientSolution[0].Vx,
+							md.results.TransientSolution[0].Vy,
+							md.results.TransientSolution[0].Vz,
+							md.results.TransientSolution[0].Vel,
+							md.results.TransientSolution[0].Pressure,
+							md.results.TransientSolution[0].Base,
+							md.results.TransientSolution[0].Surface,
+							md.results.TransientSolution[0].Thickness,
+							md.results.TransientSolution[0].Temperature,
+							md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[1].Vx,
+							md.results.TransientSolution[1].Vy,
+							md.results.TransientSolution[1].Vz,
+							md.results.TransientSolution[1].Vel,
+							md.results.TransientSolution[1].Pressure,
+							md.results.TransientSolution[1].Base,
+							md.results.TransientSolution[1].Surface,
+							md.results.TransientSolution[1].Thickness,
+							md.results.TransientSolution[1].Temperature,
+							md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[2].Vx,
+							md.results.TransientSolution[2].Vy,
+							md.results.TransientSolution[2].Vz,
+							md.results.TransientSolution[2].Vel,
+							md.results.TransientSolution[2].Pressure,
+							md.results.TransientSolution[2].Base,
+							md.results.TransientSolution[2].Surface,
+							md.results.TransientSolution[2].Thickness,
+							md.results.TransientSolution[2].Temperature,
+							md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,
+							md.results.TransientSolution[3].Vx,
+							md.results.TransientSolution[3].Vy,
+							md.results.TransientSolution[3].Vz,
+							md.results.TransientSolution[3].Vel,
+							md.results.TransientSolution[3].Pressure,
+							md.results.TransientSolution[3].Base,
+							md.results.TransientSolution[3].Surface,
+							md.results.TransientSolution[3].Thickness,
+							md.results.TransientSolution[3].Temperature,
+							md.results.TransientSolution[3].BasalforcingsGroundediceMeltingRate]