Changeset 7546

Timestamp:

Oct 9, 2009, 5:17:27 PM (15 years ago)

Author:

monica

Message:

mandurah tsunami file setup - input elevation and polygons

Location:

anuga_work/production/wa/mandurah_2009

Files:

• build_elevation.py (modified) (3 diffs)
• project.py (modified) (9 diffs)

anuga_work/production/wa/mandurah_2009/build_elevation.py

@@ -64 +64 @@
         geospatial_data[filename] = G_grid.clip(project.bounding_polygon)
 
-# Create Geospatial data from TXT files
+#Create Geospatial data from TXT files
 if not project.point_filenames == []:
     for filename in project.point_filenames:
@@ -79 +79 @@
 print 'Add geospatial objects' # except', project.offshore_name5
 G = None
+G1 = None
+G2 = None
 for key in geospatial_data:
-    #if key == project.point_filenames[4] or key == project.ascii_filenames[1]:
-    #    # Skip these files for now
-    #    continue
+   if key == 'm_peel_1m':
+        G1 = geospatial_data[key].clip(project.clip_box)
+
+   if key == 'ph10mex_clipped':
+        G2 = geospatial_data[key].clip_outside(project.clip_box)
+
+    G += geospatial_data[key]
+G = G + G1 + G2
     
-    G += geospatial_data[key]
-        
-#print 'Clip combined geospatial data'
-##G_clip = G.clip_outside(project.poly_aoi1)
-##G_all = G_clip + geospatial_data[project.point_filenames[4]] 
-#G_clipped = G_all.clip(project.poly_all)
-#G_clip = G.clip(project.bounding_polygon)
-
-
 print 'Export combined DEM file'
 G.export_points_file(project.combined_elevation + '.pts')
@@ -98 +96 @@
 # Use for comparision in ARC
 G.export_points_file(project.combined_elevation + '.txt')
-

anuga_work/production/wa/mandurah_2009/project.py

@@ -17 +17 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'western_australia'
-scenario_name = 'busselton'
-scenario_folder = 'busselton_tsunami_scenario'
+scenario_name = 'mandurah'
+scenario_folder = 'mandurah_tsunami_scenario_2009'
 
 #-------------------------------------------------------------------------------
@@ -34 +34 @@
 finaltime=80000         # final time for simulation
 
-setup = 'final'         # This can be one of three values
+setup = 'trial'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
@@ -50 +50 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number, 'Bt']
+output_comment = [setup, tide, event_number]
 
 #-------------------------------------------------------------------------------
@@ -59 +59 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['busselton_v2',   # Topo
-                        'grid_250m_clip',
-                        'grid_250m_clip2'] # Busselton Topo
+ascii_grid_filenames = [#'m_harvey_1m',   # Topo
+                        'm_peel_1m',
+                        'ph10mex_clipped'] # mandurah Topo
 
 # Format for point is x,y,elevation (with header)
-point_filenames = ['Busselton_Contour0.txt',    # Coastline
-                   'Busselton_BeachSurvey.txt', # Beach survey
-                   'Busselton_NavyFinal.txt',   # Bathymetry
-                   'Busselton_Chart.txt',       # Bathymetry Charts
-                   'Busselton_Digitised.txt',   # Digitised Fairsheet
-                   'Busselton_250m.txt',        # 250m
-                   'Bunbury_TIN.txt',           # Bunbury aoi TIN'd in ArcGIS
-                   'Busselton_TIN.txt',         # Busselton aoi TIN'd in ArcGIS
-                   'XYAHD_clip.txt']            # To extend boundary
+point_filenames = ['bathymetry_all_data.txt', # All data
+                   'coastline.txt']    # Coastline
 
 ### Add csv header list to all files in point_filenames
@@ -82 +75 @@
 # Used in build_elevation.py
 # Format for points easting,northing (no header)
-bounding_polygon_filename = 'bounding_polygon_new.csv'
+bounding_polygon_filename = 'bounding_polygon.csv'
 bounding_polygon_maxarea = 100000
 
@@ -88 +81 @@
 # Used in run_model.py
 # Format for points easting,northing (no header)                    
-interior_regions_data = [['busselton_1km.csv', 500],
-                         ['bunbury_1km.csv', 500],
-                         ['busselton_2km.csv', 10000],
-                         ['bunbury_2km.csv', 10000],
-                         ['island1.csv', 10000],
-                         ['island2.csv', 10000],
-                         ['coast_5km_d20m.csv', 40000]]
+interior_regions_data = [['area_of_interest1.csv', 500],
+                         ['area_of_interest2.csv', 500],
+                         ['area_of_interest3.csv', 500],
+                         ['bounding_area_of_interest.csv', 2500]]
 
 # LAND - used to set the initial stage/water to be offcoast only
 # Used in run_model.py.  Format for points easting,northing (no header)
-land_initial_conditions_filename = [['initial_condition_new.csv', 0],
-                                    ['initial_condition_marina.csv', 0]]
+land_initial_conditions_filename = [['initial_conditions.csv', 0]]
 
 # GAUGES - for creating timeseries at a specific point
 # Used in get_timeseries.py.  
 # Format easting,northing,name,elevation (with header)
-gauges_filename = 'MH_gauges2.csv'
+# gauges_filename = 'MH_gauges2.csv'
 
 # BUILDINGS EXPOSURE - for identifying inundated houses
 # Used in run_building_inundation.py
 # Format latitude,longitude etc (geographic)
-building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
+# building_exposure_filename = 'mandurah_res_clip.csv' # from NEXIS
 
 # AREA OF IMAGES - Extent of each image to find out highest runup
 # Header - easting,northing,id,value
 # Used in get_runup.py
-images_filename = 'define_area_clip.csv'
+# images_filename = 'define_area_clip.csv'
 
 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
@@ -123 +112 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'thinned_boundary_ordering_extend.csv'
+urs_order_filename = 'urs_order.csv'
 
 # Landward bounding points
 # Format easting,northing (no header)
-landward_boundary_filename = 'landward_boundary_new.csv'
+landward_boundary_filename = 'landward_boundary_polygon.csv'
 
 # MUX input filename.
@@ -143 +132 @@
 #-------------------------------------------------------------------------------
 
-# ASCII export grid for Busselton
-xminBusselton = 340000
-xmaxBusselton = 352000
-yminBusselton = 6271500
-ymaxBusselton = 6280000
-
-# ASCII export grid for Bunbury
-xminBunbury = 369000
-xmaxBunbury = 381000
-yminBunbury = 6308000
-ymaxBunbury = 6316500
+# This is my clipping region for topo
+clip_box = 'clip_box.csv'
+
+# ASCII export grid for mandurah
+xminmandurah = 377708
+xmaxmandurah = 380186
+yminmandurah = 6399699
+ymaxmandurah = 6401544
+
 
 ################################################################################
@@ -245 +232 @@
 # The absolute pathname for the gauges file
 # Used for get_timeseries.py
-if gauges_filename:
-    gauges = join(gauges_folder, gauges_filename)       
+#if gauges_filename:
+#    gauges = join(gauges_folder, gauges_filename)       
 
 # The absolute pathname for the building file
 # Used for run_building_inundation.py
-if building_exposure_filename:
-    building_exposure = join(gauges_folder, building_exposure_filename)
+# if building_exposure_filename:
+#    building_exposure = join(gauges_folder, building_exposure_filename)
 
 # The absolute pathname for the image file
 # Used for get_runup.py
-if images_filename:
-    images = join(polygons_folder, images_filename)
+# if images_filename:
+#    images = join(polygons_folder, images_filename)
 
 # full path to where MUX files (or meta-files) live