Changeset 4242 for anuga_work/development/dam_test_from_brad_2007
- Timestamp:
- Feb 7, 2007, 6:35:48 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/dam_test_from_brad_2007/dam_sample.py
r4219 r4242 7 7 #------------------------------------------------------------------------------ 8 8 9 import os 10 import time 11 from shutil import copy 12 from os import mkdir, access, F_OK, sep, environ, getenv, getcwd 13 from os.path import expanduser 14 import sys 9 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross 15 10 16 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross17 11 from anuga.shallow_water import Domain 18 12 from anuga.shallow_water import Reflective_boundary 19 13 from anuga.shallow_water import Dirichlet_boundary 20 from anuga.utilities.polygon import read_polygon, plot_polygons, \ 21 polygon_area, is_inside_polygon 14 15 from anuga.utilities.polygon import read_polygon, plot_polygons 22 16 from anuga.utilities.polygon import Polygon_function 23 from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files24 from anuga.geospatial_data.geospatial_data import *25 17 26 18 #------------------------------------------------------------------------------ … … 38 30 domain.set_store_vertices_uniquely(True) # Look at the 'real' triangles 39 31 domain.beta_h = 0 40 #domain.alpha_balance = 10000000041 42 43 32 44 33 … … 51 40 #dam = [[0,0], [6,0], [6,10], [0,10], [0,0]] 52 41 poly_dam = read_polygon('dam.csv') 53 #print poly_dam 42 print poly_dam 43 #plot_polygons([poly_dam]) 44 54 45 55 46 def dam1(x,y): … … 72 63 73 64 return z 74 """75 # maunal definition76 65 77 domain.set_quantity('elevation', water) # Use function for elevation 78 domain.set_quantity('friction', 0.01) # Constant friction 79 domain.set_quantity('stage', expression='elevation') # take elevation 66 ##manual definition 67 # 68 #domain.set_quantity('elevation', water) # Use function for elevation 69 #domain.set_quantity('friction', 0.01) # Constant friction 70 #domain.set_quantity('stage', expression='elevation') # take elevation 71 #domain.set_quantity('elevation', dam1) # Use function for elevation 72 73 74 # Assign values to quantities 80 75 domain.set_quantity('elevation', dam1) # Use function for elevation 81 """ 82 83 # polygon definition 84 domain.set_quantity('elevation', dam1) # Use function for elevation 85 domain.set_quantity('friction', 0.01) # Constant friction 76 domain.set_quantity('friction', 0.01) # Constant friction 86 77 domain.set_quantity('stage', Polygon_function( [(study, -1),(poly_dam, 2)] )) 87 78 … … 91 82 #------------------------------------------------------------------------------ 92 83 Bi = Dirichlet_boundary([1.5, 0, 0]) # Inflow 93 Bo = Dirichlet_boundary([-20/10, 0, 0]) 84 Bo = Dirichlet_boundary([-20/10, 0, 0]) # Outflow 94 85 Br = Reflective_boundary(domain) # Solid reflective wall 95 86 96 #domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br}) 97 domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br}) 87 domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br}) 98 88 99 89 … … 101 91 # Evolve system through time 102 92 #------------------------------------------------------------------------------ 103 for t in domain.evolve(yieldstep = 0.2, finaltime = 10.0):104 domain.write_time( )93 for t in domain.evolve(yieldstep = 0.2, finaltime = 20): 94 domain.write_time(track_speeds=True) 105 95
Note: See TracChangeset
for help on using the changeset viewer.