Changeset 5572 for anuga_work/development/landslide_4
- Timestamp:
- Jul 25, 2008, 2:29:44 PM (16 years ago)
- Location:
- anuga_work/development/landslide_4
- Files:
-
- 3 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/landslide_4/benchmark_problem4.py
r5568 r5572 15 15 from os import mkdir, access, F_OK, sep 16 16 import sys 17 import project 17 18 18 19 # Related major packages 19 20 from anuga.shallow_water import Domain, Reflective_boundary, Dirichlet_boundary 20 21 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files 22 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters 23 24 myid=0 25 numprocs=1 26 print 'output_dir',project.output_run_time_dir 27 copy_code_files(project.output_run_time_dir,'benchmark_problem4.py', 28 'project.py' ) 29 start_screen_catcher(project.output_run_time_dir, myid, numprocs) 21 30 22 31 #---------------------------------------------------------------------------- … … 35 44 bounding_polygon = [[0,0],[length,0],[length,width],[0,width]] 36 45 meshname = 'slide.msh' 37 remainder_res = 1. 46 38 47 39 48 from caching import cache … … 42 51 {'boundary_tags': {'bottom': [0], 'right': [1], 43 52 'top': [2], 'left': [3]}, 44 'maximum_triangle_area': remainder_res,53 'maximum_triangle_area': project.remainder_res, 45 54 'filename': meshname}, 46 55 #'interior_regions': interior_regions}, … … 57 66 print domain.statistics() 58 67 59 domain.set_name('slide') 68 domain.set_datadir(project.output_run_time_dir) 69 domain.set_name(project.sww_name) 60 70 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 61 71 domain.set_minimum_storable_height(0.01) … … 93 103 tsunami_source = slide_tsunami(length=0.91, 94 104 width=0.61, 95 depth=. 1, # delta from Synolakis paper105 depth=.01, # delta from Synolakis paper 96 106 slope=degrees(atan(0.5)), 97 107 thickness=a, … … 122 132 123 133 t0 = time.time() 124 125 for t in domain.evolve(yieldstep = 1, finaltime = 500):134 finaltime = 500 135 for t in domain.evolve(yieldstep = 1, finaltime = project.finaltime): 126 136 domain.write_time() 127 137 stagestep = domain.get_quantity('stage')
Note: See TracChangeset
for help on using the changeset viewer.