Changeset 4091 for anuga_work/production/perth_2006/run_perth.py

Timestamp:

Dec 19, 2006, 9:19:27 AM (18 years ago)

Author:

nick

Message:

updates to dampier and perth

File:

• anuga_work/production/perth_2006/run_perth.py (modified) (8 diffs)

anuga_work/production/perth_2006/run_perth.py

@@ -31 +31 @@
 
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
-
 from anuga.geospatial_data.geospatial_data import *
 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
 from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier
+from anuga.utilities.polygon import plot_polygons, polygon_area
 
 # Application specific imports
@@ -47 +47 @@
 # filenames
 
-build_time = '20061107_063840_build'
-bound_time = '20061102_221245_build'
-
 boundaries_name = project.scenario
 meshes_dir_name = project.meshes_dir_name+'.msh'
-#source_dir = project.boundarydir
 #boundaries_time_dir_name = project.boundaries_time_dir_name
 boundaries_dir_name = project.boundaries_dir_name
@@ -67 +63 @@
 
 print 'USER: ', project.user
+print 'min triangles', project.trigs_min,
+print 'Note: This is generally about 20% less than the final amount'
 
 #--------------------------------------------------------------------------
@@ -76 +74 @@
 
 if myid == 0:
+
     print 'start create mesh from regions'
-    interior_regions = [[project.poly_pos20_neg20, 500000],
-                    [project.poly_cbd, 50000],
-                    [project.poly_penguin, 50000]]    
+    interior_regions = [[project.poly_pos20_neg20, project.res_pos20_neg20],
+                    [project.poly_cbd, project.res_cbd],
+                    [project.poly_penguin, project.res_penguin]]    
 #                    [project.poly_interior, 1000]]       
 
@@ -89 +88 @@
 #               verbose = True)
 #    import sys; sys.exit()
+
     print 'start create mesh from regions'
-
     create_mesh_from_regions(project.bounding_polygon,
                          boundary_tags={'back': [4], 'side': [0,3],
                                         'ocean': [1, 2]},
-                         maximum_triangle_area=1000000,
+                         maximum_triangle_area=project.res_bounding,
                          interior_regions=interior_regions,
 #                         filename=meshes_time_dir_name,
@@ -157 +156 @@
 
     domain.set_quantity('elevation', 
-                    filename = project.combined_dir_name+'.xya',
+#                    filename = project.combined_dir_name+'.xya',
+                    filename = project.combined_smaller_name_dir+'.xya',
                     use_cache = True,
                     verbose = True,
@@ -201 +201 @@
 domain.set_boundary({'back': Br,
                      'side': Bd,
-                     'ocean': Bf}) 
+                     'ocean': Bd}) 
 print'finish set boundary'
 
@@ -214 +214 @@
     domain.write_boundary_statistics(tags = 'ocean')
     if allclose(t, 120):
-        domain.set_boundary({'back': Br, 'side': Bo, 'ocean': Bo})
+        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
 
     if allclose(t, 1020):
-        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
+        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})