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

Timestamp:

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

Author:

ole

Message:

John reducing domain extent and coarsening of some regions.

File:

• anuga_work/development/boxingday08/run_boxingday_polyline.py (modified) (9 diffs)

anuga_work/development/boxingday08/run_boxingday_polyline.py

@@ -40 +40 @@
 # Application specific imports
 import project                 # Definition of file names and polygons
-import create_boundary
 
 
@@ -54 +53 @@
 
 timestamp = strftime('%Y%m%d_%H%M%S',localtime())
-basename = scenario + '_polyline'
+basename = scenario[:4] + '_polyline'
 
 
@@ -71 +70 @@
 # 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
 
 
@@ -82 +81 @@
 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
-boundary_tags={'ocean': arange(1,85).tolist(), 'otherocean': [0,89], '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
 
@@ -105 +104 @@
 
 _ = cache(create_mesh_from_regions,
-          create_boundary.bounding_polygon,
+          project.bounding_polygon,
           {'boundary_tags': boundary_tags,
            'maximum_triangle_area': extent_res,
@@ -111 +110 @@
            'interior_regions': interior_regions},
           verbose = True,
-          dependencies = ['create_boundary.py']
+          dependencies = ['project.py']
           #, evaluate=True
           )
@@ -130 +129 @@
 print domain.statistics()
 
-domain.set_name(basename+'friction'+timestamp)
+domain.set_name(basename+timestamp)
 #domain.set_name(basename)
 domain.set_datadir(scenario)
@@ -168 +167 @@
 #------------------------------------------------------------------------------
 boundary_urs_in='data/boxing'
-boundary_urs_out=create_boundary.base_name
+boundary_urs_out=project.base_name
 
 print 'Available boundary tags', domain.get_boundary_tags()
 Bf = File_boundary(boundary_urs_out+'.sts', 
                    domain, time_thinning=1,
-                   use_cache=False,verbose = True,boundary_polygon=create_boundary.bounding_polygon)
+                   use_cache=True,verbose = True,boundary_polygon=project.bounding_polygon)
 
 Br = Reflective_boundary(domain)
@@ -196 +195 @@
 # 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):