# Changeset 5568

Ignore:
Timestamp:
Jul 24, 2008, 4:05:04 PM (15 years ago)
Message:

incorporated definitions from web and Synolakis paper

File:
1 edited

Unmodified
Added
Removed
• ## anuga_work/development/landslide_4/benchmark_problem4.py

 r5567 #interior_regions = [[project_slide.poly_syd1, local_res], #                    [project_slide.poly_coast, coast_res]] length = 10. width = 10. depth = 10. bounding_polygon = [[0,width],[0,0],[length,0],[length,width]] length = 104. width = 3.7 bounding_polygon = [[0,0],[length,0],[length,width],[0,width]] meshname = 'slide.msh' remainder_res = 1. from caching import cache _ = cache(create_mesh_from_regions, bounding_polygon, {'boundary_tags': {'top': [0], 'left': [1], 'bottom': [2], 'right': [3]}, {'boundary_tags': {'bottom': [0], 'right': [1], 'top': [2], 'left': [3]}, 'maximum_triangle_area': remainder_res, 'filename': meshname}, #------------------------------------------------------------------------------- from math import tan depth = 4.6 b = 9.144 # from Synolakis paper def topography(x,y): z = -depth + tan(0.5)*x z = [] for xi in x: if xi < b: z.append(-xi/2.) else: z.append(-depth) #print 'hello', xi,b, z return z domain.set_quantity('stage', 0.0) #------------------------------------------------------------------------------- from smf import slide_tsunami from math import atan, degrees tsunami_source = slide_tsunami(length=2., width=1., depth=1., slope=45, thickness=0.1, x0=length/2, a = .4572 # from Synolakis paper tsunami_source = slide_tsunami(length=0.91, width=0.61, depth=.1, # delta from Synolakis paper slope=degrees(atan(0.5)), thickness=a, x0=0.75, y0=width/2, alpha=0., #dx=0.001, domain=domain) #verbose=True) #dx=0.001, domain=domain, verbose=True) #------------------------------------------------------------------------------- Br = Reflective_boundary(domain) Bd = Dirichlet_boundary([tide,0,0]) Bd = Dirichlet_boundary([0,0,0]) domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd} ) domain.set_boundary( {'top': Bd,  'bottom': Bd, 'right': Bd, 'left': Bd} ) t0 = time.time() for t in domain.evolve(yieldstep = 1, finaltime = 10): for t in domain.evolve(yieldstep = 1, finaltime = 500): domain.write_time() stagestep = domain.get_quantity('stage')
Note: See TracChangeset for help on using the changeset viewer.