- Timestamp:
- Jul 22, 2009, 9:22:11 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/show_balanced_limiters.py
r6145 r7317 16 16 Transmissive_boundary, Time_boundary 17 17 from anuga.shallow_water.shallow_water_domain import Weir_simple as Weir 18 import anuga.utilities.log as log 18 19 19 20 from mesh_factory import rectangular … … 26 27 N = 12 27 28 28 print 'Creating domain' 29 log.critical('Creating domain') 29 30 #Create basic mesh 30 31 points, vertices, boundary = rectangular(N, N/2, len1=1.2,len2=0.6, 31 32 origin=(-0.07, 0)) 32 33 33 print 'Number of elements', len(vertices)34 log.critical('Number of elements=%d' % len(vertices)) 34 35 #Create shallow water domain 35 36 domain = Domain(points, vertices, boundary) … … 45 46 Z = Weir(inflow_stage) 46 47 47 print 'Field values' 48 log.critical('Field values') 48 49 domain.set_quantity('elevation', Z) 49 50 domain.set_quantity('friction', manning) … … 53 54 # Boundary conditions 54 55 # 55 print 'Boundaries' 56 log.critical('Boundaries') 56 57 Br = Reflective_boundary(domain) 57 58 Bt = Transmissive_boundary(domain) … … 74 75 #Initial condition 75 76 # 76 print 'Initial condition' 77 log.critical('Initial condition') 77 78 domain.set_quantity('stage', Z) 78 79 … … 82 83 domain.write_boundary_statistics(['stage'],'left') 83 84 84 print 'Done' 85 log.critical('Done') 85 86 86 87
Note: See TracChangeset
for help on using the changeset viewer.