Changeset 3661 for anuga_work/production/hobart_2006/export_results.py

Timestamp:

Sep 25, 2006, 2:26:13 PM (17 years ago)

Author:

sexton

Message:

clean up

File:

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

anuga_work/production/hobart_2006/export_results.py

@@ -3 +3 @@
 
 from anuga.shallow_water.data_manager import sww2dem
+#from anuga.pyvolution.ermapper_grids import read_ermapper_grid
+from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
 from os import sep
 
-time_dir = '20060921_083846' #MSL
+time_dir = '20060922_024842' #MSL
 directory = project.outputdir
-name = directory + time_dir + sep + 'source'
+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 = 4 # for Hamish/Alex to compare
-#which_var = 2 # for Ingo/Neil to make maps
-
+which_var = 4
+#sys.stderr.write(sys.stdout.data)
 if which_var == 0:  # Stage
     outname = name + '_stage'
@@ -33 +35 @@
 
 if which_var == 4:  # Elevation
-    outname = name + '_elevation'
+    outname = name + '_elevation_node3'
     quantityname = 'elevation'  #Elevation
 
+print 'start sww2dem'
+#sys.stderr.write(sys.stdout.data)
 sww2dem(name, basename_out = outname,
             quantity = quantityname,
-            cellsize = 25,       
+            cellsize = 25,      # Trevor would like this at 25
             # define region for viz purposes
             easting_min = project.e_min_area,
@@ -47 +51 @@
             verbose = True,
             format = 'asc')
+
+#sys.stderr.write(sys.stdout.data)
+
+#Check
+
+#data = read_ermapper_grid(name)
+#print 'Values from %s are in [%f, %f]' %(name, min(data.flat), max(data.flat))