Ignore:
Timestamp:
Apr 2, 2007, 6:23:29 PM (18 years ago)
Author:
sexton
Message:

update export_results script with updates calculation for speed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/onslow_2006/export_results.py

    r4063 r4347  
    3131#print 'Which variable do you want to export?'
    3232#which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
    33 which_var = 4
     33which_var = 3
    3434#sys.stderr.write(sys.stdout.data)
    3535if which_var == 0:  # Stage
     
    4747if which_var == 3:  # Speed
    4848    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
    5250
    5351if which_var == 4:  # Elevation
     
    6159            cellsize = 25,       
    6260            # 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,       
    6765            reduction = max, #this is because we want max quantityname
    6866            verbose = True,
Note: See TracChangeset for help on using the changeset viewer.