Changeset 4570 for anuga_work/production/shark_bay_2007/project.py

Timestamp:

Jun 28, 2007, 5:54:00 PM (18 years ago)

Author:

ole

Message:

Setup shark bay initial runs

File:

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

anuga_work/production/shark_bay_2007/project.py

@@ -32 +32 @@
 
 #tide = -3.9
-#tide = 0.0
-tide = 3.6
+tide = 0.0
+#tide = 3.6
 
 #Maybe will try to make project a class to allow these parameters to be passed in.
@@ -40 +40 @@
 finaltime=25000
 starttime=3600
-setup='final'
+setup='basic'
 source='shark_bay'
 
@@ -67 +67 @@
                    'clipped_bathymetry_final', 'coast_points_final']
 
-#final topo name
+# final topo name
 combined_name ='shark_bay_combined_elevation'
 combined_small_name = 'shark_bay_combined_elevation_small'
@@ -74 +74 @@
 anuga_dir = join(home, state, scenario, 'anuga')
 
-topographies_in_dir = join(home, state, scenario, 'elevation_final', 'test_area')
+topographies_in_dir = join(home, state, scenario,
+                           'elevation_final', 'test_area')
+
 topographies_dir = join(anuga_dir, 'topographies') # Output dir for ANUGA
 
@@ -87 +89 @@
 
 
-meshes_dir = join(anuga_dir, 'meshes')
-meshes_dir_name = join(meshes_dir, scenario_name)
+mesh_dir = join(anuga_dir, 'meshes')
+mesh_name = join(mesh_dir, scenario_name)
 
 
@@ -97 +99 @@
 # locations for boundary conditions
 if source=='shark_bay':
-    boundaries_file_name = 'shark_bay_3867_18052007'
-    boundaries_dir = join(anuga_dir, 'boundaries', 'shark_bay', '1_10000')
+    boundary_file_name = 'shark_bay_3867_18052007'
+    boundary_dir = join(anuga_dir, 'boundaries', 'shark_bay', '1_10000')
 
-boundaries_name = join(boundaries_dir, boundaries_file_name)
+boundary_name = join(boundary_dir, boundary_file_name)
 
 #output locations
-output_dir = join(anuga_dir, 'outputs')
+output_dir = join(anuga_dir, 'outputs')+sep
 output_build_time_dir = output_dir+build_time+sep
 output_run_time_dir = output_dir +run_time+dir_comment+sep
@@ -147 +149 @@
 
 #from anuga.coordinate_transforms.redfearn import redfearn
-poly_all = read_polygon(join(polygons_dir, 'boundary_extent.csv'))
-res_poly_all = 250000*res_factor
+bounding_polygon = read_polygon(join(polygons_dir, 'boundary_extent.csv'))
+res_bounding_polygon = 250000*res_factor
 
 #Interior regions
 poly_inundated_area = read_polygon(join(polygons_dir, 'inundated_area.csv'))
-res_inundated_area = 50000*res_factor
+res_inundated_area = 1000*res_factor
 
 poly_bay_area = read_polygon(join(polygons_dir, 'bay_area.csv'))                                   
-res_bay_area = 100000*res_factor                                   
+res_bay_area = 5000*res_factor                                   
 
-interior_regions = [[poly_bay_area, res_bay_area], [poly_inundated_area, res_inundated_area]]
+interior_regions = [[poly_bay_area, res_bay_area],
+                    [poly_inundated_area, res_inundated_area]]
 
-trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+trigs_min = number_mesh_triangles(interior_regions,
+                                  bounding_polygon,
+                                  res_bounding_polygon)
 
 print 'min number triangles', trigs_min
 
+onshore_polygon = read_polygon(join(topographies_in_dir,
+                                    'initial_condition.txt'))