Changeset 5633
- Timestamp:
- Aug 8, 2008, 3:43:54 PM (15 years ago)
- Location:
- anuga_work/development/landslide_4
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/landslide_4/benchmark_problem4.py
r5572 r5633 39 39 #interior_regions = [[project_slide.poly_syd1, local_res], 40 40 # [project_slide.poly_coast, coast_res]] 41 length = 104. 42 width = 3.7 41 length = 104. #wave tank length 42 width = 3.7 #wave tank width 43 43 44 44 bounding_polygon = [[0,0],[length,0],[length,width],[0,width]] … … 75 75 #------------------------------------------------------------------------------- 76 76 from math import tan 77 depth = 4.6 78 b = 9.144 # 79 def topography(x,y): 77 depth = 4.6 #depth of water domain; needs correction is actually slightly <4.6 78 b = 9.144 #from Synolakis paper 79 def topography(x,y):#sloping intially (gradient=x/2when b<9.144) then flat (depth=4.6m) 80 80 z = [] 81 81 for xi in x: … … 91 91 use_cache = True, 92 92 verbose = True, 93 alpha = 0.1 93 alpha = 0.1 #inclination of domain from strike plane 94 94 ) 95 95 … … 100 100 from math import atan, degrees 101 101 102 a = .4572 # from Synolakis paper103 tsunami_source = slide_tsunami(length=0.91, 104 width=0.61, 105 depth=.01, # de lta from Synolakis paper102 a = .4572 # slide height from Synolakis paper 103 tsunami_source = slide_tsunami(length=0.91,#slide length 104 width=0.61,#slide width 105 depth=.01, # depth below surface to top of slide (delta from Synolakis paper) 106 106 slope=degrees(atan(0.5)), 107 thickness=a, 107 thickness=a, #slide thickness 108 108 x0=0.75, 109 109 y0=width/2, -
anuga_work/development/landslide_4/project.py
r5572 r5633 35 35 sww_name = 'slide' 36 36 37 alpha = 0.1 38 friction=0.01 39 starttime=0 40 finaltime=500 41 remainder_res = 1. 42 export_cellsize=25 37 friction=0.01 #basal friction 38 starttime=0 39 finaltime=500 #duration of evolution 40 remainder_res = 1. #area of mesh triangle in m*2 41 export_cellsize=25 #length along side of mesh triangle 43 42 44 43 dir_comment='_'+'finaltime'+str(finaltime)+'_'+str(user)
Note: See TracChangeset
for help on using the changeset viewer.