Changeset 4372

Timestamp:

Apr 13, 2007, 4:59:58 PM (18 years ago)

Author:

nick

Message:

update broome

Location:

anuga_work/production/broome_2006

Files:

• project.py (modified) (5 diffs)
• run_broome.py (modified) (4 diffs)

anuga_work/production/broome_2006/project.py

@@ -28 +28 @@
 print 'gtime: ', gtime
 
-tide = 0
+#tide = -5.3
+#tide = 0
+tide = 4.9
 
 #Making assumptions about the location of scenario data
@@ -51 +53 @@
 combined_name ='broome_combined_elevation'
 combined_name1 ='broome_combined_elevation1'
+combined_name_unclipped1 ='broome_combined_elevation_unclipped1'
 combined_small_name = 'broome_combined_elevation_small'
 
@@ -79 +82 @@
 #final topo files
 combined_dir_name = topographies_dir + combined_name
+combined_dir_name_unclipped1 = topographies_dir + combined_name_unclipped1
 combined_dir_name1 = topographies_dir + combined_name1
 combined_small_name_dir = topographies_dir + combined_small_name
@@ -120 +124 @@
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 
-poly_all = read_polygon(polygons_dir+'extent.csv')
+poly_all = read_polygon(polygons_dir+'extent_small.csv')
+#poly_all = read_polygon(polygons_dir+'extent.csv')
 res_poly_all = 500000
+#res_poly_all = 150000
 
 ###############################
@@ -129 +135 @@
 poly_0 = read_polygon(polygons_dir+'neg20_coast_contour_pts.csv')
 res_0 = 100000
+#res_0 = 20000
 
 poly_1 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
 res_1 = 50000
+#res_1 = 5000
 
 poly_2 = read_polygon(polygons_dir+'broome_south_coast_inside_extent.csv')
 res_2 = 50000
+#res_2 = 5000
 
 poly_3 = read_polygon(polygons_dir+'Broome_town_pts.csv')
 res_3 = 20000
+#res_3 = 2000
 
 poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
 res_4 = 5000
+#res_4 = 500
 #assert zone == refzone
 

anuga_work/production/broome_2006/run_broome.py

@@ -32 +32 @@
 from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier
 from anuga.utilities.polygon import plot_polygons, polygon_area
+from anuga_parallel.parallel_abstraction import get_processor_name
 
 # Application specific imports
@@ -56 +57 @@
 
 # creates copy of code in output dir
-
+print "Processor Name:",get_processor_name()
 
 print 'USER: ', project.user
@@ -73 +74 @@
     print 'start create mesh from regions'
     create_mesh_from_regions(project.poly_all,
+#                         boundary_tags={'back': [1, 2], 'side': [0,3],
+#                                        'ocean': [4, 5, 6]},
                          boundary_tags={'back': [1, 2], 'side': [0,3],
-                                        'ocean': [4, 5, 6]},
+                                        'ocean': [4]},
                          maximum_triangle_area=project.res_poly_all,
                          interior_regions=project.interior_regions,
@@ -135 +138 @@
     #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name
     print 'Start Set quantity'
-
+    print 'project.combined_dir_name_unclipped1',project.combined_dir_name_unclipped1+'.txt'
     domain.set_quantity('elevation', 
-                    filename = project.combined_dir_name1+'.txt',
+                    filename = project.combined_dir_name_unclipped1+'.txt',
 #                    filename = project.combined_smaller_name_dir+'.xya',
                     use_cache = True,