Changeset 6729


Ignore:
Timestamp:
Apr 6, 2009, 1:39:35 PM (15 years ago)
Author:
ole
Message:

Cosmetics in validation test of flows

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/automated_validation_tests/flow_tests/test_inflow_using_flowline.py

    r6705 r6729  
    1818from anuga.shallow_water.shallow_water_domain import Reflective_boundary
    1919from anuga.shallow_water.shallow_water_domain import Dirichlet_boundary
    20 from anuga.shallow_water.shallow_water_domain import Transmissive_Momentum_Set_Stage_boundary
     20from anuga.shallow_water.shallow_water_domain import Transmissive_momentum_set_stage_boundary
    2121from anuga.shallow_water.shallow_water_domain import Inflow
    2222from anuga.shallow_water.data_manager import get_flow_through_cross_section
     
    4646
    4747        domain = Domain(points, vertices, boundary)   
    48         domain.set_name('inflow_flowline_test')              # Output name
     48        domain.set_name('inflow_flowline_test')
    4949
    5050
     
    5757            return z
    5858
    59         domain.set_quantity('elevation', topography)  # Use function for elevation
    60         domain.set_quantity('friction', mannings_n)   # Constant friction of conc surface   
     59        domain.set_quantity('elevation', topography)
     60        domain.set_quantity('friction', mannings_n)
    6161        domain.set_quantity('stage',
    62                             expression='elevation')   # Dry initial condition
     62                            expression='elevation') # Dry initial condition
    6363
    6464
     
    9393        #----------------------------------------------------------------------
    9494
    95         Br = Reflective_boundary(domain)              # Solid reflective wall
    96         #Bo = Dirichlet_boundary([-(slope*length)+normal_depth, 0, 0])  # Outflow stage at normal depth
     95        Br = Reflective_boundary(domain) # Solid reflective wall
     96
     97        # Define downstream boundary based on predicted depth
    9798        def normal_depth_stage_downstream(t):
    98             return (-slope*length) + normal_depth
    99         Bt = Transmissive_Momentum_Set_Stage_boundary(domain=domain,function=normal_depth_stage_downstream)
     99            return (-slope*length) + normal_depth
     100       
     101        Bt = Transmissive_momentum_set_stage_boundary(domain=domain,
     102                                                      function=normal_depth_stage_downstream)
     103       
    100104        domain.set_boundary({'left': Br, 'right': Bt, 'top': Br, 'bottom': Br})
    101105
     
    110114            if verbose :
    111115                print domain.timestepping_statistics()
     116                print domain.volumetric_balance_statistics()               
    112117
    113118
Note: See TracChangeset for help on using the changeset viewer.