Changeset 7846 for anuga_validation/convergence_study/dam_break.py
- Timestamp:
- Jun 16, 2010, 8:23:42 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_validation/convergence_study/dam_break.py
r5442 r7846 10 10 11 11 import sys 12 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross 13 from anuga.shallow_water import Domain 14 from anuga.shallow_water import Reflective_boundary 15 from anuga.shallow_water import Dirichlet_boundary 16 from anuga.shallow_water import Time_boundary 17 from anuga.shallow_water import Transmissive_boundary 18 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary 19 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files 12 import anuga 20 13 21 14 from math import cos 22 from Numeric import zeros, Float15 from numpy import zeros, float 23 16 from time import localtime, strftime, gmtime 24 17 … … 34 27 output_file = 'dam_break' 35 28 36 copy_code_files(output_dir,__file__)29 anuga.copy_code_files(output_dir,__file__) 37 30 #start_screen_catcher(output_dir+'_') 38 31 … … 47 40 48 41 # structured mesh 49 points, vertices, boundary = rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))42 points, vertices, boundary = anuga.rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2)) 50 43 51 domain = Domain(points, vertices, boundary)44 domain = anuga.Domain(points, vertices, boundary) 52 45 53 46 domain.set_name(output_file) … … 97 90 #------------------------------------------------------------------------------ 98 91 from math import sin, pi, exp 99 Br = Reflective_boundary(domain) # Solid reflective wall100 Bt = Transmissive_boundary(domain) # Continue all values on boundary101 Bd = Dirichlet_boundary([1,0.,0.]) # Constant boundary values92 Br = anuga.Reflective_boundary(domain) # Solid reflective wall 93 Bt = anuga.Transmissive_boundary(domain) # Continue all values on boundary 94 Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values 102 95 103 96
Note: See TracChangeset
for help on using the changeset viewer.