Changeset 2554
- Timestamp:
- Mar 17, 2006, 9:07:42 AM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
production/onslow_2006/run_onslow.py
r2543 r2554 12 12 13 13 14 #-------------------------------------------------------------------------------# Import necessary modules 15 #------------------------------------------------------------------------------- 14 #------------------------------------------------------------------------------ 15 # Import necessary modules 16 #------------------------------------------------------------------------------ 16 17 17 18 # Standard modules … … 31 32 32 33 33 #------------------------------------------------------------------------------ -34 #------------------------------------------------------------------------------ 34 35 # Preparation of topographic data 35 36 # … … 37 38 # Do for coarse and fine data 38 39 # Fine pts file to be clipped to area of interest 39 #------------------------------------------------------------------------------ -40 #------------------------------------------------------------------------------ 40 41 41 42 # filenames … … 64 65 65 66 # combining the coarse and fine data 66 combine_rectangular_points_files( 67 project.coarsedemname + '.pts', 67 combine_rectangular_points_files(project.coarsedemname + '.pts', 68 68 project.onshore_dem_name + '.pts', 69 69 project.combineddemname + '.pts') 70 70 71 71 72 #------------------------------------------------------------------------------ -72 #------------------------------------------------------------------------------ 73 73 # Create the triangular mesh based on overall clipping polygon with a tagged 74 74 # boundary and interior regions defined in project.py along with 75 75 # resolutions (maximal area of per triangle) for each polygon 76 #------------------------------------------------------------------------------- 77 76 #------------------------------------------------------------------------------ 78 77 from pmesh.mesh_interface import create_mesh_from_regions 79 78 … … 98 97 99 98 100 #------------------------------------------------------------------------------ -99 #------------------------------------------------------------------------------ 101 100 # Setup computational domain 102 #------------------------------------------------------------------------------ -101 #------------------------------------------------------------------------------ 103 102 104 103 domain = pmesh_to_domain_instance(meshname, Domain, … … 115 114 116 115 117 #------------------------------------------------------------------------------ -116 #------------------------------------------------------------------------------ 118 117 # Set up scenario (tsunami_source is a callable object used with set_quantity) 119 #------------------------------------------------------------------------------ -118 #------------------------------------------------------------------------------ 120 119 ''' 121 120 tsunami_source = slump_tsunami(length=30000.0, … … 131 130 132 131 ''' 133 #------------------------------------------------------------------------------ -132 #------------------------------------------------------------------------------ 134 133 # Setup initial conditions 135 #------------------------------------------------------------------------------ -134 #------------------------------------------------------------------------------ 136 135 137 136 tide = 0. … … 150 149 151 150 152 #------------------------------------------------------------------------------ -151 #------------------------------------------------------------------------------ 153 152 # Setup boundary conditions (all reflective) 154 #------------------------------------------------------------------------------ -153 #------------------------------------------------------------------------------ 155 154 156 155 from pyvolution.data_manager import ferret2sww … … 200 199 201 200 202 #------------------------------------------------------------------------------ -201 #------------------------------------------------------------------------------ 203 202 # Evolve system through time 204 #------------------------------------------------------------------------------ -203 #------------------------------------------------------------------------------ 205 204 206 205 import time
Note: See TracChangeset
for help on using the changeset viewer.