Changeset 6713 for anuga_work/development/anuga_1d/discontinuous_z.py

Timestamp:

Apr 3, 2009, 12:34:47 PM (15 years ago)

Author:

steve

Message:

Fixed pylab so that it would work on bogong, by flattening the arrays for plotting

File:

• anuga_work/development/anuga_1d/discontinuous_z.py (modified) (3 diffs)

anuga_work/development/anuga_1d/discontinuous_z.py

@@ -19 +19 @@
 domain.set_timestepping_method('rk2')
 domain.cfl = 1.0
-domain.limiter = "minmod"
+domain.limiter = "vanleer"
 
 def stage_flat(x):
@@ -85 +85 @@
 domain.set_quantity('elevation',elevation_arbitrary)
 domain.set_boundary({'exterior':Reflective_boundary(domain)})
-X=domain.vertices
-C=domain.centroids
+X=domain.vertices.flat
+C=domain.centroids.flat
 
 t0=time.time()
@@ -95 +95 @@
         print "t=",t
         print "integral", domain.quantities['height'].get_integral()
-        StageQ=domain.quantities['stage'].vertex_values
-        MomentumQ=domain.quantities['xmomentum'].vertex_values
-        ElevationQ=domain.quantities['elevation'].vertex_values
-        VelocityQ=domain.quantities['velocity'].vertex_values
+        StageQ=domain.quantities['stage'].vertex_values.flat
+        MomentumQ=domain.quantities['xmomentum'].vertex_values.flat
+        ElevationQ=domain.quantities['elevation'].vertex_values.flat
+        VelocityQ=domain.quantities['velocity'].vertex_values.flat
 
         hold(False)