Index: /issm/trunk-jpl/test/NightlyRun/test331.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test331.m	(revision 16347)
+++ /issm/trunk-jpl/test/NightlyRun/test331.m	(revision 16348)
@@ -6,5 +6,5 @@
 md.initialization.watercolumn=10^-3*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcwatercolumn=NaN*ones(md.mesh.numberofvertices,1);
-pos=find(md.mesh.y==0);
+pos=find(md.mesh.y==0 | md.mesh.x==0);
 md.hydrology.spcwatercolumn(pos)=2*10^-3;
 md=solve(md,HydrologySolutionEnum());
Index: /issm/trunk-jpl/test/NightlyRun/test331.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test331.py	(revision 16347)
+++ /issm/trunk-jpl/test/NightlyRun/test331.py	(revision 16348)
@@ -17,5 +17,5 @@
 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]
+pos=numpy.nonzero(numpy.logical_or(md.mesh.y==0.,md.mesh.x==0.))[0]
 md.hydrology.spcwatercolumn[pos]=2e-3
 md=solve(md,HydrologySolutionEnum())