Changeset 3795

Timestamp:

Oct 16, 2006, 3:59:04 PM (18 years ago)

Author:

sexton

Message:

incorporating new interior region for hobart

Location:

anuga_work/production/hobart_2006

Files:

• project.py (modified) (1 diff)
• run_hobart_clip_refine.py (modified) (1 diff)
• run_hobart_usepts.py (modified) (2 diffs)

anuga_work/production/hobart_2006/project.py

@@ -192 +192 @@
 poly_hobart3 = read_polygon(polygondir+'Hob_poly3.csv')
 poly_hobart4 = read_polygon(polygondir+'Hob_poly4.csv')
+poly_hobart5 = read_polygon(polygondir+'new_interior_region.csv')
 
 plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart4],'boundingpoly2',verbose=False)
 print 'Area of interior regions', (polygon_area(poly_hobart1)+\
-      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+polygon_area(poly_hobart4))/1000000.
+      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+\
+      polygon_area(poly_hobart4)+polygon_area(poly_hobart5))/1000000.
 
 print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.

anuga_work/production/hobart_2006/run_hobart_clip_refine.py

@@ -174 +174 @@
 
 # use 75 for onshore components (12.5m DEM)
-hobart_res = 50000
-bruny_res = 50000
+hobart_res = 7500
+bathy_res = 25000
 interior_regions = [[project.poly_hobart1, hobart_res],
                     [project.poly_hobart2, hobart_res],
                     [project.poly_hobart3, hobart_res],
-                    [project.poly_hobart4, hobart_res]]
+                    [project.poly_hobart5, bathy_res]]
 
 print 'number of interior regions', len(interior_regions)

anuga_work/production/hobart_2006/run_hobart_usepts.py

@@ -146 +146 @@
 
 # use 75 for onshore components (12.5m DEM)
-hobart_res = 50000
-bruny_res = 50000
+hobart_res = 7500
+bathy_res = 50000
 interior_regions = [[project.poly_hobart1, hobart_res],
                     [project.poly_hobart2, hobart_res],
                     [project.poly_hobart3, hobart_res],
-                    [project.poly_hobart4, hobart_res]]
+                    [project.poly_hobart5, bathy_res]]
 
 print 'number of interior regions', len(interior_regions)
@@ -164 +164 @@
                               'e12': [12], 'e13': [13], 'e14': [14],
                               'e15': [15]},
-           'maximum_triangle_area': 250000,
+           'maximum_triangle_area': 500000,
            'filename': meshname,           
            'interior_regions': interior_regions},