Changeset 2944


Ignore:
Timestamp:
May 23, 2006, 2:24:13 PM (18 years ago)
Author:
sexton
Message:

updates

Location:
production
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/damage.tex

    r2922 r2944  
    11
    2 this will be a discussion of the damage modelling
     2This section deals with modelling the damage to infrastructure as a result
     3of the inundation described in the previous sections.
     4The National Building Exposure Database (NBED) has been
     5created by Geoscience Australia so that consistent risk assessments for a range
     6of natural hazards can be conducted
     7(footnote http://www.ga.gov.au/urban/projects/ramp/NBED.jsp).
     8The NBED contains information
     9about buildings, people, infrastructure, structure value and building contents.
     10It is important to note here that the NBED contains information about
     11residential structures only.
    312
    4 need to talk about potential impact on indigeneous
    5 communities here
     13Once the maximum inundation is calculated for each building, the resultant damage
     14can then be determined as a function of its type and location from the coastline.
    615
     16reference from Ken?
     17
     18Impact on indigeneous communities are important considerations when determining
     19tsunami impact, especially as a number of communities exist in coastal regions.
     20These communities are typcially not included in national residential databases and
     21would be therefore overlooked in damage model estimates.
     22There is one indigeneous community located in this study area, as seen
     23in \ref{fig:communities}.
    724
    825\begin{figure}[hbt]
     
    1330  \label{fig:communities}
    1431\end{figure}
     32
     33\subsection{Highest Astronomical Tide}
     34\subsection{Lowest Astronomical Tide}
  • production/onslow_2006/data.tex

    r2919 r2944  
    3131\begin{tabular}{|l|l|}\hline
    3232Data & Detail \\ \hline
    33 DIGO DTED Level 2  & Onshore, 1 second $\approx$ 30m) \\ \hline
     33DIGO DTED Level 2  & Onshore, 1 second $\approx$ 30m \\ \hline
    3434DLI & Onshore, 20m DEM and orthophotography \\ \hline
    3535DPI & Offshore, fairsheet data around Onslow \\ \hline
  • production/onslow_2006/get_building_inundation.py

    r2896 r2944  
    2929
    3030# Inputs
    31 #timestampdir = '20060426_004129' # HAT
    32 timestampdir = '20060426_004237' # LAT
     31timestampdir = '20060426_004129' # HAT
     32#timestampdir = '20060426_004237' # LAT
    3333file_loc = project.outputdir + timestampdir + sep
    3434swwfile = file_loc + project.basename + '.sww'
     
    3737time_max = None
    3838
    39 def get_buildings_from_file(filename):
     39def get_buildings_from_file(filename,bounding_polygon):
    4040    from coordinate_transforms.redfearn import redfearn
    4141    fid = open(filename)
     
    5454        lon = float(fields[lon_index])
    5555        z, easting, northing = redfearn(lat,lon)
    56         buildings.append([easting, northing])
     56        utm_pt = [easting, northing]
     57        buildings.append(utm_pt)
    5758        loc = fields[name_index]
    5859   
     
    6061
    6162print '\n Buildings obtained from: %s \n' %buildings_filename
    62 buildings, lines = get_buildings_from_file(buildings_filename)
     63buildings, lines = get_buildings_from_file(buildings_filename,project.polyAll)
    6364
    6465sww_quantity = ['stage', 'elevation', 'xmomentum', 'ymomentum']
     
    8889lines[0] = lines[0].strip() +\
    8990           ',MAX INUNDATION DEPTH (m)' +\
    90            ',MAX MOMENTUM (m^2/s)' +\
    91            ',MAX SPEED (m/s) \n'
     91           ',MAX MOMENTUM (m^2/s) \n'
    9292
     93from utilities.polygon import inside_polygon
     94N = len(buildings)
    9395for k, g in enumerate(buildings):
    94     print 'Building %d of %d' %(k, len(buildings))
    95     max_depth = 0
    96     max_momentum = 0
    97     max_velocity = 0     
    98     for i, t in enumerate(T):
     96    if k%((N+10)/10)==0: print 'Building %d of %d' %(k, N)
     97    max_depth = 0.0
     98    max_momentum = 0.0
     99    max_velocity = 0.0
     100    zero_depth = 0.0
     101    zero_momentum = 0.0
     102    if inside_polygon(g,project.polyAll) == True:   
     103        for i, t in enumerate(T):
    99104            w = f(t, point_id = k)[0]
    100105            z = f(t, point_id = k)[1]
     
    103108            depth = w-z
    104109            m = sqrt(uh*uh + vh*vh)   
    105             vel = m / (depth + 1.e-30)
    106             if depth > max_depth: max_depth = w-z
     110            #vel = m / (depth + 1.e-30)
     111            if depth > max_depth: max_depth = depth
    107112            if m > max_momentum: max_momentum = m
    108             if vel > max_velocity: max_velocity = vel
     113    else:
     114        max_depth = 0.0
     115        max_momentum = 0.0
    109116
    110117    lines[k+1] = lines[k+1].strip() +\
    111118                 ',%f' %max_depth +\
    112                  ',%f' %max_momentum +\
    113                  ',%f\n' %max_velocity
     119                 ',%f\n' %max_momentum
     120
    114121
    115122augbuildingsfile = file_loc + 'augmented_buildings.csv'
  • production/onslow_2006/make_report.py

    r2900 r2944  
    5656# User defined inputs
    5757report_title = 'Tsunami impact modelling for the North West shelf: %s' %scenario_name.title()
    58 
    59 #production_dirs = {'20060424_020426_duplictate_time_steps': 'Highest Astronomical Tide',
    60 #                   '20060426_004129': 'Lowest Astronomical Tide'}
    6158
    6259production_dirs = {'20060426_004129': 'Highest Astronomical Tide',
  • production/onslow_2006/project.py

    r2922 r2944  
    5656#for MOST
    5757gauge_filename = gaugedir + 'onslow_gauges1.xya'
    58 buildings_filename = gaugedir + 'all_bld_ind.csv'
     58buildings_filename = gaugedir + 'onslow_res.csv'
    5959community_filename = gaugedir + 'CHINS_v2.csv'
    6060community_scenario = gaugedir + 'community_onslow.csv'
  • production/onslow_2006/results.tex

    r2920 r2944  
    77of the refinement is based around the important inter-tidal zones and
    88other important features such as islands and rivers.
    9 The resultant computational mesh is then seen in \ref{fig.mesh_onslow}.
     9The resultant computational mesh is then seen in \ref{fig:mesh_onslow}.
    1010
    1111\begin{figure}[hbt]
  • production/pt_hedland_2006/data.tex

    r2920 r2944  
    3333\begin{tabular}{|l|l|}\hline
    3434Data & Detail \\ \hline
    35 DIGO DTED Level 2  & Onshore, 1 second $\approx$ 30m) \\ \hline
     35DIGO DTED Level 2  & Onshore, 1 second $\approx$ 30m \\ \hline
    3636DLI & Onshore, 20m DEM and orthophotography \\ \hline
    3737\hline DPI & Offshore, fairsheet data around Onslow \\ \hline
  • production/pt_hedland_2006/project.py

    r2930 r2944  
    6666#for MOST
    6767gauge_filename = gaugedir + 'pt_hedland_gauges.xya'
     68buildings_filename = gaugedir + 'pt_hedland_res.csv'
    6869community_filename = gaugedir + 'CHINS_v2.csv'
    6970community_scenario = gaugedir + 'community_pt_hedland.csv'
  • production/pt_hedland_2006/results.tex

    r2922 r2944  
    77of the refinement is based around the important inter-tidal zones and
    88other important features such as islands and rivers.
    9 The resultant computational mesh is then seen in \ref{fig.mesh_onslow}.
     9The resultant computational mesh is then seen in \ref{fig:mesh_onslow}.
    1010
    1111
Note: See TracChangeset for help on using the changeset viewer.