Changeset 4423 for anuga_work/production/dampier_2006/export_results.py

Timestamp:

Apr 27, 2007, 1:29:59 PM (17 years ago)

Author:

sexton

Message:

updates for export scripts

File:

• anuga_work/production/dampier_2006/export_results.py (modified) (3 diffs)

anuga_work/production/dampier_2006/export_results.py

@@ -5 +5 @@
 from os import sep
 
-time_dir = '20070418_035124_run' #
+time_dir = '20070419_065050_run' # HAT
+#time_dir = '20070419_065018_run' # MSL
 cellsize = 25
 timestep = None
@@ -16 +17 @@
 print 'output dir:', name
 
-which_var = 2
-if which_var == 0:  # Stage
-    outname = name + '_stage'
-    quantityname = 'stage'
+var = [2,3] # depth and speed
 
-if which_var == 1:  # Absolute Momentum
-    outname = name + '_momentum_i1'
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
+for which_var in var:
+    if which_var == 0:  # Stage
+        outname = name + '_stage'
+        quantityname = 'stage'
 
-if which_var == 2:  # Depth
-    outname = name + '_depth'
-    quantityname = 'stage-elevation'  
+    if which_var == 1:  # Absolute Momentum
+        outname = name + '_momentum_i1'
+        quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
-if which_var == 3:  # Speed
-    outname = name + '_speed_i0'
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
+    if which_var == 2:  # Depth
+        outname = name + '_depth_facility'
+        quantityname = 'stage-elevation'  
 
-if which_var == 4:  # Elevation
-    outname = name + '_elevation'
-    quantityname = 'elevation'  #Elevation
+    if which_var == 3:  # Speed
+        outname = name + '_speed_facility'
+        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 
+    if which_var == 4:  # Elevation
+        outname = name + '_elevation'
+        quantityname = 'elevation'  #Elevation
 
-
-if is_parallel == True:
-#    print 'is_parallel',is_parallel
-    for i in range(0,nodes):
-        namei = name + '_P%d_%d' %(i,nodes)
-        outnamei = outname + '_P%d_%d' %(i,nodes)
-        print 'start sww2dem for sww file %d' %(i)
-        sww2dem(namei, basename_out = outnamei,
+    if is_parallel == True:
+    #    print 'is_parallel',is_parallel
+        for i in range(0,nodes):
+            namei = name + '_P%d_%d' %(i,nodes)
+            outnamei = outname + '_P%d_%d' %(i,nodes)
+            print 'start sww2dem for sww file %d' %(i)
+            sww2dem(namei, basename_out = outnamei,
+                        quantity = quantityname,
+                        timestep = timestep,
+                        cellsize = cellsize,      
+                        easting_min = project.e_min_area,
+                        easting_max = project.e_max_area,
+                        northing_min = project.n_min_area,
+                        northing_max = project.n_max_area,        
+                        reduction = max, 
+                        verbose = True,
+                        format = 'asc')
+    else:
+        print 'start sww2dem'
+        sww2dem(name, basename_out = outname,
                     quantity = quantityname,
                     timestep = timestep,
@@ -56 +70 @@
                     verbose = True,
                     format = 'asc')
-else:
-    print 'start sww2dem'
-    sww2dem(name, basename_out = outname,
-                quantity = quantityname,
-                timestep = timestep,
-                cellsize = cellsize,      
-                easting_min = project.e_min_area,
-                easting_max = project.e_max_area,
-                northing_min = project.n_min_area,
-                northing_max = project.n_max_area,        
-                reduction = max, 
-                verbose = True,
-                format = 'asc')