Index: /issm/trunk-jpl/test/NightlyRun/test1601.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1601.py	(revision 21631)
+++ /issm/trunk-jpl/test/NightlyRun/test1601.py	(revision 21632)
@@ -15,5 +15,5 @@
 md=parameterize(md,'../Par/SquareShelf.py')
 md=setflowequation(md,'SSA','all')
-md.stressbalance.spcvx[np.nonzero(md.mesh.y>0.)]=np.nan
+md.stressbalance.spcvx[np.where(md.mesh.y>0.)]=np.nan
 md.initialization.vx[:]=0.
 md.initialization.vy[:]=0.
@@ -23,5 +23,4 @@
 md=solve(md,'Stressbalance')
 vel0=md.results.StressbalanceSolution.Vel
-
 theta=30.*np.pi/180.
 x=md.mesh.x
@@ -35,5 +34,4 @@
 md=solve(md,'Stressbalance')
 vel1=md.results.StressbalanceSolution.Vel
-
 #plotmodel(md,'data',vel0,'data',vel1,'data',vel1-vel0,'title','Cartesian CS','title','Rotated CS','title','difference')
 print "Error between Cartesian and rotated CS: %g" % (np.max(np.abs(vel0-vel1))/(np.max(np.abs(vel0))+sys.float_info.epsilon))