Changeset 3671 for anuga_work/production/hobart_2006/export_results.py
- Timestamp:
- Sep 27, 2006, 5:55:02 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2006/export_results.py
r3668 r3671 3 3 4 4 from anuga.shallow_water.data_manager import sww2dem 5 #from anuga.pyvolution.ermapper_grids import read_ermapper_grid6 5 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher 7 6 from os import sep 8 7 9 time_dir = '2006092 5_115139' #MSL8 time_dir = '20060927_054931' #MSL 10 9 directory = project.outputdir 11 10 name = directory + time_dir +sep + 'source' … … 15 14 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3 ' )) 16 15 which_var = 4 17 #sys.stderr.write(sys.stdout.data)18 16 if which_var == 0: # Stage 19 17 outname = name + '_stage' … … 30 28 if which_var == 3: # Speed 31 29 outname = name + '_speed' 32 #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5' #Speed33 #quantityname = 'xmomentum/(stage-elevation)' #Speed34 30 quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)' #Speed 35 31 36 32 if which_var == 4: # Elevation 37 outname = name + '_elevation _node3'33 outname = name + '_elevation' 38 34 quantityname = 'elevation' #Elevation 39 35 40 36 print 'start sww2dem' 41 #sys.stderr.write(sys.stdout.data)42 37 sww2dem(name, basename_out = outname, 43 38 quantity = quantityname, 44 cellsize = 30, # would prefer this at 2539 cellsize = 25, # would prefer this at 25 45 40 # define region for viz purposes 46 41 #easting_min = project.e_min_area, … … 52 47 format = 'asc') 53 48 54 #sys.stderr.write(sys.stdout.data)55 56 #Check57 58 #data = read_ermapper_grid(name)59 #print 'Values from %s are in [%f, %f]' %(name, min(data.flat), max(data.flat))
Note: See TracChangeset
for help on using the changeset viewer.