Changeset 6911

Timestamp:

Apr 28, 2009, 11:10:30 AM (16 years ago)

Author:

jgriffin

Message:

 

Location:

anuga_work/production/australia_ph2/bamaga

Files:

• project.py (modified) (5 diffs)
• run_model.py (modified) (1 diff)

anuga_work/production/australia_ph2/bamaga/project.py

@@ -17 +17 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'australia_ph2'
-scenario_name = 'gulf'
+scenario_name = 'weipa'
 
 #-------------------------------------------------------------------------------
@@ -33 +33 @@
 ##else:
 
-event_number = 64977    # the event number or the mux file name
+event_number = 51212    # the event number or the mux file name
 
 #event_number_list = [31975, 51390, 63735] # To piggy back multiple events
@@ -41 +41 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=80000 # final time for simulation
-setup = 'final'         # This can be one of three values
+finaltime=1000 # final time for simulation
+setup = 'trial'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
@@ -48 +48 @@
 
 # index is only used when wave = Tb
-index = 1914            # index from the PTHA - Y2000 0.257m
+index = 2003            # index from the PTHA - Y2000 0.257m
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
@@ -79 +79 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['gulf_bg'] # 250m grid 2005
+ascii_grid_filenames = ['weipa_bg_pr'] # 250m grid 2005
 
 # Format for points is x,y,elevation (with header)

anuga_work/production/australia_ph2/bamaga/run_model.py

@@ -87 +87 @@
 # Number of boundary segments
 num_ocean_segments = len(event_sts) - 1
+print num_ocean_segments
 # Number of landward_boundary points
 num_land_points = file_length(project.landward_boundary)
-
+print num_land_points
 # 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],
-##               'ocean': range(num_ocean_segments)}
-boundary_tags = {'back':[3],
-                 'side':[2,4],
-                 'ocean':[0,1,5]}
+boundary_tags={'back': range(num_ocean_segments+2,
+                             num_ocean_segments+5),
+               'side': [num_ocean_segments,
+                        num_ocean_segments+1,
+                        num_ocean_segments+5,
+                        num_ocean_segments+6,
+                        num_ocean_segments+7,
+                        num_ocean_segments+8,
+                        num_ocean_segments+9],
+               'ocean': range(num_ocean_segments)}
+##boundary_tags = {'back':[3],
+##                 'side':[2,4],
+##                 'ocean':[0,1,5]}
 
 print 'boundary tags', boundary_tags