Changeset 5789 for anuga_work/production/geraldton/export_results.py

Timestamp:

Sep 25, 2008, 3:14:42 PM (16 years ago)

Author:

kristy

Message:

Updated all scripts to coincide with Perth format

File:

• anuga_work/production/geraldton/export_results.py (modified) (2 diffs)

anuga_work/production/geraldton/export_results.py

@@ -7 +7 @@
 
 
-time_dir = '20080908_173226_run_final_0_elevation_alpha0.1_kvanputt'
+time_dir = '20080911_094915_run_final_0.6_27255_alpha0.1_kvanputt'
+#time_dir = '20080911_100609_run_final_0.6_27283_alpha0.1_kvanputt'
 
 
 cellsize = 25
 #cellsize = 150
-timestep = 0
+#timestep = 0
 directory = project.output_dir
 name = directory+sep+time_dir+sep+project.scenario_name
-
-
-
+#var = [0,4] 
+var = [2] # depth and Speed
+#var = [2,3] # elevation, depth and Speed
 
 is_parallel = False
@@ -26 +27 @@
 
 
-var = [0,4] 
-#var = [2,3] # depth and Speed
-#var = [2,3] # elevation, depth and Speed
+name1 = directory+time_dir+sep+project.scenario_name
+name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_41700_0' #need to get assistance on how to make this into anything
 
+names = [name1, name2]
 
-for which_var in var:
-    if which_var == 0:  # Stage
-        outname = name + '_stage'
-        quantityname = 'stage'
+for name in names:
 
-    if which_var == 1:  # Absolute Momentum
-        outname = name + '_momentum'
-        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'
+            quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
-    if which_var == 3:  # Speed
-        outname = name + '_speed'
-        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
+        if which_var == 2:  # Depth
+            outname = name + '_depth'
+            quantityname = 'stage-elevation'  
 
-    if which_var == 4:  # Elevation
-        outname = name + '_elevation'
-        quantityname = 'elevation'  #Elevation
+        if which_var == 3:  # Speed
+            outname = name + '_speed'
+            quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 
-    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 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,
+                            quantity = quantityname,
+                            #timestep = timestep,
+                            cellsize = cellsize,      
+                            easting_min = project_grad.e_min_area,
+                            easting_max = project_grad.e_max_area,
+                            northing_min = project_grad.n_min_area,
+                            northing_max = project_grad.n_max_area,        
+                            reduction = max, 
+                            verbose = True,
+                            format = 'asc')
+        else:
+            print 'start sww2dem'
+            sww2dem(name, basename_out = outname,
                         quantity = quantityname,
                         #timestep = timestep,
-                        cellsize = cellsize,      
-                        easting_min = project_grad.e_min_area,
-                        easting_max = project_grad.e_max_area,
-                        northing_min = project_grad.n_min_area,
-                        northing_max = project_grad.n_max_area,        
+                        cellsize = cellsize,
+                        number_of_decimal_places = 4,
                         reduction = max, 
                         verbose = True,
                         format = 'asc')
-    else:
-        print 'start sww2dem'
-        sww2dem(name, basename_out = outname,
-                    quantity = quantityname,
-                    timestep = timestep,
-                    cellsize = cellsize,
-                    number_of_decimal_places = 4,
-                    reduction = max, 
-                    verbose = True,
-                    format = 'asc')