Changeset 5633 for anuga_work/development/landslide_4/benchmark_problem4.py
 Timestamp:
 Aug 8, 2008, 3:43:54 PM (15 years ago)
 File:

 1 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,
Note: See TracChangeset
for help on using the changeset viewer.