Changeset 1920 for production/karratha_2005/run_karratha.py

Timestamp:

Oct 14, 2005, 11:33:57 AM (18 years ago)

Author:

ole

Message:

Karratha

File:

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

production/karratha_2005/run_karratha.py

@@ -62 +62 @@
        'origin': project.mesh_origin,
        'mean_stage': tide,
-       'zscale': 2,                 #Enhance tsunami
+       'zscale': 1,                 #Enhance tsunami
        'fail_on_NaN': False,
        'inverted_bathymetry': True},
+      #evaluate = True,
       verbose = True)
 #FIXME: Dependencies
@@ -110 +111 @@
 domain.store = True
 
-#domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
-domain.quantities_to_be_stored = ['stage']
+domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
+#domain.quantities_to_be_stored = ['stage']
         
 print 'Number of triangles = ', len(domain)
@@ -141 +142 @@
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
-domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} )
+
+domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} ) #MOST tsunami
 
 #domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bd} )  #Nothing
@@ -150 +152 @@
 t0 = time.time()
 
-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 = 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 = 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):
+for t in domain.evolve(yieldstep = 60, finaltime = 40000):
+                       #skip_initial_step = True):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')