Changeset 6502

Timestamp:

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

Author:

myall

Message:

darwin now working; LHI still has problems with boundary

Location:

anuga_work/production/australia_ph2

Files:

• LHI/run_model.py (modified) (7 diffs)
• LHI/run_model_only_ocean.py (added)
• darwin/project.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
-
-# 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)'

anuga_work/production/australia_ph2/darwin/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 52               # specify zone of model
-event_number = 71261    # the event number or the mux file name
+event_number = 70906    # the event number or the mux file name
 # contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
 # 71261 (Timor_OTF)
@@ -34 +34 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 # index is only used when wave = Tb
@@ -40 +40 @@
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
-setup = 'trial'         # This can be one of three values
+setup = 'final'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh