Changeset 5633

Timestamp:

Aug 8, 2008, 3:43:54 PM (15 years ago)

Author:

bridgette

Message:

comments for benchmark problem 4

Location:

anuga_work/development/landslide_4

Files:

• benchmark_problem4.py (modified) (4 diffs)
• project.py (modified) (1 diff)

anuga_work/development/landslide_4/benchmark_problem4.py

@@ -39 +39 @@
 #interior_regions = [[project_slide.poly_syd1, local_res],
 #                    [project_slide.poly_coast, coast_res]]
-length = 104.
-width = 3.7
+length = 104. #wave tank length
+width = 3.7 #wave tank width
 
 bounding_polygon = [[0,0],[length,0],[length,width],[0,width]]
@@ -75 +75 @@
 #-------------------------------------------------------------------------------
 from math import tan
-depth = 4.6
-b = 9.144 # from Synolakis paper
-def topography(x,y):
+depth = 4.6 #depth of water domain; needs correction is actually slightly <4.6
+b = 9.144 #from Synolakis paper
+def topography(x,y):#sloping intially (gradient=x/2when b<9.144) then flat (depth=4.6m)
     z = []
     for xi in x:
@@ -91 +91 @@
                     use_cache = True,
                     verbose = True,
-                    alpha = 0.1
+                    alpha = 0.1 #inclination of domain from strike plane
                     )
 
@@ -100 +100 @@
 from math import atan, degrees
 
-a = .4572 # from Synolakis paper
-tsunami_source = slide_tsunami(length=0.91,
-                               width=0.61,
-                               depth=.01, # delta from Synolakis paper
+a = .4572 # slide height from Synolakis paper
+tsunami_source = slide_tsunami(length=0.91,#slide length
+                               width=0.61,#slide width
+                               depth=.01, # depth below surface to top of slide (delta from Synolakis paper)
                                slope=degrees(atan(0.5)),
-                               thickness=a, 
+                               thickness=a, #slide thickness
                                x0=0.75, 
                                y0=width/2, 

anuga_work/development/landslide_4/project.py

@@ -35 +35 @@
 sww_name = 'slide'
 
-alpha = 0.1
-friction=0.01
-starttime=0
-finaltime=500
-remainder_res = 1.
-export_cellsize=25
+friction=0.01 #basal friction
+starttime=0 
+finaltime=500 #duration of evolution
+remainder_res = 1. #area of mesh triangle in m*2
+export_cellsize=25 #length along side of mesh triangle
 
 dir_comment='_'+'finaltime'+str(finaltime)+'_'+str(user)