Ignore:
Timestamp:
Feb 7, 2007, 6:35:48 PM (18 years ago)
Author:
ole
Message:

Work on Brad's dam example

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/dam_test_from_brad_2007/dam_sample.py

    r4219 r4242  
    77#------------------------------------------------------------------------------
    88
    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
     9from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
    1510
    16 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
    1711from anuga.shallow_water import Domain
    1812from anuga.shallow_water import Reflective_boundary
    1913from anuga.shallow_water import Dirichlet_boundary
    20 from anuga.utilities.polygon import read_polygon, plot_polygons, \
    21                                     polygon_area, is_inside_polygon
     14
     15from anuga.utilities.polygon import read_polygon, plot_polygons
    2216from anuga.utilities.polygon import Polygon_function
    23 from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files
    24 from anuga.geospatial_data.geospatial_data import *
    2517
    2618#------------------------------------------------------------------------------
     
    3830domain.set_store_vertices_uniquely(True)      # Look at the 'real' triangles
    3931domain.beta_h = 0
    40 #domain.alpha_balance = 100000000
    41 
    42 
    4332
    4433
     
    5140#dam = [[0,0], [6,0], [6,10], [0,10], [0,0]]
    5241poly_dam = read_polygon('dam.csv')
    53 #print poly_dam
     42print poly_dam
     43#plot_polygons([poly_dam])
     44
    5445
    5546def dam1(x,y):
     
    7263           
    7364    return z
    74 """
    75 # maunal definition
    7665
    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
    8075domain.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
     76domain.set_quantity('friction', 0.01)                 # Constant friction
    8677domain.set_quantity('stage', Polygon_function( [(study, -1),(poly_dam, 2)] ))
    8778
     
    9182#------------------------------------------------------------------------------
    9283Bi = Dirichlet_boundary([1.5, 0, 0])         # Inflow
    93 Bo = Dirichlet_boundary([-20/10, 0, 0])         # Outflow
     84Bo = Dirichlet_boundary([-20/10, 0, 0])      # Outflow
    9485Br = Reflective_boundary(domain)             # Solid reflective wall
    9586
    96 #domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br})
    97 domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
     87domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br})
    9888
    9989
     
    10191# Evolve system through time
    10292#------------------------------------------------------------------------------
    103 for t in domain.evolve(yieldstep = 0.2, finaltime = 10.0):
    104     domain.write_time()
     93for t in domain.evolve(yieldstep = 0.2, finaltime = 20):
     94    domain.write_time(track_speeds=True)
    10595
Note: See TracChangeset for help on using the changeset viewer.