Changeset 6911


Ignore:
Timestamp:
Apr 28, 2009, 11:10:30 AM (14 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/australia_ph2/bamaga
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/bamaga/project.py

    r6895 r6911  
    1717# Note, the user needs to set up the directory system accordingly
    1818state = 'australia_ph2'
    19 scenario_name = 'gulf'
     19scenario_name = 'weipa'
    2020
    2121#-------------------------------------------------------------------------------
     
    3333##else:
    3434
    35 event_number = 64977    # the event number or the mux file name
     35event_number = 51212    # the event number or the mux file name
    3636
    3737#event_number_list = [31975, 51390, 63735] # To piggy back multiple events
     
    4141friction=0.01           # manning's friction coefficient
    4242starttime=0             # start time for simulation
    43 finaltime=80000 # final time for simulation
    44 setup = 'final'         # This can be one of three values
     43finaltime=1000 # final time for simulation
     44setup = 'trial'         # This can be one of three values
    4545                        #    trial - coarsest mesh, fast
    4646                        #    basic - coarse mesh
     
    4848
    4949# index is only used when wave = Tb
    50 index = 1914            # index from the PTHA - Y2000 0.257m
     50index = 2003            # index from the PTHA - Y2000 0.257m
    5151wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    5252
     
    7979# Used in build_elevation.py
    8080# Format for ascii grids, as produced in ArcGIS + a projection file
    81 ascii_grid_filenames = ['gulf_bg'] # 250m grid 2005
     81ascii_grid_filenames = ['weipa_bg_pr'] # 250m grid 2005
    8282
    8383# Format for points is x,y,elevation (with header)
  • anuga_work/production/australia_ph2/bamaga/run_model.py

    r6895 r6911  
    8787# Number of boundary segments
    8888num_ocean_segments = len(event_sts) - 1
     89print num_ocean_segments
    8990# Number of landward_boundary points
    9091num_land_points = file_length(project.landward_boundary)
    91 
     92print num_land_points
    9293# Boundary tags refer to project.landward_boundary
    9394# 4 points equals 5 segments start at N
    94 ##boundary_tags={'back': range(num_ocean_segments+1,
    95 ##                             num_ocean_segments+num_land_points),
    96 ##               'side': [num_ocean_segments,
    97 ##                        num_ocean_segments+num_land_points],
    98 ##               'ocean': range(num_ocean_segments)}
    99 boundary_tags = {'back':[3],
    100                  'side':[2,4],
    101                  'ocean':[0,1,5]}
     95boundary_tags={'back': range(num_ocean_segments+2,
     96                             num_ocean_segments+5),
     97               'side': [num_ocean_segments,
     98                        num_ocean_segments+1,
     99                        num_ocean_segments+5,
     100                        num_ocean_segments+6,
     101                        num_ocean_segments+7,
     102                        num_ocean_segments+8,
     103                        num_ocean_segments+9],
     104               'ocean': range(num_ocean_segments)}
     105##boundary_tags = {'back':[3],
     106##                 'side':[2,4],
     107##                 'ocean':[0,1,5]}
    102108
    103109print 'boundary tags', boundary_tags
Note: See TracChangeset for help on using the changeset viewer.