Changeset 1868 for production/karratha_2005/run_karratha.py

Timestamp:

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

Author:

ole

Message:

Ongoing work with karratha

File:

• production/karratha_2005/run_karratha.py (modified) (3 diffs)

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)