Changeset 21632

Timestamp:

03/24/17 01:21:32 (8 years ago)

Author:

bdef

Message:

CHG:Removing prints in 1601

File:

• issm/trunk-jpl/test/NightlyRun/test1601.py (modified) (3 diffs)

issm/trunk-jpl/test/NightlyRun/test1601.py

@@ -15 +15 @@
 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 +23 @@
 md=solve(md,'Stressbalance')
 vel0=md.results.StressbalanceSolution.Vel
-
 theta=30.*np.pi/180.
 x=md.mesh.x
@@ -35 +34 @@
 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))