Changeset 5714


Ignore:
Timestamp:
Aug 29, 2008, 3:25:34 PM (16 years ago)
Author:
ole
Message:

Handover from Duncan in regard the Hinwood study and the validation paper.
Added comments about locations of data

Files:
4 edited

Legend:

Unmodified
Added
Removed

  • anuga_validation/Hinwood_2008/run_dam.py

    r5713 r5714  
    196196#-------------------------------------------------------------
    197197if __name__ == "__main__":
    198    
     198
     199    # Import scenario metadata for the simulations that is to be run
    199200    from scenarios import scenarios
    200201
    201     #scenarios = [scenarios[0]]
    202    
     202    # Basic parameters of computer flume model (default: w=0.1, max=0.0001, ys=0.01, fr=0.0)
    203203    width = 0.1
    204204    maximum_triangle_area=0.01
    205     yieldstep = 0.01
    206     #yieldstep = 0.5
     205    yieldstep = 0.5
    207206    friction=0.0
    208     isTest=True
    209     #isTest=False
    210    
     207
     208    # Select coarse fast run (isTest is True) or detailed slower run (isTest is False)
     209    #isTest=True
     210    isTest=False
     211
     212    # Loop through experiments and run simulations
    211213    for run_data in scenarios:
    212214        main(run_data,

  • anuga_validation/Hinwood_2008/scenarios.py

    r5712 r5714  
    11"""
    22Scenario information, used to run simulations and create graphs.
     3
     4The Hinwood study consisted of four scenarios each conducted twice, leading to
     58 experiments in total.
    36"""
    47
    58# A list of scenario dictionaries
     9# Each element will represent one experiment (8 in the Hinwood study)
    610scenarios = []
    711
     12
     13
     14# Experiment S1R1
    815offshore_water_depth = .4
    916data = {'xleft':[-3.106,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    … …  
    2330        'break_type':['none', 'plunge', 'plunge', 'plunge', 'plunge', \
    2431                      'plunge'],
    25         #'axis':[0,80,-0.04,0.04],
    26         #'axis_maximum_x':6.0,
    27         #'ANUGA_start_time':0.38,
    28         #'band_offset':-0.5,
    29         'wave_times':[20.0,59.0] # this is in Anuga time
    30         }
    31 scenarios.append(data)
    32 
    33 
     32        'wave_times':[20.0, 59.0]
     33        }
     34scenarios.append(data)
     35
     36# Experiment S1R2
    3437data = {'xleft':[-3.106,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    3538        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    4851        'break_type':['none', 'plunge', 'plunge', 'plunge', 'plunge', \
    4952                      'plunge'],
    50         'wave_times':[23.0,59.0]  # this is in Anuga time
     53        'wave_times':[23.0, 59.0]
    5154                 }
    5255scenarios.append(data)
    5356
    5457
    55 
     58# Experiment S2R1
    5659data = {'xleft':[-4.586,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    5760        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    6770        'start_slope_x':5,
    6871        'finish_slope_x':6.5,
    69         'break_xs':[5.93,5.93,5.907,5.78], # since depth finishes at 65 sec
     72        'break_xs':[5.93, 5.93, 5.907, 5.78], # since depth finishes at 65 sec
    7073        'break_type':['none', 'front steepened', 'front steepened', 'spill'],
    71         'wave_times':[15.0,64.0] # this is in Anuga time
     74        'wave_times':[15.0, 64.0]
    7275                 }
    7376scenarios.append(data)
    7477
     78# Experiment S2R2
    7579data = {'xleft':[-4.586,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    7680        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    8690        'start_slope_x':5,
    8791        'finish_slope_x':6.5,
    88         'break_xs':[6.32,6.32,6.033,5.935],  # since depth finishes at 75 sec
     92        'break_xs':[6.32, 6.32, 6.033, 5.935],  # since depth finishes at 75 sec
    8993        'break_type':['none', 'front steepened', 'front steepened',
    9094                      'weak break'],
    91         'wave_times':[34.0,74.0] # this is in Anuga time
    92         }
    93 scenarios.append(data)
    94 
     95        'wave_times':[34.0, 74.0]
     96        }
     97scenarios.append(data)
     98
     99
     100# Experiment S3R1
     101# Now the slope has changed - and when tank was refilled the depth also changed
    95102offshore_water_depth = .336
    96103data = {'xleft':[-3.875,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    … …  
    107114                   12.569],
    108115
    109         # remove
    110116        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
    111117                               -0.143368,
    … …  
    115121        'start_slope_x':8.,
    116122        'finish_slope_x':9.5,
    117         'break_xs':[9.063,9.063,9.043,9.043],
     123        'break_xs':[9.063, 9.063, 9.043, 9.043],
    118124        'break_type':['none','collapse', 'collapse', 'collapse'],
    119         'wave_times':[30.0,74.0] # this is in Anuga time
     125        'wave_times':[30.0, 74.0]
    120126                 }
    121127scenarios.append(data)
    122128
     129# Experiment S3R2
    123130data = {'xleft':[-3.875,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    124131        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    134141                   12.569],
    135142
    136         # remove
    137143        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
    138144                               -0.143368,
    … …  
    141147        'start_slope_x':8.,
    142148        'finish_slope_x':9.5,
    143         'break_xs':[9.17,9.135,9.135,9.105],
     149        'break_xs':[9.17, 9.135, 9.135, 9.105],
    144150        'break_type':['collapse', 'collapse', 'collapse',
    145151                      'collapse'],
    146         'wave_times':[30.0,85.0] # this is in Anuga time
    147         }
    148 scenarios.append(data)
    149 
     152        'wave_times':[30.0, 85.0]
     153        }
     154scenarios.append(data)
     155
     156# Experiment S4R1
    150157data = {'xleft':[-2.43,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    151158        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    160167                   12.569],
    161168
    162         # remove
    163169        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
    164170                               -0.143368,
    … …  
    168174        'start_slope_x':6,
    169175        'finish_slope_x':9,
    170         # The ANUGA_start_time has not been subtracted from these values
    171         'break_times':[64.29,68.69,73.49,79.39,84.69], # last wave removed
    172         'break_xs':[7.673,7.673,7.43,7.448,7.448], # Pressure info ends early
     176
     177        'break_xs':[7.673, 7.673, 7.43, 7.448, 7.448], # Pressure info ends early
    173178        'break_type':['none','spill', 'spill', 'spill', 'spill'],
    174         'wave_times':[34.0,75.0] # this is in Anuga time
    175         }
    176 scenarios.append(data)
    177 
    178 
     179        'wave_times':[34.0, 75.0]
     180        }
     181scenarios.append(data)
     182
     183# Experiment S4R2
    179184data = {'xleft':[-2.43,0.0 - offshore_water_depth],  # Av' of ADV and Gauge A
    180185        'xtoe':[0.0,0.0 - offshore_water_depth],
    … …  
    190195                   12.569],
    191196
    192         # remove
    193197        'gauge_bed_elevation':[-0.336000, -0.237263, -0.213789, -0.190315,
    194198                               -0.143368,
    … …  
    198202        'start_slope_x':6,
    199203        'finish_slope_x':9,
    200         'break_xs':[8.53,7.46,7.492,7.492,7.444],
     204        'break_xs':[8.53, 7.46, 7.492, 7.492, 7.444],
    201205        'break_type':['collapse', 'spill', 'spill', 'spill', 'spill'],
    202         'wave_times':[34.0,75.0] # this is in Anuga time
    203         }
    204 scenarios.append(data)
    205 
     206        'wave_times':[34.0, 75.0]
     207        }
     208scenarios.append(data)
     209

  • anuga_work/development/Hinwood_2008/run_dam.py

    r5711 r5714  
    373373                             use_limits=False,
    374374                             friction=friction,
    375                              end_tag='_I')
     375                             end_tag='_Z')
    376376        #gauges_for_slope(pro_instance.outputdir,[run_data])

  • anuga_work/development/Hinwood_2008/validation_graphs.py

    r5709 r5714  
    2222
    2323END_TAG = '_I'
    24 #END_TAG = '_Z'
    25 #END_TAG = '_I'
     24END_TAG = '_Ole'
     25
    2626
    2727WIDTH = 0.1
Note: See TracChangeset for help on using the changeset viewer.