Index: /issm/trunk-jpl/test/NightlyRun/test331.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test331.m	(revision 14556)
+++ /issm/trunk-jpl/test/NightlyRun/test331.m	(revision 14557)
@@ -14,6 +14,6 @@
 field_tolerances={1e-09,1e-09,1e-09};
 field_values={...
-	(md.results.HydrologySolution.HydrologyWaterVx),...
-	(md.results.HydrologySolution.HydrologyWaterVy),...
-	(md.results.HydrologySolution.Watercolumn),...
+	(md.results.HydrologySolution(1).HydrologyWaterVx),...
+	(md.results.HydrologySolution(1).HydrologyWaterVy),...
+	(md.results.HydrologySolution(1).Watercolumn),...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test331.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test331.py	(revision 14556)
+++ /issm/trunk-jpl/test/NightlyRun/test331.py	(revision 14557)
@@ -14,10 +14,9 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.py')
 md=setflowequation(md,'macayeal','all')
-md.surfaceforcings.issmbgradients=1
-md.initialization.watercolumn=10^-3*numpy.ones((md.mesh.numberofvertices,1))
-md.hydrology.spcwatercolumn=float('nan')*numpy.ones(md.mesh.numberofvertices,1)
-md.nodeonicefront[pos]=1 
-md.hydrology.spcwatercolumn[pos]=2e-0-3
 md.cluster=generic('name',oshostname(),'np',3)
+md.initialization.watercolumn=1e-3*numpy.ones((md.mesh.numberofvertices,1))
+md.hydrology.spcwatercolumn=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
+pos=numpy.nonzero(md.mesh.y==0)[0]
+md.hydrology.spcwatercolumn[pos]=2e-3
 md=solve(md,HydrologySolutionEnum())
 
@@ -26,6 +25,6 @@
 field_tolerances=[1e-09,1e-09,1e-09]
 field_values=[\
-	md.results.HydrologySolution.HydrologyWaterVx,\
-	md.results.HydrologySolution.HydrologyWaterVy,\
-	md.results.HydrologySolution.Watercolumn,\
+	md.results.HydrologySolution[0].HydrologyWaterVx,\
+	md.results.HydrologySolution[0].HydrologyWaterVy,\
+	md.results.HydrologySolution[0].Watercolumn,\
 	]