Changeset 7907

Timestamp:

Aug 2, 2010, 2:11:03 PM (14 years ago)

Author:

steve

Message:

Problem with plotting vertices. Had to use flat structure

File:

• trunk/anuga_work/development/sudi/sw_1d/shock_detector_shv/dam_h_2sides_RUN.py (modified) (2 diffs)

trunk/anuga_work/development/sudi/sw_1d/shock_detector_shv/dam_h_2sides_RUN.py

@@ -57 +57 @@
 C=domain.centroids
 import time
-yieldstep=finaltime=1.0
+yieldstep=finaltime=5.0
 t0=time.time()
 i=1
 
-while finaltime < 1.01:
+while finaltime < 5.01:
     for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
         domain.write_time()
@@ -78 +78 @@
     #print "u_error %.10f" %(u_error)
     #print 'That took %.2f seconds' %(time.time()-t0)
-    X = domain.vertices
-    StageQ = domain.quantities['stage'].vertex_values
-    XmomQ = domain.quantities['xmomentum'].vertex_values
-    VelQ = domain.quantities['velocity'].vertex_values
-    BedQ = domain.quantities['elevation'].vertex_values
+    X = domain.vertices.flat
+    StageQ = domain.quantities['stage'].vertex_values.flat
+    XmomQ = domain.quantities['xmomentum'].vertex_values.flat
+    VelQ = domain.quantities['velocity'].vertex_values.flat
+    BedQ = domain.quantities['elevation'].vertex_values.flat
     #h, uh, u = analytical_sol(X.flat, domain.time)