Changeset 7658 for anuga_work


Ignore:
Timestamp:
Mar 9, 2010, 7:28:28 PM (14 years ago)
Author:
fountain
Message:

updates to bunbury storm surge model

Location:
anuga_work/production/bunbury_storm_surge_2009
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/bunbury_storm_surge_2009/build_elevation.py

    r7657 r7658  
    9494G9 = geospatial_data[keylist[8]]
    9595G10 = geospatial_data[keylist[9]]
    96 G11 = geospatial_data[keylist[10]]
     96#G11 = geospatial_data[keylist[10]]
    9797
    9898
     
    120120#~ print 'G10', G10.attributes
    121121
    122 G = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8 + G9 + G10 + G11
     122G = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8 + G9 + G10 #+ G11
    123123
    124124# G = None
     
    144144
    145145print 'Export individual text files'
    146 for i in xrange(len(keylist)):
    147     G[i+1].export_points_file(join(project.topographies_folder, keylist[i] +'_export.txt'))
     146#for i in xrange(len(keylist)):
     147#    G[i+1].export_points_file(join(project.topographies_folder, keylist[i] +'_export.txt'))
    148148 
    149 # G1.export_points_file(join(project.topographies_folder, keylist[0] +'_export.txt'))
    150 # G2.export_points_file(join(project.topographies_folder, keylist[1] +'_export.txt'))
    151 # G3.export_points_file(join(project.topographies_folder, keylist[2] +'_export.txt'))
    152 # G4.export_points_file(join(project.topographies_folder, keylist[3] +'_export.txt'))
    153 # G5.export_points_file(join(project.topographies_folder, keylist[4] +'_export.txt'))
    154 # G6.export_points_file(join(project.topographies_folder, keylist[5] +'_export.txt'))
    155 # G7.export_points_file(join(project.topographies_folder, keylist[6] +'_export.txt'))
    156 # G8.export_points_file(join(project.topographies_folder, keylist[7] +'_export.txt'))
    157 # G9.export_points_file(join(project.topographies_folder, keylist[8] +'_export.txt'))
    158 # G10.export_points_file(join(project.topographies_folder, keylist[9] +'_export.txt'))
    159 # G11.export_points_file(join(project.topographies_folder, keylist[10] +'_export.txt'))
     149G1.export_points_file(join(project.topographies_folder, keylist[0] +'_export.txt'))
     150G2.export_points_file(join(project.topographies_folder, keylist[1] +'_export.txt'))
     151G3.export_points_file(join(project.topographies_folder, keylist[2] +'_export.txt'))
     152G4.export_points_file(join(project.topographies_folder, keylist[3] +'_export.txt'))
     153G5.export_points_file(join(project.topographies_folder, keylist[4] +'_export.txt'))
     154G6.export_points_file(join(project.topographies_folder, keylist[5] +'_export.txt'))
     155G7.export_points_file(join(project.topographies_folder, keylist[6] +'_export.txt'))
     156G8.export_points_file(join(project.topographies_folder, keylist[7] +'_export.txt'))
     157G9.export_points_file(join(project.topographies_folder, keylist[8] +'_export.txt'))
     158G10.export_points_file(join(project.topographies_folder, keylist[9] +'_export.txt'))
     159#G11.export_points_file(join(project.topographies_folder, keylist[10] +'_export.txt'))
  • anuga_work/production/bunbury_storm_surge_2009/project.py

    r7657 r7658  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 import anuga.utilities.log as log
     10#import anuga.utilities.log as log
    1111
    1212
     
    3838friction = 0.01           # manning's friction coefficient
    3939starttime = 86400          # start time for simulation - -equivalent to 0000h 4 April 1978
    40 finaltime = 87840 #172800         # final time for simulation - 24 hours for TC Alby
     40finaltime = 172800         # final time for simulation - 24 hours for TC Alby
    4141
    4242setup = 'storm_surge_final'      # This can be one of four values
     
    8585                   'Busselton_Chart_Clip_ss.txt',               # Clipped from Busselton_Chart  - see Busselton Tsunami Scenario 2009
    8686                   'Busselton_NavyFinal_Clip_ss.txt',   # Clipped from Busselton_NavyFinal - see Busselton Tsunami Scenario 2009
    87                    'DPI5U1A02_01a_edited.txt',          # Bathymetric LiDAR from DPI - split into manageable pieces and edited so
    88                    'DPI5U1A02_01b_edited.txt',          # depths below 0 m are negative, and all soundings on land (ie positive)
    89                    'DPI5U1A02_01c_edited.txt',          # are removed as these are not corrected to "bare earth".
    90                    'DPI5U1A02_01d_edited.txt',
    91                    'DPI5U1A02_01e_edited.txt',
    92                    'Leschenault_TIN.txt']               # TIN created over the Leschenault Estuary and Inlet]                                   
     87                   'DPI5U1A02_01a_edited_v2.txt',          # Bathymetric LiDAR from DPI - split into manageable pieces and edited so
     88                   'DPI5U1A02_01b_edited_v2.txt',               # depths below 0 m are negative, and all soundings on land (ie positive)
     89                   'DPI5U1A02_01c_edited_v2.txt',          # are removed as these are not corrected to "bare earth".
     90                   'DPI5U1A02_01d_edited_v2.txt',
     91                   'DPI5U1A02_01e_edited_v2.txt']
     92#                   'Leschenault_TIN.txt']               # TIN created over the Leschenault Estuary and Inlet]                                 
    9393
    9494# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     
    233233
    234234# Log file name
    235 log.log_filename = join(output_run, scenario_name) + '.log'
     235#log.log_filename = join(output_run, scenario_name) + '.log'
    236236
    237237# The absolute pathname for the storm gate file
  • anuga_work/production/bunbury_storm_surge_2009/run_model.py

    r7657 r7658  
    4949from setup_model import project
    5050##import build_urs_boundary as bub
    51 
    52 # Imports for floodgates
    53 from project import floodgate_boundary, floodgate_resolution
    5451
    5552
Note: See TracChangeset for help on using the changeset viewer.