Changeset 1920

Timestamp:

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

Author:

ole

Message:

Karratha

Location:

production/karratha_2005

Files:

• export_results.py (modified) (1 diff)
• project.py (modified) (1 diff)
• run_karratha.py (modified) (4 diffs)

production/karratha_2005/export_results.py

@@ -5 +5 @@
 
 #name = project.outputname + '_0.0tide_notsunami'
-name = project.outputname + '_0.0tide' 
+#name = project.outputname + '_0.0tide'
 
-sww2dem(name, basename_out = name+'_crop',
-        quantity = 'depth',
-        cellsize = 10,
-        easting_min = 465646,
-        easting_max = 480286,
-        northing_min = 7712865,
-        northing_max = 7724233,        
+name = project.outputname
+
+sww2dem(name, basename_out = name + '_momentum',
+        #quantity = 'stage-elevation',  #Depth
+        quantity = '(xmomentum**2 + ymomentum**2)**0.5',  #Absolute momentum
+        cellsize = 50,
+        #cellsize = 10,        
+        #easting_min = 465646,
+        #easting_max = 480286,
+        #northing_min = 7712865,
+        #northing_max = 7724233,        
         reduction = max,
         verbose = True,

production/karratha_2005/project.py

@@ -7 +7 @@
 import sys
 
-#Making assumptions abot the location of scenario data
+#Making assumptions about the location of scenario data
 scenario_dir_name = 'karratha_tsunami_scenario_2005'
-basename = 'karratha_100m'
+#basename = 'karratha_100m'
+basename = 'elevation50m'
 boundary_basename = 'SU-AU'
 

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')