Changeset 2954

Timestamp:

May 24, 2006, 9:47:45 AM (18 years ago)

Author:

nick

Message:

update pt hedland

Location:

production/pt_hedland_2006

Files:

• project.py (modified) (4 diffs)
• run_pt_hedland.py (modified) (6 diffs)

production/pt_hedland_2006/project.py

@@ -10 +10 @@
 from coordinate_transforms.redfearn import degminsec2decimal_degrees
 from time import localtime, strftime
+
+from geospatial_data import *
 
 
@@ -43 +45 @@
     user = getenv('LOGNAME')
     print 'USER:', user
-
-#    home = sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation'
-#    home = expanduser('~')
 
 #Derive subdirectories and filenames
@@ -130 +129 @@
 d2 = [658264.0, 7926314.0]
 d3 = [552686.0, 7871580.0]
-d4 = [604415.81, 7733013.56]
-d5 = [656561.15, 7732615.11]
-d6 = [708940.32, 7750510.33]
+#d4 = [604415.81, 7733013.56]
+d4 = [630000.0, 7733013.56]
+#d5 = [656561.15, 7732615.11]
+d5 = [657000.0, 7732615.11]
+#d6 = [708940.32, 7750510.33]
+d6 = [682000.0, 7740510.33]
 
 polyAll = [d0, d1, d2, d3, d4, d5, d6]
@@ -153 +155 @@
 poly_region = [j0, j1, j2, j3]
 
+coast_buffer_file = datadir+'buff_poly_pts_test.xya'
+G = Geospatial_data(file_name=coast_buffer_file,delimiter=' ')
+poly_coast = G.get_data_points()
+#print 'get_data_points()',G.get_data_points()
+
+

production/pt_hedland_2006/run_pt_hedland.py

@@ -29 +29 @@
 from pyvolution.combine_pts import combine_rectangular_points_files 
 from pyvolution.pmesh2domain import pmesh_to_domain_instance
-from geospatial_data import add_points_files
+#from geospatial_data import add_points_files
 
 # Application specific imports
@@ -121 +121 @@
 # derive poly_coast from project.coast_name using alpha_shape
 interior_regions = [[project.poly_pt_hedland, pt_hedland_res],
-                    [project.poly_region, region_res]]
+                    [project.poly_coast, coast_res]]
 
 print 'number of interior regions', len(interior_regions)
@@ -137 +137 @@
         count += 1
 
-figname = 'pt_hedland_polys'
-plot_polygons([project.polyAll,project.poly_pt_hedland,project.poly_region],
-              figname,
-              verbose = True)
+#figname = 'pt_hedland_polys'
+#plot_polygons([project.polyAll,project.poly_pt_hedland,project.poly_region],
+#              figname,
+#              verbose = True)
 
 if count == 0:
@@ -149 +149 @@
     import sys; sys.exit()
     
-
+print 'start create mesh from regions'
 from caching import cache
 _ = cache(create_mesh_from_regions,
@@ -157 +157 @@
           {'boundary_tags': {'topright': [0], 'top': [1],'topleft': [2], 'left': [3],
                              'bottomleft': [4], 'bottomright': [5], 'right': [6]},
-           'maximum_triangle_area': 500000,
+           'maximum_triangle_area': 10000000,
            'filename': meshname,           
            'interior_regions': interior_regions},
@@ -199 +199 @@
                     )
 
-print 'hi1'
 #print 'have sent quantities OK - now exiting'
 #import sys; sys.exit()