Changeset 4242 for anuga_work/development/dam_test_from_brad_2007

Timestamp:

Feb 7, 2007, 6:35:48 PM (18 years ago)

Author:

ole

Message:

Work on Brad's dam example

File:

• anuga_work/development/dam_test_from_brad_2007/dam_sample.py (modified) (6 diffs)

anuga_work/development/dam_test_from_brad_2007/dam_sample.py

@@ -7 +7 @@
 #------------------------------------------------------------------------------
 
-import os
-import time
-from shutil import copy
-from os import mkdir, access, F_OK, sep, environ, getenv, getcwd
-from os.path import expanduser
-import sys
+from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
 
-from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
 from anuga.shallow_water import Domain
 from anuga.shallow_water import Reflective_boundary
 from anuga.shallow_water import Dirichlet_boundary
-from anuga.utilities.polygon import read_polygon, plot_polygons, \
-                                    polygon_area, is_inside_polygon
+
+from anuga.utilities.polygon import read_polygon, plot_polygons
 from anuga.utilities.polygon import Polygon_function
-from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files 
-from anuga.geospatial_data.geospatial_data import *
 
 #------------------------------------------------------------------------------
@@ -38 +30 @@
 domain.set_store_vertices_uniquely(True)      # Look at the 'real' triangles
 domain.beta_h = 0
-#domain.alpha_balance = 100000000
-
-
 
 
@@ -51 +40 @@
 #dam = [[0,0], [6,0], [6,10], [0,10], [0,0]]
 poly_dam = read_polygon('dam.csv')
-#print poly_dam
+print poly_dam
+#plot_polygons([poly_dam])
+
 
 def dam1(x,y):
@@ -72 +63 @@
             
     return z
-"""
-# maunal definition
 
-domain.set_quantity('elevation', water)                # Use function for elevation
-domain.set_quantity('friction', 0.01)                       # Constant friction 
-domain.set_quantity('stage', expression='elevation')        # take elevation
+##manual definition
+#
+#domain.set_quantity('elevation', water)                # Use function for elevation
+#domain.set_quantity('friction', 0.01)                       # Constant friction 
+#domain.set_quantity('stage', expression='elevation')        # take elevation
+#domain.set_quantity('elevation', dam1)                # Use function for elevation
+
+
+# Assign values to quantities
 domain.set_quantity('elevation', dam1)                # Use function for elevation
-"""
-
-# polygon definition
-domain.set_quantity('elevation', dam1)                # Use function for elevation
-domain.set_quantity('friction', 0.01)                       # Constant friction 
+domain.set_quantity('friction', 0.01)                 # Constant friction 
 domain.set_quantity('stage', Polygon_function( [(study, -1),(poly_dam, 2)] )) 
 
@@ -91 +82 @@
 #------------------------------------------------------------------------------
 Bi = Dirichlet_boundary([1.5, 0, 0])         # Inflow
-Bo = Dirichlet_boundary([-20/10, 0, 0])         # Outflow
+Bo = Dirichlet_boundary([-20/10, 0, 0])      # Outflow
 Br = Reflective_boundary(domain)             # Solid reflective wall
 
-#domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br})
-domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
+domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br})
 
 
@@ -101 +91 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-for t in domain.evolve(yieldstep = 0.2, finaltime = 10.0):
-    domain.write_time()
+for t in domain.evolve(yieldstep = 0.2, finaltime = 20):
+    domain.write_time(track_speeds=True)