Changeset 9022 for trunk/anuga_work/development/gareth/tests/parabolic

Timestamp:

Nov 13, 2013, 9:44:05 PM (11 years ago)

Author:

davies

Message:

Removing most significant velocity artefacts from dev_audusse

Location:

trunk/anuga_work/development/gareth/tests/parabolic

Files:

• parabolaplot.py (modified) (2 diffs)
• parabolic.py (modified) (1 diff)

trunk/anuga_work/development/gareth/tests/parabolic/parabolaplot.py

@@ -31 +31 @@
         pyplot.draw()
         #pyplot.plot( (0,1),(0,0), 'r' )
-        pyplot.plot(p2.x[v],p2.elev[1,v]) 
+        #pyplot.plot(p2.x[v],p2.elev[v]) 
         pyplot.title(str(i)+'/200') # : velocity does not converge to zero' )
         pyplot.xlabel('x')
@@ -63 +63 @@
 w = D0 + 2*A*D0/(L**2)*numpy.cos(omega*p2.time)*( (p2.x[ppp]-4.*L/2.) -A/2.*numpy.cos(omega*p2.time))
 w2 = D0 + 2*A*D0/(L**2)*numpy.cos(omega*p2.time)*( (p2.x[ppp2]-4.*L/2.) -A/2.*numpy.cos(omega*p2.time))
-w2 = w2*(w2>p2.elev[1,ppp2])+p2.elev[1,ppp2]*(w2<=p2.elev[1,ppp2])
+w2 = w2*(w2>p2.elev[ppp2])+p2.elev[ppp2]*(w2<=p2.elev[ppp2])
 
 pyplot.clf()

trunk/anuga_work/development/gareth/tests/parabolic/parabolic.py

@@ -23 +23 @@
 domain.set_name('parabola_v2')                         # Output to file runup.sww
 domain.set_datadir('.')                          # Use current folder
+domain.set_store_centroids(True)
 #domain.set_minimum_allowed_height(0.01)