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

Timestamp:

Jan 24, 2007, 3:03:26 PM (18 years ago)

Author:

nick

Message:

updates

File:

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

anuga_work/production/dampier_2006/run_dampier.py

@@ -67 +67 @@
 # resolutions (maximal area of per triangle) for each polygon
 #--------------------------------------------------------------------------
+
+poly = [[0,0],[0,100],[100,100],[100,0]]
+
+create_mesh_from_regions(poly,
+                             boundary_tags={'back': [0], 'side': [1,3],
+                                            'ocean': [2]},
+                             maximum_triangle_area=1,
+                             interior_regions=None,
+                             filename=meshes_dir_name,
+                             use_cache=True,
+                             verbose=True)
+
+sys.exit()
 
 if myid == 0:
@@ -90 +103 @@
 print domain.statistics()
 
-"""
+
 print 'starting to create boundary conditions'
 boundaries_in_dir_name = project.boundaries_in_dir_name
@@ -102 +115 @@
     cache(urs2sww,
           (boundaries_in_dir_name,
-    #       boundaries_time_dir_name), 
            boundaries_dir_name), 
           {'verbose': True,
@@ -109 +121 @@
            'minlon': project.west_boundary,
            'maxlon': project.east_boundary,
-#           'minlat': project.south,
-#           'maxlat': project.north,
-#           'minlon': project.west,
-#           'maxlon': project.east,
            'mint': 0, 'maxt': 35100,
            'origin': domain.geo_reference.get_origin(),
@@ -121 +129 @@
            )
 barrier()
-"""
+
 
 #-------------------------------------------------------------------------
@@ -131 +139 @@
 
     from polygon import *
+    #following sets the stage/water to be offcoast only
     IC = Polygon_function( [(project.poly_bathy, tide)], default = 0.)
     domain.set_quantity('stage', IC)
@@ -170 +179 @@
 print 'domain id', id(domain)
 print 'Reading Boundary file'
-#Bf = File_boundary(boundaries_dir_name + '.sww',
-#                  domain, time_thinning=5, use_cache=True, verbose=True)
+Bf = File_boundary(boundaries_dir_name + '.sww',
+                  domain, time_thinning=10, use_cache=True, verbose=True)
 
 print 'finished reading boundary file'
@@ -184 +193 @@
 domain.set_boundary({'back': Br,
                      'side': Bd,
-                     'ocean': Bd}) 
+                     'ocean': Bf}) 
 print'finish set boundary'