- Timestamp:
- Mar 2, 2010, 2:31:30 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/bunbury_storm_surge_2009/project.py
r7613 r7647 46 46 # storm_surge_final - as 'final' but with a longer yieldstep (12 mins) 47 47 48 #------------------------------------------------------------------------------49 #Floodgate parameters50 #------------------------------------------------------------------------------51 52 #changeable parameters specific to the floodgates53 54 floodgate_height = 4.4 #the relief of the floodgates above adjacent topography55 floodgate_thickness = 1 #the thickness of the floodgates in meters56 floodgate_LHS_x = 373546.02 #the x coordinate of the left hand side of the floodgates57 floodgate_LHS_y = 6312332.62 #the y coordinate of the left hand side of the floodgates58 floodgate_RHS_x = 373553.97 #the x coordinate of the right hand side of the floodgates59 floodgate_RHS_y = 6312331.97 #the y coordinate of the right hand side of the floodgates60 61 floodgate_start = 86400 #the time in seconds at which the floodgates begin to close62 floodgate_duration = 10 #the number of seconds that it takes the floodgates to close63 64 #polygon around the floodgates used to define an interior region when building mesh65 floodgate_boundary = [[floodgate_LHS_x,(floodgate_LHS_y - floodgate_thickness)],66 [floodgate_RHS_x,(floodgate_RHS_y - floodgate_thickness)],67 [floodgate_RHS_x,(floodgate_RHS_y + floodgate_thickness)],68 [floodgate_LHS_x,(floodgate_LHS_y + floodgate_thickness)]]69 floodgate_resolution = 1.0 #max triangle area to use when building mesh around floodgates70 48 71 49 #------------------------------------------------------------------------------- … … 91 69 # Used in build_elevation.py 92 70 # Format for ascii grids, as produced in ArcGIS + a projection file 93 ascii_grid_filenames = ['swwa_10m', # LiDAR mosaiced and resampled to 10m 94 'bunbury_5m'] # Bunbury LiDar grid resampled to 5 m - 1m data caused memory errors 95 #~ 'bunbury_nth_a', # 1m LiDAR from WA Dot - nothern quarter (file split as too large for dem2pts) 96 #~ 'bunbury_nth_b', # 1m LiDAR from WA Dot - 2nd northmost quarter (file split as too large for dem2pts) 97 #~ 'bunbury_sth_a', # 1m LiDAR from WA Dot - 2nd southmost quarter (file split as too large for dem2pts) 98 #~ 'bunbury_sth_b'] # 1m LiDAR from WA Dot - southern quarter (file split as too large for dem2pts) 99 # 'bunbury_aoi'] # 1m LiDAR from WA DoT - clipped to area_of_interest.csv 71 ascii_grid_filenames = ['swwa_10m_IC', # LiDAR mosaiced and resampled to 10m, clipped to IC to remove values in water 72 'bunbury_5m_IC'] # Bunbury LiDar grid resampled to 5 m (1m data caused memory errors) 73 # and clipped to Initial Conditions (to remove values in water) 74 # 'bunbury_nth_a', # 1m LiDAR from WA Dot - nothern quarter (file split as too large for dem2pts) 75 # 'bunbury_nth_b', # 1m LiDAR from WA Dot - 2nd northmost quarter (file split as too large for dem2pts) 76 # 'bunbury_sth_a', # 1m LiDAR from WA Dot - 2nd southmost quarter (file split as too large for dem2pts) 77 # 'bunbury_sth_b'] # 1m LiDAR from WA Dot - southern quarter (file split as too large for dem2pts) 78 # 'bunbury_aoi'] # 1m LiDAR from WA DoT - clipped to area_of_interest.csv 100 79 # 'bunbury_nth', # 1m LiDAR from WA Dot - nothern half (file split as too large for dem2pts) 101 80 # 'bunbury_sth'] # 1m LiDAR data from WA DoT - southern half … … 103 82 104 83 # Format for point is x,y,elevation (with header) 105 point_filenames = ['DPI.txt', # Bathymetry data from DPI 106 'Busselton_Chart_Clip_ss.txt', # Clipped from Busselton_Chart - see Busselton Tsunami Scenario 2009 107 'Busselton_NavyFinal_Clip_ss.txt', # Clipped from Busselton_NavyFinal - see Busselton Tsunami Scenario 2009 108 'Leschenault_TIN.txt', # TIN created over the Leschenault Estuary 109 'DPI5U1A02_01a_edited.txt', 'DPI5U1A02_01b_edited.txt', 110 'DPI5U1A02_01c_edited.txt','DPI5U1A02_01d_edited.txt', 111 'DPI5U1A02_01e_edited.txt'] # Bathymetric LiDAR from DPI - split into manageable pieces 84 point_filenames = ['DPI.txt', # Bathymetry data from DPI 85 'Busselton_Chart_Clip_ss.txt', # Clipped from Busselton_Chart - see Busselton Tsunami Scenario 2009 86 'Busselton_NavyFinal_Clip_ss.txt', # Clipped from Busselton_NavyFinal - see Busselton Tsunami Scenario 2009 87 'Leschenault_TIN.txt', # TIN created over the Leschenault Estuary 88 'Leschenault_inlet_TIN.txt' # TIN created over the Leschenault Inlet 89 'DPI5U1A02_01a_edited.txt', # Bathymetric LiDAR from DPI - split into manageable pieces 90 'DPI5U1A02_01b_edited.txt', 91 'DPI5U1A02_01c_edited.txt', 92 'DPI5U1A02_01d_edited.txt', 93 'DPI5U1A02_01e_edited.txt'] 112 94 113 95 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh … … 122 104 interior_regions_data = [['intermediate.csv', 2500], 123 105 ['area_of_interest.csv', 100], 124 [floodgate_boundary, floodgate_resolution]] 106 ['storm_gate_area.csv', 1], 107 ['floodgates.csv', floodgate_resolution]] 125 108 126 109 # LAND - used to set the initial stage/water to be offcoast only
Note: See TracChangeset
for help on using the changeset viewer.