Changeset 5572

Timestamp:

Jul 25, 2008, 2:29:44 PM (16 years ago)

Author:

bridgette

Message:

update script for bridgettes use

Location:

anuga_work/development

Files:

• landslide_3 (moved) (moved from anuga_work/development/landslide_benchmark_problem)
• landslide_3/bench3A_t05.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3A_t05.txt)
• landslide_3/bench3A_t10.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3A_t10.txt)
• landslide_3/bench3A_t15.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3A_t15.txt)
• landslide_3/bench3B_t05.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3B_t05.txt)
• landslide_3/bench3B_t10.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3B_t10.txt)
• landslide_3/bench3B_t25.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3B_t25.txt)
• landslide_3/bench3B_t45.txt (copied) (copied from anuga_work/development/landslide_benchmark_problem/bench3B_t45.txt)
• landslide_3/slide.py (copied) (copied from anuga_work/development/landslide_benchmark_problem/slide.py)
• landslide_4/benchmark_problem4.py (modified) (6 diffs)
• landslide_4/export_results.py (added)
• landslide_4/get_timeseries.py (added)
• landslide_4/project.py (added)

anuga_work/development/landslide_4/benchmark_problem4.py

@@ -15 +15 @@
 from os import mkdir, access, F_OK, sep
 import sys
+import project
 
 # Related major packages
 from anuga.shallow_water import Domain, Reflective_boundary, Dirichlet_boundary 
 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
+from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
+
+myid=0
+numprocs=1 
+print 'output_dir',project.output_run_time_dir
+copy_code_files(project.output_run_time_dir,'benchmark_problem4.py', 
+             'project.py' )
+start_screen_catcher(project.output_run_time_dir, myid, numprocs)
 
 #----------------------------------------------------------------------------
@@ -35 +44 @@
 bounding_polygon = [[0,0],[length,0],[length,width],[0,width]]
 meshname = 'slide.msh'
-remainder_res = 1.
+
 
 from caching import cache
@@ -42 +51 @@
            {'boundary_tags': {'bottom': [0], 'right': [1], 
                               'top': [2], 'left': [3]},
-           'maximum_triangle_area': remainder_res,
+           'maximum_triangle_area': project.remainder_res,
            'filename': meshname},
            #'interior_regions': interior_regions},
@@ -57 +66 @@
 print domain.statistics()
 
-domain.set_name('slide')
+domain.set_datadir(project.output_run_time_dir)
+domain.set_name(project.sww_name)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 domain.set_minimum_storable_height(0.01)
@@ -93 +103 @@
 tsunami_source = slide_tsunami(length=0.91,
                                width=0.61,
-                               depth=.1, # delta from Synolakis paper
+                               depth=.01, # delta from Synolakis paper
                                slope=degrees(atan(0.5)),
                                thickness=a, 
@@ -122 +132 @@
 
 t0 = time.time()
-
-for t in domain.evolve(yieldstep = 1, finaltime = 500): 
+finaltime = 500
+for t in domain.evolve(yieldstep = 1, finaltime = project.finaltime): 
     domain.write_time()
     stagestep = domain.get_quantity('stage')