Changeset 4601 for anuga_work/production/shark_bay_2007/run_shark_bay.py

Timestamp:

Jul 6, 2007, 2:40:57 PM (18 years ago)

Author:

ole

Message:

Work on shark bay

File:

• anuga_work/production/shark_bay_2007/run_shark_bay.py (modified) (5 diffs)

anuga_work/production/shark_bay_2007/run_shark_bay.py

@@ -74 +74 @@
     if myid == 0:
         copy_code_files(kwargs['output_dir'],__file__, 
-                 dirname(project.__file__)+sep+ project.__name__+'.py' )
+                        dirname(project.__file__)+sep+project.__name__+'.py' )
 
         store_parameters(**kwargs)
@@ -97 +97 @@
         print 'start create mesh from regions'
         
-        # FIXME (Ole): What if tags are wrong?
         create_mesh_from_regions(project.bounding_polygon,
-                                 boundary_tags={'back': [1, 2, 3, 4, 5],
-                                                'side': [0, 6],
-                                                'ocean': [7, 8, 9, 10, 11]},
+                                 boundary_tags=project.boundary_tags,
                                  maximum_triangle_area=project.res_bounding_polygon,
                                  interior_regions=project.interior_regions,
@@ -185 +182 @@
                         time_thinning=time_thinning,
                         mean_stage=tide, 
-                        use_cache=True,
+                        use_cache=False,
                         verbose=True)
                     
@@ -197 +194 @@
     print'set_boundary'
 
-    domain.set_boundary({'back': Br,
-                         'side': Bd,
-                         'ocean': Bf}) 
+    #domain.set_boundary({'back': Br,
+    #                     'side': Bd,
+    #                     'ocean': Bf})
+    domain.set_boundary({'tide': Bd,
+                         'ocean': Bf})     
     print'finish set boundary'
 
@@ -211 +210 @@
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')     
+
+
+    # Reset stage in study area
+    print 'Reset initial condition in study area'
+    domain.set_quantity('stage', 0.0,
+                        polygon=project.poly_tsunami_approach_area)
+    
 
     for t in domain.evolve(yieldstep = kwargs['yieldstep'], finaltime = 21600