Changeset 2682 for production/onslow_2006/export_results.py

Timestamp:

Apr 10, 2006, 9:39:39 AM (19 years ago)

Author:

nick

Message:

updates to onslow
add print outs to least_squares and mesh

File:

• production/onslow_2006/export_results.py (modified) (2 diffs)

production/onslow_2006/export_results.py

@@ -1 +1 @@
 import project, os
+from os import sep
 
 from pyvolution.data_manager import sww2dem
 from pyvolution.ermapper_grids import read_ermapper_grid
 
-#name = project.outputname 
-name = project.outputname
+#name = project.outputname
 
+time_dir = "20060328_094539_alpha_0.1"
+directory = project.outputdir
+name = directory + time_dir +sep + "source"
+
+print'output dir:', name
 #print 'Which variable do you want to export?'
 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
-which_var = 2
+which_var = 4
 
 if which_var == 0:  # Stage
@@ -29 +34 @@
     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
 
+if which_var == 4:  # Elevation
+    outname = name + '_elevation'
+    quantityname = 'elevation'  #Elevation
+
     
 sww2dem(name, basename_out = outname,
             quantity = quantityname,
-            cellsize = 25,       # Trevor would like this at 25
+            cellsize = 30,       # Trevor would like this at 25
             # define region for viz purposes
-            easting_min = project.eminviz,
-            easting_max = project.emaxviz,
-            northing_min = project.nminviz,
-            northing_max = project.nmaxviz,        
+            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, #this is because we want max quantityname
             verbose = True,