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

Timestamp:

Sep 26, 2008, 3:46:52 PM (16 years ago)

Author:

kristy

Message:

 

File:

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

anuga_work/production/geraldton/export_results.py

@@ -1 @@
-import project, os
+import project_250m, os
 import sys
 
@@ -7 +7 @@
 
 
-time_dir = '20080911_094915_run_final_0.6_27255_alpha0.1_kvanputt'
-#time_dir = '20080911_100609_run_final_0.6_27283_alpha0.1_kvanputt'
-
+time_dir = '20080925_162836_run_final_0_27283_250m_none_kvanputt'
 
 cellsize = 25
 #cellsize = 150
 #timestep = 0
-directory = project.output_dir
-name = directory+sep+time_dir+sep+project.scenario_name
-#var = [0,4] 
-var = [2] # depth and Speed
+directory = project_250m.output_dir
+name = directory+sep+time_dir+sep+project_250m.scenario_name
+
+var = [0,4] 
+#var = [2] # depth and Speed
 #var = [2,3] # elevation, depth and Speed
 
@@ -26 +25 @@
 print 'output dir:', name
 
+name = directory+time_dir+sep+project_250m.scenario_name
 
-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
+for which_var in var:
+    if which_var == 0:  # Stage
+        outname = name + '_stage'
+        quantityname = 'stage'
 
-names = [name1, name2]
+    if which_var == 1:  # Absolute Momentum
+        outname = name + '_momentum'
+        quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
-for name in names:
+    if which_var == 2:  # Depth
+        outname = name + '_depth'
+        quantityname = 'stage-elevation'  
 
-    for which_var in var:
-        if which_var == 0:  # Stage
-            outname = name + '_stage'
-            quantityname = 'stage'
+    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 == 1:  # Absolute Momentum
-            outname = name + '_momentum'
-            quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
+    if which_var == 4:  # Elevation
+        outname = name + '_elevation'
+        quantityname = 'elevation'  #Elevation
 
-        if which_var == 2:  # Depth
-            outname = name + '_depth'
-            quantityname = 'stage-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 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,
+    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,
-                        number_of_decimal_places = 4,
+                        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,
+                    number_of_decimal_places = 4,
+                    reduction = max, 
+                    verbose = True,
+                    format = 'asc')