Changeset 8078

Timestamp:

Nov 19, 2010, 6:20:08 PM (13 years ago)

Author:

habili

Message:

updated code so that wave.py runs

File:

• trunk/anuga_validation/validation_tests/periodic_forcing/wave.py (modified) (4 diffs)

trunk/anuga_validation/validation_tests/periodic_forcing/wave.py

@@ -11 +11 @@
 
 import sys
-from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
-from anuga.shallow_water import Domain
-from anuga.shallow_water import Reflective_boundary
-from anuga.shallow_water import Dirichlet_boundary
-from anuga.shallow_water import Time_boundary
-from anuga.shallow_water import Transmissive_boundary
-from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
-from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
-
+import anuga
 from math import cos
-from Numeric import zeros, Float
 from time import localtime, strftime, gmtime
 from os import sep
@@ -36 +27 @@
 output_file = 'wave'
 
-copy_code_files(output_dir,__file__)
+#copy_code_files(output_dir,__file__)
 #start_screen_catcher(output_dir+sep)
 
@@ -50 +41 @@
 
 # structured mesh
-points, vertices, boundary = rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
+points, vertices, boundary = anuga.rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
 
-domain = Domain(points, vertices, boundary) 
+domain = anuga.Domain(points, vertices, boundary) 
 
 domain.set_name(output_file)                
@@ -90 +81 @@
 #------------------------------------------------------------------------------
 from math import sin, pi, exp
-Br = Reflective_boundary(domain)      # Solid reflective wall
-Bt = Transmissive_boundary(domain)    # Continue all values on boundary 
-Bd = Dirichlet_boundary([1,0.,0.]) # Constant boundary values
+Br = anuga.Reflective_boundary(domain)      # Solid reflective wall
+Bt = anuga.Transmissive_boundary(domain)    # Continue all values on boundary 
+Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
 amplitude = 1
-Bw = Time_boundary(domain=domain,     # Time dependent boundary  
+Bw = anuga.Time_boundary(domain=domain,     # Time dependent boundary  
 ## Sine wave
                   f=lambda t: [(-amplitude*sin((1./300.)*t*2*pi)), 0.0, 0.0])