Changeset 3775

Timestamp:

Oct 13, 2006, 2:16:35 PM (18 years ago)

Author:

ole

Message:

Work on simple parallel test

File:

• anuga_core/source/anuga_parallel/test_parallel_sw_runup.py (modified) (1 diff)

anuga_core/source/anuga_parallel/test_parallel_sw_runup.py

@@ -68 +68 @@
 #------------------------------------------------------------------------------
 
-for t in domain.evolve(yieldstep = 0.1, finaltime = 10.0):
+interpolation_points = [[0.4,0.5], [0.6,0.5], [0.8,0.5], [0.9,0.5]]
+gauge_values = []
+for _ in interpolation_points:
+    gauge_values.append([])
+
+time = []
+
+for t in domain.evolve(yieldstep = 0.1, finaltime = 5.0):
     domain.write_time()
+
     
+    # Record time series at known points
+    time.append(domain.get_time())
+    
+    stage = domain.get_quantity('stage')
+    w = stage.get_values(interpolation_points=interpolation_points)
+    
+    for i, _ in enumerate(interpolation_points):
+        gauge_values[i].append(w[i])
 
+
+print
+print time
+print
+for i, _ in enumerate(interpolation_points):
+    print i, gauge_values[i]
+    print 
+
+