Changeset 6511 for anuga_work/production/australia_ph2/LHI

Timestamp:

Mar 13, 2009, 1:25:37 PM (16 years ago)

Author:

jgriffin

Message:

Fized default boundary in Field_boundary

File:

• anuga_work/production/australia_ph2/LHI/run_model.py (modified) (3 diffs)

anuga_work/production/australia_ph2/LHI/run_model.py

@@ -80 +80 @@
 # Reading the landward defined points, this incorporates the original clipping
 # polygon minus the 100m contour
-##landward_boundary = read_polygon(project.landward_boundary)
+landward_boundary = read_polygon(project.landward_boundary)
 
 # Combine sts polyline with landward points
-bounding_polygon_sts = event_sts #+ landward_boundary
+bounding_polygon_sts = event_sts + landward_boundary
 
 # Number of boundary segments
 num_ocean_segments = len(event_sts) - 1
 # Number of landward_boundary points
-num_land_points = 0 #file_length(project.landward_boundary)
+num_land_points = file_length(project.landward_boundary)
 
 # Boundary tags refer to project.landward_boundary
@@ -151 +151 @@
                         domain, mean_stage=project.tide,
                         time_thinning=1,
-                        default_boundary=Bd,
+                        default_boundary=Dirichlet_boundary([0, 0, 0]),
                         boundary_polygon=bounding_polygon_sts,                    
                         use_cache=True,
@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd,
-                       verbose=True)
+    Tb = Time_boundary(domain,f,default_boundary=Bd )
 
     domain.set_boundary({'back': Br,