Changeset 4066

Timestamp:

Dec 7, 2006, 9:49:27 AM (18 years ago)

Author:

nick

Message:

updates to dampier

Location:

anuga_work/production/dampier_2006

Files:

• project.py (modified) (2 diffs)
• run_dampier.py (modified) (3 diffs)

anuga_work/production/dampier_2006/project.py

@@ -160 +160 @@
 #west_boundary = degminsec2decimal_degrees(116,0,0)
 #east_boundary = degminsec2decimal_degrees(118,00,0)
-west_boundary = degminsec2decimal_degrees(116,10,0)
-east_boundary = degminsec2decimal_degrees(117,15,0)
+west_boundary = degminsec2decimal_degrees(116,00,0)
+east_boundary = degminsec2decimal_degrees(117,20,0)
 
 
@@ -178 +178 @@
 #bounding_polygon, zone =\
 #                  convert_from_latlon_to_utm([p1, p2, p3, p4, p5, p6, p7])
+print bounding_polygon
 refzone = zone
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon

anuga_work/production/dampier_2006/run_dampier.py

@@ -70 +70 @@
                  dirname(project.__file__)+sep+ project.__name__+'.py' )
 barrier()
-start_screen_catcher(project.output_run_time_dir, myid, numprocs)
+#start_screen_catcher(project.output_run_time_dir, myid, numprocs)
 
 print 'USER: ', project.user
@@ -89 +89 @@
                     [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([project.poly_coast,project.poly_pipeline,project.poly_facility,
+               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)
+    import sys; sys.exit()
+    print 'start create mesh from regions'
+    interior_regions = [#[project.karratha_polygon, 25000],
+#                    [project.cipma_polygon, 1000],
+#                    [project.poly_pipeline, 5000],
+#                    [project.poly_facility, 500]]    
+                    [project.poly_interior, 1000]]    
+#    meshes_dir_name = project.meshes_dir_name + '.msh'
 
     create_mesh_from_regions(project.bounding_polygon,
@@ -194 +211 @@
 domain.set_store_vertices_uniquely(False)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
-domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
+domain.set_maximum_allowed_speed(0.0) # Allow a little runoff (0.1 is OK)
 
 #-------------------------------------------------------------------------