Changeset 6293
- Timestamp:
- Feb 9, 2009, 9:30:25 AM (16 years ago)
- Location:
- anuga_work/production/australia_ph2/ceduna
- Files:
-
- 2 edited
-
build_elevation.py (modified) (1 diff)
-
project.py (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/ceduna/build_elevation.py
r6291 r6293 61 61 verbose=True) 62 62 63 # Create Geospatial data from TXT files64 for filename in project.point_filenames:65 absolute_filename = join(project.topographies_folder, filename)66 geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,67 verbose=True)63 ### Create Geospatial data from TXT files 64 ##for filename in project.point_filenames: 65 ## absolute_filename = join(project.topographies_folder, filename) 66 ## geospatial_data[filename] = Geospatial_data(file_name=absolute_filename, 67 ## verbose=True) 68 68 69 69 -
anuga_work/production/australia_ph2/ceduna/project.py
r6291 r6293 31 31 # this section needs to be updated to reflect the modelled community. 32 32 # Note, the user needs to set up the directory system accordingly 33 state = 'western_australia' 34 scenario_name = 'busselton' 35 scenario_folder = 'busselton_tsunami_scenario' 33 state = 'australia_ph2' 34 scenario_name = 'ceduna' 36 35 37 36 #------------------------------------------------------------------------------ … … 40 39 # Model specific parameters. One or all can be changed each time the 41 40 # run_scenario script is executed 42 tide = 0 #0.641 tide = 0 43 42 event_number = 27255 # Java 9.3 worst case for Perth 44 #event_number = 68693 # Sumatra 9.245 #event_number = 27283 # Java 9.3 original46 tide = 0 #0.647 #event_number = 27255 # Java 9.3 worst case for Perth48 #event_number = 68693 # Sumatra 9.249 event_number = 27283 # Java 9.3 original50 43 alpha = 0.1 # smoothing parameter for mesh 51 friction =0.01 # manning's friction coefficient52 starttime =053 finaltime =80000 # final time for simulation44 friction = 0.01 # manning's friction coefficient 45 starttime = 0 46 finaltime = 80000 # final time for simulation 54 47 55 setup ='final' # Final can be replaced with trial or basic.48 setup = 'trial' # Final can be replaced with trial or basic. 56 49 # Either will result in a coarser mesh that will allow a 57 50 # faster, but less accurate, simulation. … … 88 81 # Used in build_elevation.py 89 82 # Format for ascii grids, as produced in ArcGIS + a projection file 90 ascii_grid_filenames = ['busselton_v2', # Topo 91 'grid_250m_clip'] # Busselton Topo 83 ascii_grid_filenames = ['grid250m'] # 250m grid 2005 92 84 93 85 # Format for point is x,y,elevation (with header) 94 point_filenames = ['Busselton_Contour0.txt', # Coastline 95 'Busselton_BeachSurvey.txt', # Beach survey 96 'Busselton_NavyFinal.txt', # Bathymetry 97 'Busselton_Chart.txt', # Bathymetry Charts 98 'Busselton_Digitised.txt', # Digitised Fairsheet 99 'Busselton_250m.txt', # 250m 100 'Bunbury_TIN.txt', # Bunbury aoi TIN'd in ArcGIS 101 'Busselton_TIN.txt', # Busselton aoi TIN'd in ArcGIS 102 'XYAHD_clip.txt'] # To extend boundary 103 104 105 # LAND - used to set the initial stage/water to be offcoast only 106 # Used in run_model,py 107 # Format for points easting,northing (no header) 108 land_initial_conditions_filename = [['initial_condition_extend.csv', 0], 109 ['initial_condition_marina.csv', 0]] 86 ##point_filenames = ['Busselton_Contour0.txt', # Coastline 87 ## 'Busselton_BeachSurvey.txt', # Beach survey 88 ## 'Busselton_NavyFinal.txt', # Bathymetry 89 ## 'Busselton_Chart.txt', # Bathymetry Charts 90 ## 'Busselton_Digitised.txt', # Digitised Fairsheet 91 ## 'Busselton_250m.txt', # 250m 92 ## 'Bunbury_TIN.txt', # Bunbury aoi TIN'd in ArcGIS 93 ## 'Busselton_TIN.txt', # Busselton aoi TIN'd in ArcGIS 94 ## 'XYAHD_clip.txt'] # To extend boundary 110 95 111 96 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh … … 114 99 bounding_polygon_filename = 'bounding_polygon.csv' 115 100 116 # INTERIOR REGIONS - for designing the mesh117 # Used in run_model.py118 # Format for points easting,northing (no header)119 interior_regions_data = [['busselton_1km.csv', 500],120 ['bunbury_1km.csv', 500],121 ['busselton_2km.csv', 10000],122 ['bunbury_2km.csv', 10000],123 ['island1.csv', 10000],124 ['island2.csv', 10000],125 ['coast_5km_d20m.csv', 40000]]126 127 101 # GAUGES - for creating timeseries at a specific point 128 102 # Used in get_timeseries.py 129 103 # Format easting,northing,name,elevation (with header) 130 gauges_filename = 'gauges.csv' 131 132 # BUILDINGS EXPOSURE - for identifying inundated houses 133 # Used in run_building_inundation.py 134 # Format latitude,longitude etc (geographic) 135 building_exposure_filename = 'busselton_res_clip.csv' #from NEXIS 104 ##gauges_filename = 'gauges.csv' 136 105 137 106 # BOUNDING POLYGON … … 143 112 # Thinned ordering file from Hazard Map (geographic) 144 113 # Format is index,latitude,longitude (with header) 145 urs_order_filename = ' thinned_boundary_ordering_extend.csv'114 urs_order_filename = 'urs_order.csv' 146 115 147 116 # Landward bounding points 148 117 # Format easting,northing (no header) 149 landward_boundary_filename = 'landward_boundary_extend.csv' 150 151 #------------------------------------------------------------------------------ 152 # Clipping regions for export to asc and regions for clipping data 153 # Final inundation maps should only be created in regions of the finest mesh 154 #------------------------------------------------------------------------------ 155 156 # ASCII export grid for Busselton 157 xminBusselton = 340000 158 xmaxBusselton = 352000 159 yminBusselton = 6271500 160 ymaxBusselton = 6280000 161 162 # ASCII export grid for Bunbury 163 xminBunbury = 369000 164 xmaxBunbury = 381000 165 yminBunbury = 6308000 166 ymaxBunbury = 6316500 167 118 landward_boundary_filename = 'landward_boundary.csv' 168 119 169 120 #------------------------------------------------------------------------------ … … 177 128 # Directory Structure 178 129 #------------------------------------------------------------------------------ 179 anuga_folder = join(home, state, scenario_ folder, 'anuga')130 anuga_folder = join(home, state, scenario_name, 'anuga') 180 131 topographies_folder = join(anuga_folder, 'topographies') 181 132 polygons_folder = join(anuga_folder, 'polygons') … … 191 142 # The absolute pathname of the all elevation, generated in build_elevation.py 192 143 combined_elevation = join(topographies_folder, combined_elevation_basename) 193 194 144 195 145 # The absolute pathname of the mesh, generated in run_model.py … … 216 166 # The absolute pathname for the gauges file 217 167 # Used for get_timeseries.py 218 gauges = join(gauges_folder, gauges_filename) 219 220 # The absolute pathname for the building file 221 # Used for run_building_inundation.py 222 building_exposure = join(gauges_folder, building_exposure_filename) 168 ##gauges = join(gauges_folder, gauges_filename) 223 169 224 170 #------------------------------------------------------------------------------ … … 226 172 #------------------------------------------------------------------------------ 227 173 228 # Create list of land polygons with initial conditions229 land_initial_conditions = []230 for filename, MSL in land_initial_conditions_filename:231 polygon = read_polygon(join(polygons_folder, filename))232 land_initial_conditions.append([polygon, MSL])233 234 # Create list of interior polygons with scaling factor235 interior_regions = []236 for filename, maxarea in interior_regions_data:237 polygon = read_polygon(join(polygons_folder, filename))238 interior_regions.append([polygon, maxarea*scale_factor])239 240 174 # Initial bounding polygon for data clipping 241 175 bounding_polygon = read_polygon(join(polygons_folder, 242 176 bounding_polygon_filename)) 243 177 bounding_maxarea = 100000*scale_factor 178 179 interior_regions = [] 244 180 245 181 # Estimate the number of triangles
Note: See TracChangeset
for help on using the changeset viewer.
