Changeset 3809

Timestamp:

Oct 17, 2006, 12:21:30 PM (19 years ago)

Author:

sexton

Message:

introduce refined grids around specified sites for Hobart

Location:

anuga_work/production/hobart_2006

Files:

• project.py (modified) (5 diffs)
• run_hobart_clipdata_refine.py (added)

anuga_work/production/hobart_2006/project.py

@@ -124 +124 @@
 #outputname = outputtimedir + basename  #Used by post processing
 
-
 ###############################
 # Domain definitions
@@ -146 +145 @@
 polyAll = read_polygon(polygondir+'new_extent_2.csv')
 plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
-print 'Area of bounding polygon', polygon_area(polyAll)
+#print 'Area of bounding polygon', polygon_area(polyAll)
 polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
-print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
-print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
+#print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
+#print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
 
 polyAll_refine2 = read_polygon(polygondir+'export_refine_2.csv')
-print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
-print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
+#print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
+#print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
 
 ###################################################################
@@ -164 +163 @@
 northingmin = 5245000
 northingmax = 5260000
-
+# for original model setup
 # clipping 25m data set - Hobart
 eastingmin25 = 524208.387
@@ -183 +182 @@
 northingmax25_3 = 5258511.857
 
+# export for refined grids at paleo sites
+eastingmin_site13 = 532835.94
+eastingmax_site13 = 533848.23
+northingmin_site13 = 5237384.18
+northingmax_site13 = 5238195.03
+
+eastingmin_kingston = 525905.49
+eastingmax_kingston = 527773.42
+northingmin_kingston = 5240256.71
+northingmax_kingston = 5242148.58
+
+eastingmin_bruny = 523638.9
+eastingmax_bruny = 533000.88
+northingmin_bruny = 5204967.41
+northingmax_bruny = 5212177.69
+
 ###############################
 # Interior region definitions
@@ -194 +209 @@
 poly_hobart5 = read_polygon(polygondir+'new_interior_region.csv')
 
-plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'boundingpoly2',verbose=False)
-print 'Area of interior regions', (polygon_area(poly_hobart1)+\
-      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.
-print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.
+plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'orig_model',verbose=False)
+
+# to refine around paleo sites
+poly_kingston = read_polygon(polygondir+'Kingston.csv')
+poly_bruny = read_polygon(polygondir+'bruny_refine.csv')
+poly_site13 = read_polygon(polygondir+'Site13.csv')
+
+plot_polygons([polyAll, poly_kingston,poly_bruny,poly_site13,poly_hobart5],'orig_model',verbose=False)
+#print 'Area of interior regions', (polygon_area(poly_hobart1)+\
+#      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.
+#print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.