Ignore:
Timestamp:
Sep 1, 2010, 6:11:25 PM (12 years ago)
Author:
steve
Message:

Added enquiry points to culverts

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py

    r7980 r7984  
    4646
    4747points, vertices, boundary = rectangular_cross(int(length/dx), int(width/dy),
    48                                                len1=length, len2=width)
     48                                                    len1=length, len2=width)
    4949domain = Domain(points, vertices, boundary)   
    5050domain.set_name('Test_WIDE_BRIDGE')                 # Output name
     
    168168"""
    169169from anuga.structures.culvert_operator import Culvert_operator
     170#culvert0 = Culvert_operator(domain,
     171#                            end_point0=[40.0, 75.0],
     172#                            end_point1=[50.0, 75.0],
     173#                            width=50.0,
     174#                            height=10.0,
     175#                            apron=5.0,
     176#                            verbose=False)
     177
     178
    170179culvert1 = Culvert_operator(domain,
    171                             end_point0=[40.0, 75.0],
    172                             end_point1=[50.0, 75.0],
    173                             width=50.0,
     180                            end_point0=[40.0, 87.5],
     181                            end_point1=[50.0, 87.5],
     182                            width=25.0,
    174183                            height=10.0,
     184                            apron=5.0,
     185                            verbose=False)
     186
     187
     188culvert2 = Culvert_operator(domain,
     189                            end_point0=[40.0, 62.5],
     190                            end_point1=[50.0, 62.5],
     191                            width=25.0,
     192                            height=10.0,
     193                            apron=5.0,
    175194                            verbose=False)
    176195
     
    220239for t in domain.evolve(yieldstep = 1, finaltime = 100):
    221240    print domain.timestepping_statistics()
     241    print domain.volumetric_balance_statistics()
    222242   
    223243
Note: See TracChangeset for help on using the changeset viewer.