Changeset 8248 for trunk/anuga_work

Timestamp:

Nov 15, 2011, 4:36:00 PM (13 years ago)

Author:

steve

Message:

commiting at the anuga_core level

Location:

trunk/anuga_work/shallow_water_balanced_steve

Files:

• run_step.py (modified) (2 diffs)
• run_wave.py (modified) (3 diffs)
• swb_domain.py (modified) (1 diff)

trunk/anuga_work/shallow_water_balanced_steve/run_step.py

@@ -12 +12 @@
 import sys
 import anuga
+#from anuga import Domain
 from swb_domain import *
 
@@ -46 +47 @@
 # structured mesh
 points, vertices, boundary = anuga.rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
+
+
+create_mesh_from_regions(bounding_polygon,
+                             boundary_tags,
+                             maximum_triangle_area=maximum_triangle_area,
+                             interior_regions=interior_regions,
+                             filename=mesh_filename,
+                             interior_holes=interior_holes,
+                             hole_tags=hole_tags,
+                             poly_geo_reference=poly_geo_reference,
+                             mesh_geo_reference=mesh_geo_reference,
+                             minimum_triangle_angle=minimum_triangle_angle,
+                             fail_if_polygons_outside=fail_if_polygons_outside,
+                             use_cache=False,
+                             verbose=verbose)
 
 domain = Domain(points, vertices, boundary) 

trunk/anuga_work/shallow_water_balanced_steve/run_wave.py

@@ -12 +12 @@
 import sys
 import anuga
-from swb_domain import *
+from anuga import Domain
+#from swb_domain import *
 
 from math import cos
@@ -56 +57 @@
 domain.set_timestepping_method('rk2')
 domain.set_default_order(2)
+domain.set_beta(2.0)
 
 print domain.get_timestepping_method()
@@ -79 +81 @@
 Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
 amplitude = 1
-wave_length = 1000.0
+wave_length = 300.0
 Bw = anuga.Time_boundary(domain=domain,     # Time dependent boundary
 ## Sine wave

trunk/anuga_work/shallow_water_balanced_steve/swb_domain.py

@@ -77 +77 @@
         self.forcing_terms[1] = swb_gravity
 
+        print "Swb_domain"
 
     def check_integrity(self):