Changeset 4347


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

update export_results script with updates calculation for speed

Location:
anuga_work/production
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/broome_2006/export_results.py

    r4063 r4347  
    2626if which_var == 3:  # Speed
    2727    outname = name + '_speed'
    28     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     28    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    2929
    3030if which_var == 4:  # Elevation

  • anuga_work/production/dampier_2006/export_results.py

    r4069 r4347  
    2929if which_var == 3:  # Speed
    3030    outname = name + '_speed_i0'
    31     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    32     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     31    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    3332
    3433if which_var == 4:  # Elevation

  • anuga_work/production/hobart_2006/export_results.py

    r3944 r4347  
    2929if which_var == 3:  # Speed
    3030    outname = name + '_speed'
    31     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     31    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    3232
    3333if which_var == 4:  # Elevation

  • anuga_work/production/newcastle_2006/export_results.py

    r4346 r4347  
    6666if which_var == 3:  # Speed
    6767    outname = name + '_speed_i0'
    68     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    69     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     68    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    7069
    7170if which_var == 4:  # Elevation

  • 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,

  • anuga_work/production/pt_hedland_2006/export_results.py

    r4063 r4347  
    4444if which_var == 3:  # Speed
    4545    outname = name + '_speed'
    46     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    47     #quantityname = 'xmomentum/(stage-elevation)'  #Speed
    48     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     46    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    4947
    5048if which_var == 4:  # Elevation

  • anuga_work/production/sydney_2006/export_results.py

    r4346 r4347  
    4040if which_var == 3:  # Speed
    4141    outname = name + '_speed_i0'
    42     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    43     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     42    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    4443
    4544if which_var == 4:  # Elevation

  • anuga_work/production/wollongong_2006/export_results.py

    r4346 r4347  
    4949if which_var == 3:  # Speed
    5050    outname = name + '_speed_i0'
    51     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    52     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     51    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    5352
    5453if which_var == 4:  # Elevation
Note: See TracChangeset for help on using the changeset viewer.