Changeset 5401 for anuga_work/development/boxingday08/run_boxingday_urs.py

Timestamp:

Jun 13, 2008, 10:29:14 AM (15 years ago)

Author:

ole

Message:

John reducing domain extent and coarsening of some regions.

File:

• anuga_work/development/boxingday08/run_boxingday_urs.py (modified) (8 diffs)

anuga_work/development/boxingday08/run_boxingday_urs.py

@@ -40 +40 @@
 # Application specific imports
 import project                 # Definition of file names and polygons
-import create_boundary
 
 #------------------------------------------------------------------------------
@@ -53 +52 @@
 
 timestamp = strftime('%Y%m%d_%H%M%S',localtime())
-basename = scenario + '_urs'
+basename = scenario[:4] + '_urs'
 
 
@@ -70 +69 @@
 # resolutions (maximal area of per triangle) for each polygon
 #------------------------------------------------------------------------------
-extent_res     = 10000000.0
+extent_res     = 1000000.0
 contour20m_res = 50000.0
 island_res     = 5000.0
 bay_res        = 2000.0
-patong_res     = 100.0
+patong_res     = 400.0
 
 # make polygon that contains land that does not affect result.
@@ -80 +79 @@
 interior_regions = [[project.patong, patong_res],
                     [project.bay, bay_res],
-                    [project.contour20m, contour20m_res],
-                    [project.island_north, island_res],
-                    [project.island_south, island_res],
-                    [project.island_south2, island_res]]
+                    [project.contour20m, contour20m_res]]#,
+                    #[project.island_north, island_res],
+                    #[project.island_south, island_res],
+                    #[project.island_south2, island_res]]
 
 #for coarse run to test gauges
-extent_res = 10000000.0
-interior_regions=None
+#extent_res = 10000000.0
+#interior_regions=None
 
 from Numeric import arange,allclose
-print 'Beware not all boundary segments that are tagged as ocean in polyline are tagged here. missing are 79,80,81,82,83,84 which are set to other ocean'
-boundary_tags={'ocean': arange(1,79).tolist(), 'otherocean': [0,89,84,83,82,81,80,79], 'land': [86, 87], 'both': [85,88]}
-
-
-#trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)
+boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]}
+
+
+#trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res)
 #print 'Minimum number of traingles ', trigs_min
 
@@ -104 +102 @@
 
 _ = cache(create_mesh_from_regions,
-          create_boundary.bounding_polygon,
+          project.bounding_polygon,
           {'boundary_tags': boundary_tags,
            'maximum_triangle_area': extent_res,
@@ -110 +108 @@
            'interior_regions': interior_regions},
           verbose = True,
-          #dependencies = ['create_boundary.py']
+          #dependencies = ['project.py']
           dependencies = ['run_boxingday_urs.py']
           #, evaluate=True
@@ -131 +129 @@
 print domain.statistics()
 
-domain.set_name(basename+'friction'+timestamp)
+domain.set_name(basename+timestamp)
 #domain.set_name(basename)
 domain.set_datadir(scenario)
@@ -208 +206 @@
 # start when boundary stage first becomes non-zero.
 
-for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0, 
+for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0, 
                        skip_initial_step = False):
     if allclose(t,10800.0):