Changeset 4069

Timestamp:

Dec 7, 2006, 5:02:56 PM (18 years ago)

Author:

sexton

Message:

nsw slide modelling updates

Location:

anuga_work/production

Files:

• dampier_2006/export_results.py (modified) (5 diffs)
• newcastle_2006/project_slide.py (modified) (3 diffs)
• newcastle_2006/run_newcastle_slide.py (modified) (1 diff)
• wollongong_2006/run_gong_slide.py (modified) (1 diff)

anuga_work/production/dampier_2006/export_results.py

@@ -20 +20 @@
 
 if which_var == 1:  # Absolute Momentum
-    outname = name + '_momentum_i50'
+    outname = name + '_momentum_i1'
     quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
@@ -28 +28 @@
 
 if which_var == 3:  # Speed
-    outname = name + '_speed_all'
+    outname = name + '_speed_i0'
     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
@@ -37 +37 @@
 
 if is_parallel == True:
-    for i in range(2,nodes):
+    for i in range(0,nodes):
         namei = name + '_P%d_%d' %(i,nodes)
         outnamei = outname + '_P%d_%d' %(i,nodes)
@@ -43 +43 @@
         sww2dem(namei, basename_out = outnamei,
                     quantity = quantityname,
-                    timestep = 50,
+                    timestep = 1,
                     cellsize = 100,      
                     #easting_min = project.e_min_area,
@@ -49 +49 @@
                     #northing_min = project.n_min_area,
                     #northing_max = project.n_max_area,        
-                    #reduction = max, 
+                    reduction = max, 
                     verbose = True,
                     format = 'asc')

anuga_work/production/newcastle_2006/project_slide.py

@@ -126 +126 @@
 bulli_density = 1.46
 bulli_slope = 4.0
-bulli_alpha = 126.0 + 90.0
+bulli_alpha = 126.0 - 90.0
 
 shovel_depth = 968.0
@@ -134 +134 @@
 shovel_density = 1.49
 shove_slope = 4.0
-shovel_alpha = 118.0 + 90.0
+shovel_alpha = 118.0 - 90.0
 
 yacaaba_depth = 1119.0
@@ -142 +142 @@
 yacaaba_density = 1.48
 yacaaba_slope = 3.0
-yacaaba_alpha = 133.0 + 90.0
+yacaaba_alpha = 133.0 - 90.0

anuga_work/production/newcastle_2006/run_newcastle_slide.py

@@ -169 +169 @@
 import time
 t0 = time.time()
+from Numeric import allclose
+from anuga.abstract_2d_finite_volumes.quantity import Quantity
 
-for t in domain.evolve(yieldstep = 30, finaltime = 480): 
+for t in domain.evolve(yieldstep = 30, finaltime = 5000): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'e14')

anuga_work/production/wollongong_2006/run_gong_slide.py

@@ -172 +172 @@
 import time
 t0 = time.time()
+from Numeric import allclose
+from anuga.abstract_2d_finite_volumes.quantity import Quantity
 
-for t in domain.evolve(yieldstep = 30, finaltime = 480): 
+for t in domain.evolve(yieldstep = 30, finaltime = 5000): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'e14')