source: anuga_work/development/alpha_validation/find_alpha.py @ 4072

Last change on this file since 4072 was 3970, checked in by nick, 18 years ago

start development of automatic alpha validation

File size: 1.0 KB
Line 
1
2from abstract_2d_finite_volumes import domain
3
4
5
6create_mesh_from_regions(project.bounding_polygon,
7                         boundary_tags={'back': [7, 8], 'side': [0, 6],
8                                        'ocean': [1, 2, 3, 4, 5]},
9                         maximum_triangle_area=100000,
10                         interior_regions=interior_regions,
11#                         filename=meshes_time_dir_name,
12                         filename=meshes_dir_name,
13                         use_cache=True,
14                         verbose=True)
15
16
17
18alphas = [0.001, 0.01, 0.1, 1.0, 10.0, 100.0]
19domains = {}
20print 'Setup computational domain'
21for alpha in alphas:
22   
23    domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
24    print domain.statistics()
25    domain.set_quantity('elevation', 
26                    filename = project.topographies_dir + build_time + sep + project.combined_final_name + '.pts',
27                    use_cache = True,
28                    verbose = True,
29                    alpha = alpha)
Note: See TracBrowser for help on using the repository browser.