Changeset 6502 for anuga_work/production/australia_ph2/LHI/run_model.py

Timestamp:

Mar 13, 2009, 9:02:45 AM (14 years ago)

Author:

myall

Message:

darwin now working; LHI still has problems with boundary

File:

• anuga_work/production/australia_ph2/LHI/run_model.py (modified) (7 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
-
-# NO LANDWARD BOUNDARY FOR ISLAND MODEL
 ##landward_boundary = read_polygon(project.landward_boundary)
 
@@ -88 +86 @@
 
 # Number of boundary segments
-num_ocean_segments = len(event_sts) #- 1
+num_ocean_segments = len(event_sts) - 1
 # Number of landward_boundary points
 num_land_points = 0 #file_length(project.landward_boundary)
@@ -94 +92 @@
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={#'back': range(num_ocean_segments+1,
-               #              num_ocean_segments+num_land_points),
-               #'side': [num_ocean_segments,
-               #         num_ocean_segments+num_land_points],
+boundary_tags={'back': range(num_ocean_segments+1,
+                             num_ocean_segments+num_land_points),
+               'side': [num_ocean_segments,
+                        num_ocean_segments+num_land_points],
                'ocean': range(num_ocean_segments)}
 
@@ -106 +104 @@
                                     interior_regions=project.interior_regions,
                                     mesh_filename=project.meshes,
-                                    use_cache=False,
+                                    use_cache=True,
                                     verbose=True)
 print domain.statistics()
@@ -127 +125 @@
 else:
     IC = 0
-domain.set_quantity('stage', IC, use_cache=False, verbose=True)
+domain.set_quantity('stage', IC, use_cache=True, verbose=True)
 domain.set_quantity('friction', project.friction) 
 domain.set_quantity('elevation', 
                     filename=project.combined_elevation+'.pts',
-                    use_cache=False,
+                    use_cache=True,
                     verbose=True,
                     alpha=project.alpha)
@@ -155 +153 @@
                         default_boundary=Bd,
                         boundary_polygon=bounding_polygon_sts,                    
-                        use_cache=False,
+                        use_cache=True,
                         verbose=True)
-    domain.set_boundary({'ocean': Bf}) 
+    domain.set_boundary({'back': Br,
+                     'side': Bd,
+                     'ocean': Bf}) 
     
 elif project.wave == 'Tb':
@@ -163 +163 @@
                        verbose=True)
 
-    domain.set_boundary({'ocean': Tb})
+    domain.set_boundary({'back': Br,
+                         'side': Bd,
+                         'ocean': Tb})
 else:
     print 'No wave specified in project script (Bf or Tb)'