Changeset 3037
- Timestamp:
- Jun 1, 2006, 10:03:28 AM (18 years ago)
- File:
-
- 1 edited
-
production/pt_hedland_2006/run_pt_hedland.py (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
production/pt_hedland_2006/run_pt_hedland.py
r2955 r3037 1 """Script for running a tsunami inundation scenario for Onslow, WA, Australia.1 """Script for running a tsunami inundation scenario for Port Hedland, WA, Australia. 2 2 3 3 Source data such as elevation and boundary data is assumed to be available in … … 5 5 The output sww file is stored in project.outputtimedir 6 6 7 The scenario is defined by a triangular mesh created from project.polygon ,8 the elevation data and a simulated submarine landslide.9 10 Ole Nielsen and Duncan Gray, GA - 2005 and Nick Bartzis, GA - 20067 The scenario is defined by a triangular mesh created from project.polygon and 8 the elevation data. 9 10 Ole Nielsen and Duncan Gray, GA - 2005, Nick Bartzis and Jane Sexton, GA - 2006 11 11 """ 12 12 13 13 14 #-------------------------------------------------------------------------------# Import necessary modules 15 #------------------------------------------------------------------------------- 14 #------------------------------------------------------------------------------ 15 # Import necessary modules 16 #------------------------------------------------------------------------------ 16 17 17 18 # Standard modules … … 42 43 from pyvolution.util import Screen_Catcher 43 44 44 #------------------------------------------------------------------------------ -45 #------------------------------------------------------------------------------ 45 46 # Preparation of topographic data 46 47 # … … 48 49 # Do for coarse and fine data 49 50 # Fine pts file to be clipped to area of interest 50 #------------------------------------------------------------------------------ -51 #------------------------------------------------------------------------------ 51 52 52 53 # filenames … … 108 109 109 110 ''' 110 #------------------------------------------------------------------------------ -111 #------------------------------------------------------------------------------ 111 112 # Create the triangular mesh based on overall clipping polygon with a tagged 112 113 # boundary and interior regions defined in project.py along with 113 114 # resolutions (maximal area of per triangle) for each polygon 114 #------------------------------------------------------------------------------ -115 #------------------------------------------------------------------------------ 115 116 116 117 from pmesh.mesh_interface import create_mesh_from_regions … … 166 167 167 168 168 #------------------------------------------------------------------------------ -169 #------------------------------------------------------------------------------ 169 170 # Setup computational domain 170 #------------------------------------------------------------------------------- 171 #------------------------------------------------------------------------------ 172 171 173 ''' 172 174 domain = pmesh_to_domain_instance(meshname, Domain, … … 185 187 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 186 188 187 #------------------------------------------------------------------------------ -189 #------------------------------------------------------------------------------ 188 190 # Setup initial conditions 189 #------------------------------------------------------------------------------ -191 #------------------------------------------------------------------------------ 190 192 191 193 tide = 0. … … 205 207 #import sys; sys.exit() 206 208 207 #------------------------------------------------------------------------------ -209 #------------------------------------------------------------------------------ 208 210 # Setup boundary conditions (all reflective) 209 #------------------------------------------------------------------------------ -211 #------------------------------------------------------------------------------ 210 212 print 'start ferret2sww' 211 213 # skipped as results in file SU-AU_clipped is correct for all WA … … 253 255 'bottomleft': Br, 'bottomright': Br, 'right': Br}) 254 256 255 #------------------------------------------------------------------------------ -257 #------------------------------------------------------------------------------ 256 258 # Evolve system through time 257 #------------------------------------------------------------------------------ -259 #------------------------------------------------------------------------------ 258 260 import time 259 261 t0 = time.time()
Note: See TracChangeset
for help on using the changeset viewer.
