Changeset 6298 for anuga_work/production/australia_ph2/sydney/project.py

Timestamp:

Feb 9, 2009, 4:26:06 PM (16 years ago)

Author:

jgriffin

Message:

scripts modified to include file_length function

File:

• anuga_work/production/australia_ph2/sydney/project.py (modified) (7 diffs)

anuga_work/production/australia_ph2/sydney/project.py

@@ -41 +41 @@
 # run_scenario script is executed
 tide = 0                #0.6
-#event_number = 27255   # Java 9.3 worst case for Perth
-#event_number = 68693   # Sumatra 9.2
-event_number = 27283    # Java 9.3 original
+event_number = 7875    # Java 9.3 original
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -95 +93 @@
 # Used in run_busselton,py
 # Format for points easting,northing (no header)
-##land_initial_conditions_filename = [['initial_condition_extend.csv', 0],
-##                                    ['initial_condition_comerong_island.csv', 0],
-##                                    ['initial_condition_gabo_island.csv', 0],
-##                                    ['initial_condition_montague_island.csv', 0]]
+land_initial_conditions_filename = [['initial_condition_extend.csv', 0],
+                                    ['initial_condition_comerong_island.csv', 0],
+                                    ['initial_condition_gabo_island.csv', 0],
+                                    ['initial_condition_montague_island.csv', 0]]
 
 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
@@ -105 +103 @@
 bounding_polygon_filename = 'bounding_polygon.csv'
 
+# INTERIOR REGIONS -  for designing the mesh
+# Used in run_model.py
+# Format for points easting,northing (no header)                    
+##interior_regions_data = [['coast_3km_buffer.csv', 50000]]
+                         
 
 # BOUNDING POLYGON
@@ -119 +122 @@
 # Format easting,northing (no header)
 landward_boundary_filename = 'landward_boundary_extend.csv'
+
+
 
 
@@ -148 +153 @@
 
 # The absolute pathname of the mesh, generated in run_busselton.py
-meshes = join(meshes_folder, scenario_name, '.msh')
+meshes = meshes_folder + sep + scenario_name + '.msh'
 
 # The absolute pathname for the urs order points, used within build_boundary.py
@@ -158 +163 @@
 
 # The absolute pathname for the .sts file, generated in build_boundary.py
-event_sts = join(boundaries_folder, str(event_number), scenario_name)
+event_sts = join(boundaries_folder, str(event_number))
 
 # The absolute pathname for the output folder names
@@ -174 +179 @@
 
 # Create list of land polygons with initial conditions
-##land_initial_conditions = []
-##for filename, MSL in land_initial_conditions_filename:
+land_initial_conditions = []
+for filename, MSL in land_initial_conditions_filename:
+    polygon = read_polygon(join(polygons_folder, filename))
+    land_initial_conditions.append([filename, MSL])
+
+### Create list of interior polygons with scaling factor
+interior_regions = []
+##for filename, maxarea in interior_regions_data:
 ##    polygon = read_polygon(join(polygons_folder, filename))
-##    land_initial_conditions.append([filename, MSL])
-
-# Create list of interior polygons with scaling factor
-interior_regions = []
-#for filename, maxarea in interior_regions_data:
-#    polygon = read_polygon(join(polygons_folder, filename))
-#    interior_regions.append([polygon, maxarea*scale_factor])
+##    interior_regions.append([polygon, maxarea*scale_factor])
 
 # Initial bounding polygon for data clipping 
 bounding_polygon = read_polygon(join(polygons_folder,
                                      bounding_polygon_filename))
-bounding_maxarea = 100000*scale_factor
+bounding_maxarea = 125000*scale_factor
 
 # Estimate the number of triangles