Changeset 5633


Ignore:
Timestamp:
Aug 8, 2008, 3:43:54 PM (16 years ago)
Author:
bridgette
Message:

comments for benchmark problem 4

Location:
anuga_work/development/landslide_4
Files:
2 edited

Legend:

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

    r5572 r5633  
    3939#interior_regions = [[project_slide.poly_syd1, local_res],
    4040#                    [project_slide.poly_coast, coast_res]]
    41 length = 104.
    42 width = 3.7
     41length = 104. #wave tank length
     42width = 3.7 #wave tank width
    4343
    4444bounding_polygon = [[0,0],[length,0],[length,width],[0,width]]
     
    7575#-------------------------------------------------------------------------------
    7676from math import tan
    77 depth = 4.6
    78 b = 9.144 # from Synolakis paper
    79 def topography(x,y):
     77depth = 4.6 #depth of water domain; needs correction is actually slightly <4.6
     78b = 9.144 #from Synolakis paper
     79def topography(x,y):#sloping intially (gradient=x/2when b<9.144) then flat (depth=4.6m)
    8080    z = []
    8181    for xi in x:
     
    9191                    use_cache = True,
    9292                    verbose = True,
    93                     alpha = 0.1
     93                    alpha = 0.1 #inclination of domain from strike plane
    9494                    )
    9595
     
    100100from math import atan, degrees
    101101
    102 a = .4572 # from Synolakis paper
    103 tsunami_source = slide_tsunami(length=0.91,
    104                                width=0.61,
    105                                depth=.01, # delta from Synolakis paper
     102a = .4572 # slide height from Synolakis paper
     103tsunami_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)
    106106                               slope=degrees(atan(0.5)),
    107                                thickness=a,
     107                               thickness=a, #slide thickness
    108108                               x0=0.75,
    109109                               y0=width/2,
  • anuga_work/development/landslide_4/project.py

    r5572 r5633  
    3535sww_name = 'slide'
    3636
    37 alpha = 0.1
    38 friction=0.01
    39 starttime=0
    40 finaltime=500
    41 remainder_res = 1.
    42 export_cellsize=25
     37friction=0.01 #basal friction
     38starttime=0
     39finaltime=500 #duration of evolution
     40remainder_res = 1. #area of mesh triangle in m*2
     41export_cellsize=25 #length along side of mesh triangle
    4342
    4443dir_comment='_'+'finaltime'+str(finaltime)+'_'+str(user)
Note: See TracChangeset for help on using the changeset viewer.