- Timestamp:
- Apr 5, 2012, 12:37:54 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_work/development/gareth/tests/wave/run_wave.py
r8354 r8384 13 13 import anuga 14 14 from anuga import Domain 15 #from balanced_dev import * 15 16 from balanced_dev import * 17 from balanced_dev import Domain as Domain 16 18 17 19 from math import cos … … 39 41 # Setup domain 40 42 #------------------------------------------------------------------------------ 41 dx = 500.43 dx = 1000. 42 44 dy = dx 43 45 L = 100000. … … 55 57 # Setup Algorithm 56 58 #------------------------------------------------------------------------------ 57 domain.set_timestepping_method('rk2')58 domain.set_default_order(2)59 domain.set_beta(2.0)59 #domain.set_timestepping_method('rk2') 60 #domain.set_default_order(2) 61 #domain.set_beta(2.0) 60 62 61 print domain.get_timestepping_method()63 #print domain.get_timestepping_method() 62 64 63 65 #domain.use_edge_limiter = True … … 81 83 Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values 82 84 amplitude = 1 83 wave_length = 300 0.085 wave_length = 300.0 84 86 Bw = anuga.Time_boundary(domain=domain, # Time dependent boundary 85 87 ## Sine wave … … 116 118 #------------------------------------------------------------------------------ 117 119 118 for t in domain.evolve(yieldstep = 10.0, finaltime = 60*60.* 10.):120 for t in domain.evolve(yieldstep = 10.0, finaltime = 60*60.*3.): 119 121 domain.write_time() 120 122 if interactive_visualisation:
Note: See TracChangeset
for help on using the changeset viewer.