Changeset 4347 for anuga_work/production/onslow_2006
- Timestamp:
- Apr 2, 2007, 6:23:29 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/onslow_2006/export_results.py
r4063 r4347 31 31 #print 'Which variable do you want to export?' 32 32 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3 ' )) 33 which_var = 433 which_var = 3 34 34 #sys.stderr.write(sys.stdout.data) 35 35 if which_var == 0: # Stage … … 47 47 if which_var == 3: # Speed 48 48 outname = name + '_speed' 49 #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5' #Speed 50 #quantityname = 'xmomentum/(stage-elevation)' #Speed 51 quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)' #Speed 49 quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))' #Speed 52 50 53 51 if which_var == 4: # Elevation … … 61 59 cellsize = 25, 62 60 # define region for viz purposes 63 #easting_min = project.e_min_area,64 #easting_max = project.e_max_area,65 #northing_min = project.n_min_area,66 #northing_max = project.n_max_area,61 easting_min = project.e_min_area, 62 easting_max = project.e_max_area, 63 northing_min = project.n_min_area, 64 northing_max = project.n_max_area, 67 65 reduction = max, #this is because we want max quantityname 68 66 verbose = True,
Note: See TracChangeset
for help on using the changeset viewer.