Changeset 6293


Ignore:
Timestamp:
Feb 9, 2009, 9:30:25 AM (15 years ago)
Author:
myall
Message:
 
Location:
anuga_work/production/australia_ph2/ceduna
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/ceduna/build_elevation.py

    r6291 r6293  
    6161                                                verbose=True)
    6262
    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)
     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)
    6868
    6969
  • anuga_work/production/australia_ph2/ceduna/project.py

    r6291 r6293  
    3131# this section needs to be updated to reflect the modelled community.
    3232# 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'
     33state = 'australia_ph2'
     34scenario_name = 'ceduna'
    3635
    3736#------------------------------------------------------------------------------
     
    4039# Model specific parameters. One or all can be changed each time the
    4140# run_scenario script is executed
    42 tide = 0              #0.6
     41tide = 0
    4342event_number = 27255 # Java 9.3 worst case for Perth
    44 #event_number = 68693 # Sumatra 9.2
    45 #event_number = 27283  # Java 9.3 original
    46 tide = 0                #0.6
    47 #event_number = 27255   # Java 9.3 worst case for Perth
    48 #event_number = 68693   # Sumatra 9.2
    49 event_number = 27283    # Java 9.3 original
    5043alpha = 0.1             # smoothing parameter for mesh
    51 friction=0.01           # manning's friction coefficient
    52 starttime=0             
    53 finaltime=80000         # final time for simulation
     44friction = 0.01           # manning's friction coefficient
     45starttime = 0             
     46finaltime = 80000         # final time for simulation
    5447
    55 setup='final'  # Final can be replaced with trial or basic.
     48setup = 'trial'  # Final can be replaced with trial or basic.
    5649               # Either will result in a coarser mesh that will allow a
    5750               # faster, but less accurate, simulation.
     
    8881# Used in build_elevation.py
    8982# Format for ascii grids, as produced in ArcGIS + a projection file
    90 ascii_grid_filenames = ['busselton_v2',   # Topo
    91                         'grid_250m_clip'] # Busselton Topo
     83ascii_grid_filenames = ['grid250m'] # 250m grid 2005
    9284
    9385# 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
    11095
    11196# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     
    11499bounding_polygon_filename = 'bounding_polygon.csv'
    115100
    116 # INTERIOR REGIONS -  for designing the mesh
    117 # Used in run_model.py
    118 # 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 
    127101# GAUGES - for creating timeseries at a specific point
    128102# Used in get_timeseries.py
    129103# 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'
    136105
    137106# BOUNDING POLYGON
     
    143112# Thinned ordering file from Hazard Map (geographic)
    144113# Format is index,latitude,longitude (with header)
    145 urs_order_filename = 'thinned_boundary_ordering_extend.csv'
     114urs_order_filename = 'urs_order.csv'
    146115
    147116# Landward bounding points
    148117# 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 
     118landward_boundary_filename = 'landward_boundary.csv'
    168119
    169120#------------------------------------------------------------------------------
     
    177128# Directory Structure
    178129#------------------------------------------------------------------------------
    179 anuga_folder = join(home, state, scenario_folder, 'anuga')
     130anuga_folder = join(home, state, scenario_name, 'anuga')
    180131topographies_folder = join(anuga_folder, 'topographies')
    181132polygons_folder = join(anuga_folder, 'polygons')
     
    191142# The absolute pathname of the all elevation, generated in build_elevation.py
    192143combined_elevation = join(topographies_folder, combined_elevation_basename)
    193 
    194144
    195145# The absolute pathname of the mesh, generated in run_model.py
     
    216166# The absolute pathname for the gauges file
    217167# 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)       
    223169
    224170#------------------------------------------------------------------------------
     
    226172#------------------------------------------------------------------------------
    227173
    228 # Create list of land polygons with initial conditions
    229 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 factor
    235 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 
    240174# Initial bounding polygon for data clipping
    241175bounding_polygon = read_polygon(join(polygons_folder,
    242176                                     bounding_polygon_filename))
    243177bounding_maxarea = 100000*scale_factor
     178
     179interior_regions = []
    244180
    245181# Estimate the number of triangles                     
Note: See TracChangeset for help on using the changeset viewer.