Changeset 5162 for anuga_work/development

Timestamp:

Mar 11, 2008, 8:43:22 PM (17 years ago)

Author:

steve

Message:

Updated some methods for quantity. Looks like we can use old
limiting system with larger values of beta.

Location:

anuga_work/development

Files:

• analytical_solutions/oblique_shock.py (modified) (5 diffs)
• near_shore_PMD/run_beach.py (modified) (2 diffs)

anuga_work/development/analytical_solutions/oblique_shock.py

@@ -29 +29 @@
 
 
-from shallow_water import Domain, Constant_height
-from shallow_water import Transmissive_boundary, Reflective_boundary,\
+from anuga.shallow_water import Domain
+from anuga.shallow_water import Transmissive_boundary, Reflective_boundary,\
      Dirichlet_boundary
-
+from anuga.visualiser import RealtimeVisualiser
 from math import sqrt, cos, sin, pi
-from mesh_factory import oblique
+from mesh_factory import oblique_cross
 
 
@@ -44 +44 @@
 leny = 30.
 lenx = 40.
-n = 50
-m = 60
+n = 100
+m = 120
 
-points, elements, boundary = oblique(m, n, lenx, leny)
+points, elements, boundary = oblique_cross(m, n, lenx, leny)
 domain = Domain(points, elements, boundary)
 
@@ -54 +54 @@
 
 # Store output
-domain.store=True
+#domain.store=True
 
 # Output format
-domain.format="sww" #NET.CDF binary format
+#domain.format="sww" #NET.CDF binary format
                     # "dat" for ASCII
 
@@ -87 +87 @@
 #Initial condition
 h = 0.5
-domain.set_quantity('stage', Constant_height(x_slope, h) )
+domain.set_quantity('stage', expression = 'elevation + %d'%h )
 
+#---------------------------------
+# Setup visualization
+#---------------------------------
+vis = RealtimeVisualiser(domain)
+vis.render_quantity_height("elevation", dynamic=False)
+vis.render_quantity_height("stage", zScale=10, dynamic=True)
+vis.colour_height_quantity('stage', (0.75, 0.5, 0.5))
+vis.start()
 
 ######################
@@ -96 +104 @@
 for t in domain.evolve(yieldstep = 0.5, finaltime = 50):
     domain.write_time()
+    vis.update()
 
 print 'That took %.2f seconds' %(time.time()-t0)
+
+
+vis.evolveFinished()
+vis.join()
 
 #FIXME: Compute average water depth on either side of shock and compare

anuga_work/development/near_shore_PMD/run_beach.py

@@ -167 +167 @@
     domain.set_boundary( {'wall': Br, 'wave': Bwp_velocity} )
 
+
     #-------------------------------------------------------------------------
     # Evolve system through time
@@ -174 +175 @@
     for t in domain.evolve(yieldstep, finaltime):    
         domain.write_time()
+
     print 'That took %.2f seconds' %(time.time()-t0)
     print 'finished'