Changeset 4282 for anuga_work/production/dampier_2006/run_dampier.py

Timestamp:

Feb 27, 2007, 9:39:30 AM (17 years ago)

Author:

nick

Message:

update to dampier

File:

• anuga_work/production/dampier_2006/run_dampier.py (modified) (11 diffs)

anuga_work/production/dampier_2006/run_dampier.py

@@ -34 +34 @@
 from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier
 from anuga_parallel.parallel_abstraction import get_processor_name
+from anuga.caching import myhash
 # Application specific imports
 import project                 # Definition of file names and polygons
@@ -81 +82 @@
 
 sys.exit()
-'''
+
 if myid == 0:
     
     print 'start create mesh from regions'
+#    print 'project.res_poly_all',project.res_poly_all
+#    print 'project.interior_regions_test',project.interior_regions_test
+#    print 'meshes_dir_name',meshes_dir_name
     create_mesh_from_regions(project.poly_all,
                              boundary_tags={'back': [2,3], 'side': [0, 1, 4],
                                             'ocean': [5]},
                              maximum_triangle_area=project.res_poly_all,
-                             interior_regions=project.interior_regions,
+                             interior_regions=project.interior_regions_test,
                              filename=meshes_dir_name,
                              use_cache=True,
@@ -96 +100 @@
 # to sync all processors are ready 
 barrier()
+'''
+create_mesh_from_regions(project.poly_all,
+                             boundary_tags={'back': [2,3], 'side': [0, 1, 4],
+                                            'ocean': [5]},
+                             maximum_triangle_area=project.res_poly_all,
+                             interior_regions=project.interior_regions_test,
+                             filename=meshes_dir_name,
+                             use_cache=True,
+                             verbose=True)
 
 #-------------------------------------------------------------------------
@@ -101 +114 @@
 #-------------------------------------------------------------------------
 print 'Setup computational domain'
-domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
+from caching import cache
+domain = cache(Domain, (meshes_dir_name), {'use_cache':True, 'verbose':True}, verbose=True)
+print 'domain id', id(domain) 
+print 'myhash', myhash(domain)      
+       
 print domain.statistics()
+
 boundaries_dir_name=project.boundaries_dir_name
 
@@ -108 +126 @@
 
 from anuga.shallow_water.data_manager import urs2sww
-
+'''
 # put above distribute
 print 'boundary file is: ',project.boundaries_dir_name
@@ -129 +147 @@
            )
 barrier()
-
+'''
 
 #-------------------------------------------------------------------------
@@ -148 +166 @@
 #                    filename = project.combined_dir_name + '.pts',
 # MUST USE TXT FILES FOR CACHING TO WORK!
-                    filename = project.combined_dir_name + '.txt',
-                    use_cache = False,
+                    filename = project.combined_smaller_dir_name + '.txt',
+                    use_cache = True,
                     verbose = True,
                     alpha = 0.1)
@@ -166 +184 @@
 # Set domain parameters
 #------------------------------------------------------
-
+print 'domain id', id(domain)
 domain.set_name(project.scenario_name)
 domain.set_datadir(project.output_run_time_dir)
@@ -174 +192 @@
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
+print 'domain id', id(domain)
+domain.beta_h = 0
+domain.limit2007 = 1
 
 #-------------------------------------------------------------------------
@@ -180 +201 @@
 print 'Available boundary tags', domain.get_boundary_tags()
 print 'domain id', id(domain)
-print 'Reading Boundary file'
-Bf = File_boundary(boundaries_dir_name + '.sww',
+print 'Reading Boundary file',project.boundaries_dir_name4 + '.sww'
+#Bf = File_boundary(boundaries_dir_name + '.sww',
+Bf = File_boundary(project.boundaries_dir_name4 + '.sww',
                   domain, time_thinning=4, use_cache=True, verbose=True)
 
@@ -204 +226 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 120, finaltime = 9000):
+for t in domain.evolve(yieldstep = 60, finaltime = 9000):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'ocean')