Changeset 6597


Ignore:
Timestamp:
Mar 24, 2009, 2:15:16 PM (15 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales/batemans_bay
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r6580 r6597  
    1717# Note, the user needs to set up the directory system accordingly
    1818state = 'new_south_wales'
    19 scenario_name = 'batemans_bay'
    20 scenario_folder = 'batemans_bay_tsunami_scenario_2009'
     19scenario_name = 'gosford'
     20scenario_folder = 'gosford_tsunami_scenario_2009'
    2121
    2222#-------------------------------------------------------------------------------
     
    2626# Model specific parameters.
    2727# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT
    29 event_number = 58368    # the event number or the mux file name
     28tide = 0            # difference between MSL and HAT
     29
     30# the event number or the mux file name
     31##event_number =     #1 in 200 yr Puyesgur
     32##event_number =     #1 in 500 yr
     33##event_number =     #1 in 1000 yr
     34##event_number =     #1 in 2000 yr
     35##event_number =     #1 in 5000 yr
     36event_number = 58242    #1 in 10000 yr Puysegur
     37event_number = 51436    #1 in 10000 yr New Hebrides
     38
    3039alpha = 0.1             # smoothing parameter for mesh
    3140friction=0.01           # manning's friction coefficient
    3241starttime=0             # start time for simulation
    33 finaltime=6000         # final time for simulation
     42finaltime=60000         # final time for simulation
    3443
    3544setup = 'final'         # This can be one of three values
     
    5867# Used in build_elevation.py
    5968# Format for ascii grids, as produced in ArcGIS + a projection file
    60 ascii_grid_filenames = ['1a',# Topographic data
    61                         '1b',
    62                         '2b',
    63                         '2a_3',
    64                         '3b',
    65                         '3a',
    66                         '4a_2',
    67                         '4b',
    68                         'off1',
    69                         'off2',
    70                         'off3',
    71                         'bbhd',
    72                         'sd100031996_p',
    73                         'sd100031996_p2',
    74                         'sd100031996_p3',
    75                         'sd100031996_p4']
     69ascii_grid_filenames = []
    7670                   
    7771# Format for point is x,y,elevation (with header)
    78 point_filenames = ['SD100031996_jgriffin_clip.csv',
    79                    'tomaga_offshore_AHD_MGA_1997.csv',
    80                    'Batemans_BBHD_MGA_1995.csv',
    81                    'moruya_AHD_MGA_2000.csv']
     72point_filenames = []
    8273         
    8374
     
    110101# Used in get_timeseries.py. 
    111102# Format easting,northing,name,elevation (with header)
    112 gauges_filename = 'gauges.csv'
     103##gauges_filename = 'gauges.csv'
    113104
    114105# BUILDINGS EXPOSURE - for identifying inundated houses
     
    245236# The absolute pathname for the gauges file
    246237# Used for get_timeseries.py
    247 gauges = join(gauges_folder, gauges_filename)       
     238##gauges = join(gauges_folder, gauges_filename)       
    248239
    249240# The absolute pathname for the building file
  • anuga_work/production/new_south_wales/batemans_bay/run_model.py

    r6580 r6597  
    161161    print domain.boundary_statistics(tags='ocean')
    162162
    163 ##for t in domain.evolve(yieldstep=5,
    164 ##                       finaltime=project.finaltime,
    165 ##                       skip_initial_step=True):
    166 ##    print domain.timestepping_statistics()
    167 ##    print domain.boundary_statistics(tags='ocean')
     163for t in domain.evolve(yieldstep=20,
     164                       finaltime=6000,
     165                       skip_initial_step=True):
     166    print domain.timestepping_statistics()
     167    print domain.boundary_statistics(tags='ocean')
     168
     169
     170for t in domain.evolve(yieldstep=project.yieldstep,
     171                       finaltime=project.finaltime,
     172                       skip_initial_step=True):
     173    print domain.timestepping_statistics()
     174    print domain.boundary_statistics(tags='ocean')
    168175
    169176print 'Simulation took %.2f seconds' % (time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.