Changeset 5449


Ignore:
Timestamp:
Jun 30, 2008, 9:35:07 AM (14 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/plot.py

    r5426 r5449  
    7777    return file_exp, file_sim, location_sims, location_exps, outputdir_name, \
    7878           save_as_list
    79 
     79def plot():
     80    from scenarios import scenarios
     81    outputdir_tag = "_good_tri_area_0.01_A"
     82    outputdir_tag = "_test"
     83    #scenarios = scenarios[1] # !!!!!!!!!!!!!!!!!!!!!!
     84    #scenarios = [scenarios[5]] # !!!!!!!!!!!!!!!!!!!!!!
     85    # Supported formats: emf, eps, pdf, png, ps, raw, rgba, svg, svgz.
     86    plot_type = ".pdf"
     87    for run_data in scenarios:
     88       
     89        temp = Hinwood_files_locations(run_data, outputdir_tag, plot_type)
     90        file_exp, file_sim, location_sims, location_exps, outputdir_name, \
     91                  save_as_list = temp
     92        print "run_data['scenario_id']", run_data['scenario_id']
     93        #location_sims = [location_sims[0]]
     94        #location_exps = [location_exps[0]]
     95        #save_as_list = [save_as_list[0]]
     96        for loc_sim, loc_exp, save_as, gauge in map(None, location_sims,
     97                                                    location_exps,
     98                                                    save_as_list,
     99                                                    run_data['gauge_x']):
     100            time_date = strftime('plot date: %d/%m/%Y Time: %H:%M:%S',
     101                                      localtime())
     102            plot_title = "Scenario: " + outputdir_name + "\n" + \
     103                         "X Gauge (m):" + str(gauge) + "    " + time_date
     104                         
     105            print "Doing ", plot_title
     106            plot_compare_csv(location_sim=loc_sim,
     107                             file_sim=file_sim,
     108                             location_exp=loc_exp,
     109                             file_exp=file_exp,
     110                             plot_title=plot_title,
     111                             y_label='Water height (m)',
     112                             is_interactive=False,
     113                             save_as=save_as)
     114   
    80115#-------------------------------------------------------------
    81116if __name__ == "__main__":
  • anuga_work/development/Hinwood_2008/run_dam.py

    r5447 r5449  
    250250   
    251251    from scenarios import scenarios
     252    from slope import gauges_for_slope
     253    #from plot import plot
    252254
    253255
    254256    #4 is 0.01 5 is 0.001
    255     run_type = 1
     257    run_type = 5
    256258    #for run_data in [scenarios[5]]:
    257259    for run_data in scenarios:
     
    259261              run_type = run_type,
    260262              outputdir_name=run_data['scenario_id'])
    261    
     263    gauges_for_slope()
     264    #plot()
  • anuga_work/development/Hinwood_2008/slope.py

    r5447 r5449  
    104104   
    105105    outputdir_tag = "_good_tri_area_0.01_A"
     106    outputdir_tag = "_test"
    106107    scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
    107108    for run_data in scenarios:
     
    117118def gauges_for_slope():
    118119    from scenarios import scenarios
    119     from scenarios import scenarios
    120120   
    121121    outputdir_tag = "_good_tri_area_0.01_A"
    122     scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
     122    outputdir_tag = "_test"
     123    #scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
    123124    dx = 0.005
    124125    for run_data in scenarios:
Note: See TracChangeset for help on using the changeset viewer.