Changeset 4141

Timestamp:

Jan 8, 2007, 10:49:08 AM (17 years ago)

Author:

nick

Message:

update perth_2006

File:

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

anuga_work/production/perth_2006/run_perth.py

@@ -44 +44 @@
 #------------------------------------------------------------------------------
 
+start_screen_catcher(project.output_run_time_dir, myid, numprocs)
 
 # filenames
-
 boundaries_name = project.scenario
 meshes_dir_name = project.meshes_dir_name+'.msh'
@@ -59 +59 @@
                  dirname(project.__file__)+sep+ project.__name__+'.py' )
 barrier()
-
-start_screen_catcher(project.output_run_time_dir, myid, numprocs)
 
 print 'USER: ', project.user
@@ -76 +74 @@
 
     print 'start create mesh from regions'
-    interior_regions = [[project.poly_pos20_neg20, project.res_pos20_neg20],
-                    [project.poly_cbd, project.res_cbd],
-                    [project.poly_penguin, project.res_penguin]]    
+#    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]]       
 
@@ -90 +88 @@
 
     print 'start create mesh from regions'
-    create_mesh_from_regions(project.bounding_polygon,
+    create_mesh_from_regions(project.poly_all,
                          boundary_tags={'back': [4], 'side': [0,3],
                                         'ocean': [1, 2]},
-                         maximum_triangle_area=project.res_bounding,
-                         interior_regions=interior_regions,
+                         maximum_triangle_area=project.res_poly_all,
+                         interior_regions=project.interior_regions,
 #                         filename=meshes_time_dir_name,
                          filename=meshes_dir_name,
@@ -151 +149 @@
 
     domain.set_quantity('stage', tide)
-    domain.set_quantity('friction', 0.017) 
+    domain.set_quantity('friction', 0.01) 
     #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name
     print 'Start Set quantity'
 
     domain.set_quantity('elevation', 
-#                    filename = project.combined_dir_name+'.xya',
-                    filename = project.combined_smaller_name_dir+'.xya',
+                    filename = project.combined_dir_name+'.txt',
+#                    filename = project.combined_smaller_name_dir+'.xya',
                     use_cache = True,
                     verbose = True,