Changeset 3997

Timestamp:

Nov 16, 2006, 9:27:20 AM (18 years ago)

Author:

nick

Message:

updates to dampier

Location:

anuga_work/production/dampier_2006

Files:

• project.py (modified) (1 diff)
• run_dampier.py (modified) (3 diffs)

anuga_work/production/dampier_2006/project.py

@@ -204 +204 @@
 poly_facility = read_polygon(polygons_dir+'facility.csv')
 
-poly_pipeline = read_polygon(polygons_dir+'pipeline.csv')
+poly_pipeline = read_polygon(polygons_dir+'pipeline2.csv')
+
+poly_interior = read_polygon(polygons_dir+'interior.csv')
+
+poly_coast = read_polygon(polygons_dir+'coast_final.csv')
 
 clip_poly_e = read_polygon(polygons_dir+'gap_e.csv')

anuga_work/production/dampier_2006/run_dampier.py

@@ -81 +81 @@
 
 if myid == 0:
+    
+    print 'start create mesh from regions'
+    interior_regions = [#[project.karratha_polygon, 25000],
+#                    [project.poly_coast, 10000],
+#                    [project.poly_pipeline, 2000],
+#                    [project.poly_facility, 500]]    
+                    [project.poly_interior, 1000]]       
+    
+    from anuga.utilities.polygon import plot_polygons
+    figname = project.output_run_time_dir + 'poly_pic'
+    plot_polygons([interior_regions,
+               project.bounding_polygon],
+               figname,
+               verbose = True)
 #    if access(project.meshes_time_dir,F_OK) == 0:
 #        mkdir(project.meshes_time_dir)
-    if access(project.meshes_dir,F_OK) == 0:
-        mkdir(project.meshes_dir)
-    print 'start create mesh from regions'
-    interior_regions = [#[project.karratha_polygon, 25000],
-#                    [project.cipma_polygon, 1000],
-                    [project.poly_pipeline, 5000],
-                    [project.poly_facility, 500]]    
+#    if access(project.meshes_dir,F_OK) == 0:
+#        mkdir(project.meshes_dir)
+ 
 #    meshes_dir_name = project.meshes_dir_name + '.msh'
 
@@ -112 +122 @@
 domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
 print domain.statistics()
-
-
 
 print 'starting to create boundary conditions'
@@ -238 +246 @@
     domain.write_boundary_statistics(tags = 'ocean')
     print 'time: ',time.time()
+ 
+    if allclose(t, 6000):
+        domain.set_quantity('xmomentum', 0)
+        domain.set_quantity('ymomentum', 0)
+#import sys; sys.exit()
      
 for t in domain.evolve(yieldstep = 60, finaltime = 28800