Changeset 1868

Timestamp:

Oct 5, 2005, 5:40:34 PM (19 years ago)

Author:

ole

Message:

Ongoing work with karratha

Location:

production/karratha_2005

Files:

• create_mesh.py (modified) (2 diffs)
• project.py (modified) (2 diffs)
• run_karratha.py (modified) (3 diffs)

production/karratha_2005/create_mesh.py

@@ -92 +92 @@
     """
 
-    from pmesh.mesh import *
+    from pmesh.mesh import Mesh
     from pyvolution.coordinate_transforms.redfearn import redfearn
     from pyvolution.util import populate_polygon
@@ -107 +107 @@
 
 
+
+    
     print "***********************"
     print "geo ref", geo

production/karratha_2005/project.py

@@ -63 +63 @@
 
 #Interior regions
-karratha_south = degminsec2decimal_degrees(-20,45,0)
-karratha_north = degminsec2decimal_degrees(-20,37,0)
-karratha_west = degminsec2decimal_degrees(116,45,0)
-karratha_east = degminsec2decimal_degrees(116,55,0)
+karratha_south = degminsec2decimal_degrees(-20,44,0)
+karratha_north = degminsec2decimal_degrees(-20,42,0)
+karratha_west = degminsec2decimal_degrees(116,48,0)
+karratha_east = degminsec2decimal_degrees(116,53,30)
 
 k0 = [karratha_south, karratha_west]
@@ -75 +75 @@
 karratha_polygon = [k0, k1, k2, k3]
 
+
+#Interior regions
+dampier_south = degminsec2decimal_degrees(-20,40,0)
+dampier_north = degminsec2decimal_degrees(-20,38,10)
+dampier_west = degminsec2decimal_degrees(116,43,0)
+dampier_east = degminsec2decimal_degrees(116,45,0)
+
+d0 = [dampier_south, dampier_west]
+d1 = [dampier_south, dampier_east]
+d2 = [dampier_north, dampier_east]
+d3 = [dampier_north, dampier_west]    
+
+dampier_polygon = [d0, d1, d2, d3]
+
+
+#Interior regions
+refinery_south = degminsec2decimal_degrees(-20,38,0)
+refinery_north = degminsec2decimal_degrees(-20,36,0)
+refinery_west = degminsec2decimal_degrees(116,44,0)
+refinery_east = degminsec2decimal_degrees(116,46,10)
+
+d0 = [refinery_south, refinery_west]
+d1 = [refinery_south, refinery_east]
+d2 = [refinery_north, refinery_east]
+d3 = [refinery_north, refinery_west]    
+
+refinery_polygon = [d0, d1, d2, d3]

production/karratha_2005/run_karratha.py

@@ -11 +11 @@
 """
 
-tide = 0.75   #HMWS estimate by Colin French, GA
+#tide = 0.75   #HMWS estimate by Colin French, GA
+tide = 0       #HMW 
 
 
@@ -83 +84 @@
 from create_mesh import create_mesh
 
-interior_regions = [[project.karratha_polygon, 10000]]
+interior_regions = [[project.karratha_polygon, 25000],
+                    [project.dampier_polygon, 8000],
+                    [project.refinery_polygon, 8000]]
+
 m = cache(create_mesh,
           project.polygon,
           {'boundary_tags': {'back': [7, 8], 'side': [0, 6], 'ocean': [1, 2, 3, 4, 5]},
-          'resolution': 80000,
+          'resolution': 100000,
           'filename': project.meshname + '.msh',
           'interior_regions': interior_regions},
@@ -142 +146 @@
 
 
-#Run
-#for t in domain.evolve(yieldstep = 600, finaltime = 15000):
-#    domain.write_time()
-#    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
-#
-#for t in domain.evolve(yieldstep = 10, finaltime = 35000):
-#    domain.write_time()
-#    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
+#Evolve
+import time
+t0 = time.time()
 
-for t in domain.evolve(yieldstep = 60, finaltime = 40000):
+for t in domain.evolve(yieldstep = 60, finaltime = 15000):
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
+
+for t in domain.evolve(yieldstep = 20, finaltime = 35000,
+                       skip_initial_step = True):
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
+
+for t in domain.evolve(yieldstep = 60, finaltime = 40000,
+                       skip_initial_step = True):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')       
+
+print 'That took %.2f seconds' %(time.time()-t0)