Changeset 6511
- Timestamp:
- Mar 13, 2009, 1:25:37 PM (15 years ago)
- Location:
- anuga_work/production/australia_ph2
- Files:
-
- 14 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/LHI/run_model.py
r6502 r6511 80 80 # Reading the landward defined points, this incorporates the original clipping 81 81 # polygon minus the 100m contour 82 ##landward_boundary = read_polygon(project.landward_boundary)82 landward_boundary = read_polygon(project.landward_boundary) 83 83 84 84 # Combine sts polyline with landward points 85 bounding_polygon_sts = event_sts #+ landward_boundary85 bounding_polygon_sts = event_sts + landward_boundary 86 86 87 87 # Number of boundary segments 88 88 num_ocean_segments = len(event_sts) - 1 89 89 # Number of landward_boundary points 90 num_land_points = 0 #file_length(project.landward_boundary)90 num_land_points = file_length(project.landward_boundary) 91 91 92 92 # Boundary tags refer to project.landward_boundary … … 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 164 163 165 164 domain.set_boundary({'back': Br, -
anuga_work/production/australia_ph2/adelaide/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, -
anuga_work/production/australia_ph2/albany/run_model.py
r6484 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 164 163 165 164 domain.set_boundary({'back': Br, -
anuga_work/production/australia_ph2/broome/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, -
anuga_work/production/australia_ph2/carnarvon/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, -
anuga_work/production/australia_ph2/ceduna/run_model.py
r6476 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 164 163 165 164 domain.set_boundary({'back': Br, -
anuga_work/production/australia_ph2/darwin/run_model.py
r6476 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 164 163 165 164 domain.set_boundary({'back': Br, -
anuga_work/production/australia_ph2/esperance/run_model.py
r6476 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 164 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 163 165 164 domain.set_boundary({'back': Br, 166 165 'side': Bd, -
anuga_work/production/australia_ph2/eucla_motel/run_model.py
r6476 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, … … 160 160 161 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd, 163 verbose=True) 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 164 163 165 164 domain.set_boundary({'back': Br, -
anuga_work/production/australia_ph2/gold_coast/project.py
r6496 r6511 28 28 tide = 0 # difference between MSL and HAT 29 29 zone = 56 # specify zone of model 30 event_number = 5 1453 # the event number or the mux file name30 event_number = 50883 # the event number or the mux file name 31 31 alpha = 0.1 # smoothing parameter for mesh 32 32 friction=0.01 # manning's friction coefficient … … 220 220 # full path to where MUX files (or meta-files) live 221 221 mux_input = join(event_folder, mux_input_filename) 222 222 #full path to boundary polygon csv 223 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv') 224 -
anuga_work/production/australia_ph2/hobart/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, -
anuga_work/production/australia_ph2/perth/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True, -
anuga_work/production/australia_ph2/sydney/project.py
r6498 r6511 28 28 # One or all can be changed each time the run_model script is executed 29 29 tide = 0 # difference between MSL and HAT 30 zone = 5 0# specify zone of model31 event_number = 58 326# the event number or the mux file name30 zone = 56 # specify zone of model 31 event_number = 58222 # the event number or the mux file name 32 32 alpha = 0.1 # smoothing parameter for mesh 33 33 friction=0.01 # manning's friction coefficient … … 39 39 wave = 'Bf' # Bf (sts wave) Tb (index wave) 40 40 41 setup = ' final' # This can be one of three values41 setup = 'trial' # This can be one of three values 42 42 # trial - coarsest mesh, fast 43 43 # basic - coarse mesh … … 88 88 # LAND - used to set the initial stage/water to be offcoast only 89 89 # Used in run_model.py. Format for points easting,northing (no header) 90 land_initial_conditions_filename = [['initial_condition_extend.csv', 0], 91 ['initial_condition_comerong_island.csv', 0], 92 ['initial_condition_gabo_island.csv', 0], 93 ['initial_condition_montague_island.csv', 0]] 90 land_initial_conditions_filename = [] 94 91 95 92 # GAUGES - for creating timeseries at a specific point … … 101 98 # Used in run_building_inundation.py 102 99 # Format latitude,longitude etc (geographic) 103 building_exposure_filename = ' busselton_res_clip.csv' # from NEXIS100 building_exposure_filename = 'sydney_res_clip.csv' # from NEXIS 104 101 105 102 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively -
anuga_work/production/australia_ph2/sydney/run_model.py
r6460 r6511 151 151 domain, mean_stage=project.tide, 152 152 time_thinning=1, 153 default_boundary= Bd,153 default_boundary=Dirichlet_boundary([0, 0, 0]), 154 154 boundary_polygon=bounding_polygon_sts, 155 155 use_cache=True,
Note: See TracChangeset
for help on using the changeset viewer.