Changeset 8067 for trunk/anuga_validation/convergence_study/wave.py
- Timestamp:
- Nov 4, 2010, 5:23:35 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_validation/convergence_study/wave.py
r5442 r8067 11 11 12 12 import sys 13 from anuga .abstract_2d_finite_volumes.mesh_factoryimport rectangular_cross14 from anuga .shallow_waterimport Domain15 from anuga .shallow_waterimport Reflective_boundary16 from anuga .shallow_waterimport Dirichlet_boundary17 from anuga .shallow_waterimport Time_boundary18 from anuga .shallow_waterimport Transmissive_boundary19 from anuga .shallow_water import Transmissive_Momentum_Set_Stage_boundary20 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files13 from anuga import rectangular_cross 14 from anuga import Domain 15 from anuga import Reflective_boundary 16 from anuga import Dirichlet_boundary 17 from anuga import Time_boundary 18 from anuga import Transmissive_boundary 19 from anuga import Transmissive_momentum_set_stage_boundary 20 #from anuga.data_manager import start_screen_catcher, copy_code_files 21 21 22 22 from math import cos 23 from Numeric import zeros, Float 23 import numpy as num 24 24 from time import localtime, strftime, gmtime 25 25 from os import sep … … 36 36 output_file = 'wave' 37 37 38 copy_code_files(output_dir,__file__)38 #copy_code_files(output_dir,__file__) 39 39 #start_screen_catcher(output_dir+sep) 40 40 41 interactive_visualisation = False41 interactive_visualisation = True 42 42 43 43 #------------------------------------------------------------------------------ … … 60 60 # Setup Algorithm 61 61 #------------------------------------------------------------------------------ 62 domain.set_timestepping_method('rk2') 62 63 domain.set_timestepping_method(2) 63 64 domain.set_default_order(2) 64 65 65 66 print domain.get_timestepping_method() 66 67 67 domain.use_edge_limiter = True 68 domain.tight_slope_limiters = False 69 domain.use_centroid_velocities = False 68 #domain.use_edge_limiter = True 69 #domain.tight_slope_limiters = False 70 #domain.use_centroid_velocities = False 71 72 domain.set_beta(1.9) 70 73 71 74 domain.CFL = 1.0 72 75 73 domain.beta_w = 1.074 domain.beta_w_dry = 0.075 domain.beta_uh = 1.076 domain.beta_uh_dry = 0.077 domain.beta_vh = 1.078 domain.beta_vh_dry = 0.076 #domain.beta_w = 1.0 77 #domain.beta_w_dry = 0.0 78 #domain.beta_uh = 1.0 79 #domain.beta_uh_dry = 0.0 80 #domain.beta_vh = 1.0 81 #domain.beta_vh_dry = 0.0 79 82 80 83
Note: See TracChangeset
for help on using the changeset viewer.