Changeset 3284

Timestamp:

Jul 6, 2006, 2:39:35 PM (18 years ago)

Author:

sexton

Message:

yet more updates

Location:

production/pt_hedland_2006

Files:

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

production/pt_hedland_2006/project.py

@@ -136 +136 @@
 d6 = [690000.0, 7740510.33]
 #polyAll = [d0, d1, d2, d3, d4, d5, d6]
-polyAll = [d0_50, d1_50, d2_50, d4, d5, d6]
+#polyAll = [d0_50, d1_50, d2_50, d4, d5, d6]
+# from Hamish
+h0=[629262.17, 7747205.47]
+h1=[552686.00, 7871579.99] #d3
+h2=[658264.00, 7926314.00] #d2
+h3=[710986.99, 7925796.99] #d1
+h4=[763851.99, 7934357.99] #d0
+h5=[701485.21, 7770656.86]
+h6=[698273.75, 7762227.38]
+h7=[698194.23, 7762018.65]
+h8=[691627.41, 7744781.98]
+h9=[679220.75, 7743604.59]
+h10=[653512.59, 7740528.56]
+h11=[634777.71, 7738247.17]
+h12=[629443.86, 7746910.37]
+h13=[629396.84, 7746986.75]
+h14=[629352.32, 7747059.06]
+h15=[629276.24, 7747182.63]
+h16=[629262.17, 7747205.47] #repeat of h0
+# using Hamish's new bounding polygon
+polyAll = [d0_50, d1_50, d2_50, h16,h15,h14,h13,h12,h11,h10,h9,h8,h7,h6,h5] 
 
 #Interior region - Pt Hedland town

production/pt_hedland_2006/run_pt_hedland.py

@@ -100 +100 @@
 from pmesh.mesh_interface import create_mesh_from_regions
 
-region_res = 100000
-coast_res = 2500
-pt_hedland_res = 1500
+region_res = 50000
+coast_res = 500
+pt_hedland_res = 500
 interior_regions = [[project.poly_pt_hedland, pt_hedland_res],
                     [project.poly_region, region_res]]
@@ -123 +123 @@
                              'left': [2], 'bottomleft': [3],
                              'bottomright': [4], 'right': [5]},
-           'maximum_triangle_area': 500000,
+           'maximum_triangle_area': 250000,
            'filename': meshname,           
            'interior_regions': interior_regions},
@@ -209 +209 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 240, finaltime = 12240): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'topleft')     
-
-for t in domain.evolve(yieldstep = 120, finaltime = 15600
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'topleft')     
-
-for t in domain.evolve(yieldstep = 60, finaltime = 22020
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'topleft')     
+for t in domain.evolve(yieldstep = 240, finaltime = 10800): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'topright')     
+
+for t in domain.evolve(yieldstep = 120, finaltime = 16200
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'topright')     
+
+for t in domain.evolve(yieldstep = 60, finaltime = 21600
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'topright')     
     
-for t in domain.evolve(yieldstep = 120, finaltime = 27060
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'topleft')     
+for t in domain.evolve(yieldstep = 120, finaltime = 27000
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'topright')     
 
 for t in domain.evolve(yieldstep = 240, finaltime = 36000
                        ,skip_initial_step = True): 
     domain.write_time()
-    domain.write_boundary_statistics(tags = 'topleft')     
+    domain.write_boundary_statistics(tags = 'topright')   
   
 print 'That took %.2f seconds' %(time.time()-t0)