Changeset 5604 for anuga_work/development/convergence_okushiri_2008

Timestamp:

Aug 5, 2008, 9:50:22 AM (17 years ago)

Author:

Leharne

Message:

Update project_truescale.py to include x&y extents for export_results.py

Location:

anuga_work/development/convergence_okushiri_2008

Files:

• export_results.py (added)
• project_truescale.py (modified) (3 diffs)

anuga_work/development/convergence_okushiri_2008/project_truescale.py

@@ -43 +43 @@
 run_time = time+'_run'
 
+# Set anuga input directory names
+
+anuga_dir = join(home,'anuga')+sep
+
+mesh_dir = join(anuga_dir, 'meshes')+sep
+mesh_name = join(mesh_dir, 'okushiri_truescale')
+
+polygons_dir = join(anuga_dir, 'polygons')+sep # Created with ArcGIS (csv files)
+
+#------------------------
 # Run parameters
+#------------------------
+
 finaltime=450
-setup='original'
-
+setup='sixteen'
+polygons = 'contour_polygons'
+
+if setup =='no polygons':
+    print 'no interior polygons'
+    base_resolution=1
+    yieldstep=1
+if setup =='octuple':
+    print '8 times original resolution'
+    base_resolution=0.125
+    yieldstep=1
+if setup =='sixteen':
+    print '16 times original resolution'
+    base_resolution=0.0625
+    yieldstep=1
+if setup =='quadruple':
+    print '4 times original resolution'
+    base_resolution=0.25
+    yieldstep=1
+if setup =='double':
+    print 'double original resolution'
+    base_resolution=0.5
+    yieldstep=1
+if setup =='1.5':
+    print '1.5 times original resolution'
+    base_resolution=0.66
+    yieldstep=1
 if setup =='original':
     print 'original resolution'
     base_resolution=1
     yieldstep=1
-if setup =='double':
-    print 'double original resolution'
-    base_resolution=0.5
+if setup =='0.75':
+    print '0.75 times original resolution'
+    base_resolution=1.33
     yieldstep=1
 if setup =='half':
@@ -59 +96 @@
     base_resolution=2
     yieldstep=1
-if setup =='no polygons':
-    print 'half original resolution'
-    base_resolution=2
+if setup =='quarter':
+    print '1/4 original resolution'
+    base_resolution=4
+    yieldstep=1
+if setup =='eighth':
+    print '1/8 original resolution'
+    base_resolution=8
+    yieldstep=1
+if setup =='sixteenth':
+    print '1/16 original resolution'
+    base_resolution=16
+    yieldstep = 1
+if setup =='1-32':
+    print '1/32 original resolution'
+    base_resolution=32
+    yieldstep=1
+if setup =='1-64':
+    print '1/64 original resolution'
+    base_resolution=64
     yieldstep=1
   
+  
+#------------------------------
+# Polygon definitions
+#------------------------------
     
-# Set anuga directories
-anuga_dir = join(home,'anuga')+sep
-
-dir_comment='_'+setup+'_'+str(user)
-
-mesh_dir = join(anuga_dir, 'meshes')+sep
-mesh_name = join(mesh_dir, 'okushiri_truescale')
-
-polygons_dir = join(anuga_dir, 'polygons')+sep # Created with ArcGIS (csv files)
-
-# Output locations
+poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
+res_poly_all = 16000*base_resolution
+
+# Original polygon definitions
+
+poly_gulleys = read_polygon(polygons_dir+'gulleys_polygon.csv')
+res_gulleys = 3.2*base_resolution
+
+poly_island = read_polygon(polygons_dir+'island_polygon.csv')
+res_island = 32*base_resolution
+
+poly_rhs = read_polygon(polygons_dir+'rhs_polygon.csv')
+res_rhs = 80*base_resolution
+
+                    
+# Contour-based polygon definitions
+
+poly_25 = read_polygon(polygons_dir+'polygon_25m.csv')
+res_poly_25 = 800*base_resolution
+
+poly_10 = read_polygon(polygons_dir+'polygon_10m.csv')
+res_poly_10 = 80*base_resolution
+
+poly_5 = read_polygon(polygons_dir+'polygon_5m.csv')
+res_poly_5 = 32*base_resolution
+
+poly_1 = read_polygon(polygons_dir+'polygon_1m.csv')
+res_poly_1 = 10*base_resolution
+
+                    
+if polygons =='original_polygons':
+    print 'original polygon definition'
+    interior_regions = [[poly_gulleys,res_gulleys],[poly_island,res_island],
+                    [poly_rhs,res_rhs]]
+
+if polygons =='contour_polygons':
+    print 'contour-based polygon definition'
+    interior_regions = [[poly_25,res_poly_25],[poly_10,res_poly_10],
+                    [poly_5,res_poly_5],[poly_1,res_poly_1] ]
+
+
+trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+print 'min number triangles', trigs_min
+
+this_area = polygon_area(poly_all)
+print this_area
+
+
+# Polygon QC
+#polygon_plot='polygons.png'
+#plot_polygons([poly_all, poly_25, poly_10, poly_5, poly_1], style=None, figname=polygon_plot)
+
+
+#--------------------------------
+# Set anuga output locations
+#--------------------------------
+
+# Directory comment
+dir_comment='_'+setup+'_'+polygons+'_'+str(user)
+
 output_dir = join(anuga_dir, 'outputs')+sep
 output_run_time_dir = output_dir+run_time+dir_comment+sep
@@ -86 +192 @@
 
 # Vertex coordinates
-vertex_filename = output_run_time_dir+setup+'vertex_coordinates.txt'
-
-#------------------------------
-# Polygon definitions
-#------------------------------
-
-poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
-res_poly_all = 16000*base_resolution
-
-#print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0
-
-poly_gulleys = read_polygon(polygons_dir+'gulleys_polygon.csv')
-res_gulleys = 3.2*base_resolution
-
-poly_island = read_polygon(polygons_dir+'island_polygon.csv')
-res_island = 32*base_resolution
-
-poly_rhs = read_polygon(polygons_dir+'rhs_polygon.csv')
-res_rhs = 80*base_resolution
-
-interior_regions = [[poly_gulleys,res_gulleys],[poly_island,res_island],
-                    [poly_rhs,res_rhs]]
-                    
-trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
-
-print 'min number triangles', trigs_min
-
-
-
-
-
-
+vertex_filename = output_run_time_dir+'vertex_coordinates.txt'
+
+
+#--------------------------------
+# Area definitions for sww2dem
+#--------------------------------
+
+xminDeep = 0
+xmaxDeep = 1000
+yminDeep = 0
+ymaxDeep = 1360.8
+
+xminMid = 1000
+xmaxMid = 1700
+yminMid = 0
+ymaxMid = 1360.8
+
+xminShallow = 1700
+xmaxShallow = 2179.2
+yminShallow = 0
+ymaxShallow = 1360.8
+
+
+
+
+