Changeset 5590


Ignore:
Timestamp:
Aug 1, 2008, 10:05:08 AM (11 years ago)
Author:
duncan
Message:

Current Hinwood scenario

Location:
anuga_work/development/Hinwood_2008
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/Hinwood_2008/calc_norm.py

    r5577 r5590  
    1414from os import sep, path
    1515from csv import writer
     16from time import localtime, strftime
     17
     18# Related major packages
    1619from Numeric import arange, array, zeros, Float, where, greater, less, \
    1720     compress, argmin, choose, searchsorted, sqrt, sum
    1821
    19 # Related major packages
    2022import project         
    2123from os import sep
     
    142144    if is_interactive:
    143145        ion()
    144        
     146    time_date = strftime('plot date: %d/%m/%Y Time: %H:%M:%S',
     147                                      localtime())
    145148    subplot(212)   
    146     plot_title = id + " Root Mean Square Deviation comparison"
     149    plot_title = id + " Root Mean Square Deviation comparison" + '\n' \
     150                 + time_date
     151                 
    147152    title(plot_title)
    148153    y_label = "RMSD"
    149154    ylabel(y_label)
     155    xlabel("x location, m")
    150156    grid(True)
    151157
     
    164170        rmsd_list = [float(x) for x in simulation['rmsd']]
    165171        lines.append(plot(locations, rmsd_list))
     172       
     173
     174    for break_x in run_data['break_xs']:
     175        axvspan(break_x-0.001,break_x+0.001, facecolor='g')
    166176       
    167177    figlegend(lines, outputdir_tags,'upper left')
     
    273283    #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
    274284    #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
    275     for outputdir_tag in outputdir_tags:
    276         auto_rrms(outputdir_tag, scenarios, "stage")
     285    #for outputdir_tag in outputdir_tags:
     286     #   auto_rrms(outputdir_tag, scenarios, "stage")
    277287   
    278288    #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
    279     #auto_plot_rrms(outputdir_tags, scenarios, "stage")
    280     plot_settings(outputdir_tags, scenarios, "stage")
     289    auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
     290    #plot_settings(outputdir_tags, scenarios, "stage")
  • anuga_work/development/Hinwood_2008/run_dam.py

    r5577 r5590  
    112112        finaltime = None       
    113113        maximum_triangle_area=0.0001
     114        #outputdir_name += '_good'
     115       
     116    elif run_type == 7:
     117        # this is not a test
     118        # Output will go to a file
     119        # The sww file will be interpolated
     120        yieldstep = 0.01
     121        finaltime = None       
     122        maximum_triangle_area=0.00001
     123        #outputdir_name += '_good'
     124       
     125    elif run_type == 8:
     126        # this is not a test
     127        # Output will go to a file
     128        # The sww file will be interpolated
     129        yieldstep = 0.05
     130        finaltime = None       
     131        maximum_triangle_area=0.00001
    114132        #outputdir_name += '_good'
    115133       
     
    232250    Bd = Dirichlet_boundary([0.3,0,0])
    233251    Bts = Time_boundary(domain, function)
     252    #Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
    234253    domain.set_boundary( {'wall': Br, 'wave': Bts} )
    235254    #domain.set_boundary( {'wall': Br, 'wave': Bd} )
     
    287306    # 5 is 0.001
    288307    # 6 is 0.0001
     308    # 7 is 0.00001
    289309   
    290310    #run_type = 1
    291     run_type = 5
     311    run_type = 8
    292312    #for run_data in [scenarios[5]]:
    293313    #scenarios = scenarios[2:]
     
    295315    width = 1.0
    296316    width = 0.1
     317    #width = 0.01
    297318    for run_data in scenarios:
    298319        pro_instance = main( run_data['scenario_id'] + '_boundary.tsm'  ,
     
    302323                             outputdir_name=run_data['scenario_id'],
    303324                             use_limits=False,
    304                              friction=0.012, 
     325                             friction=0.012,
    305326                             end_tag='_G')
    306327        #gauges_for_slope(pro_instance.outputdir,[run_data])
  • anuga_work/development/Hinwood_2008/slope.py

    r5577 r5590  
    674674    #scenarios = [scenarios[0]]
    675675    outputdir_tag = "_good_tri_area_0.01_limiterD"
    676     outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_F"
     676    outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.0001_F"
    677677    slope_tag = ""
    678678    #outputdir_tag = "_test_limiterC"
    679     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
     679    scenarios = [scenarios[4]] # !!!!!!!!!!!!!!!!!!!!!!
    680680    #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
    681681   
    682     gauges_for_slope(slope_tag, outputdir_tag, scenarios)
     682    #gauges_for_slope(slope_tag, outputdir_tag, scenarios)
    683683    #auto_graph_slopes(outputdir_tag, scenarios) #, is_interactive=True)
    684     auto_find_min_slopes(slope_tag, outputdir_tag, scenarios)
     684    #auto_find_min_slopes(slope_tag, outputdir_tag, scenarios)
    685685    #auto_graph_froudes(outputdir_tag, scenarios)
    686     #auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios)
     686    auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios)
    687687    #g = Get_file_name(scenarios[0], outputdir_tag, slope_tag)
    688688    #for wave_file, save_as, wave_number in Get_file_name(
Note: See TracChangeset for help on using the changeset viewer.