Changeset 5702


Ignore:
Timestamp:
Aug 28, 2008, 9:40:45 AM (16 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/get_timeseries.py

    r5669 r5702  
    1818#timestamp='20080724_121200_run_trial_0.6_polyline_alpha0.1_kvanputt'
    1919#timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
    20 timestamp='20080807_102253_run_final_0.0_polyline_newExtent_kvanputt'
     20timestamp1='20080815_103818_run_final_0_polyline_newExtent_kvanputt'
     21timestamp2='20080815_103708_run_final_0.6_polyline_newExtent_kvanputt'
    2122
     23timestamps = [timestamp1, timestamp2]
     24for timestamp in timestamps:
    2225
    23 filename=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
    24 print 'Hello', filename
     26    filename1=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
     27    filename2=project.output_dir+timestamp+sep+'busselton_time_38340'+sep+project.scenario_name+'_time_38340_0.sww'
    2528
    26 sww2csv_gauges(filename,
    27                 #project.gauges_dir_name,
    28                 project.gauges_dir_name2,
    29                 quantities = ['stage','speed','depth','elevation'],
    30                 verbose=True)
    31                
     29    filenames = [filename1, filename2]
     30    for filename in filenames:
     31
     32        sww2csv_gauges(filename,
     33                        project.gauges_dir_name,
     34                        #project.gauges_dir_name2,
     35                        quantities = ['stage','speed','depth','elevation'],
     36                        verbose=True)
     37                   
    3238
    3339   
  • anuga_work/production/geraldton/build_geraldton.py

    r5655 r5702  
    9595
    9696print'clip combined geospatial object by bounding polygon'
    97 G_clip = G.clip(project.poly_all, verbose=True)
     97G_clip = G.clip(project.poly_all, verbose=False)
    9898
    9999print'export combined DEM file'
  • anuga_work/production/geraldton/project.py

    r5669 r5702  
    9797offshore_in_dir_name3 = topographies_in_dir + offshore_name3
    9898
     99##output topo file location
    99100onshore_dir_name = topographies_dir + onshore_name
    100101island_dir_name = topographies_dir + island_name
  • anuga_work/production/perth/get_timeseries.py

    r5669 r5702  
    2222
    2323filename1=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
    24 #filename2=project.output_dir+timestamp+sep+'perth_time_39900'+sep+project.scenario_name+'_time_39900_0.sww'
     24filename2=project.output_dir+timestamp+sep+'perth_time_39900'+sep+project.scenario_name+'_time_39900_0.sww'
    2525
    2626
    2727
    28 filenames = [filename1]
     28filenames = [filename1, filename2]
    2929for filename in filenames:
    3030
  • anuga_work/production/perth/plot_timeseries.py

    r5430 r5702  
    77import pylab
    88
    9 timestamp='20080619_115643_run_trial_0.6_exmouth_kvanputt'
     9timestamp1='20080815_103336_run_final_0.6_polyline_alpha0.1_kvanputt'
     10timestamp2='20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'
    1011
    11 input_dir_name=project.output_dir+timestamp+sep
    12 print input_dir_name
     12timestamps = [timestamp1, timestamp2]
     13for timestamp in timestamps:
    1314
    14 csv2timeseries_graphs(directories_dic={input_dir_name:['test',0, 0]},
    15                         output_dir=input_dir_name,
    16                         base_name='gauge_',
    17                         plot_numbers='',
    18                         quantities=['stage','speed','depth'],
    19                         extra_plot_name='',
    20                         assess_all_csv_files=True,                           
    21                         create_latex=False,
    22                         verbose=True)
     15    input_dir_name1=project.output_dir+timestamp+sep
     16    input_dir_name2=project.output_dir+timestamp+sep+'perth_time_39900'+sep
    2317
     18    input_dir_names = [input_dir_name1, input_dir_name2]
     19    for input_dir_name in input_dir_names:
     20
     21        print input_dir_name
     22
     23        csv2timeseries_graphs(directories_dic={input_dir_name:['test',0, 0]},
     24                                output_dir=input_dir_name,
     25                                base_name='gauge_',
     26                                plot_numbers='',
     27                                quantities=['stage','speed','depth'],
     28                                extra_plot_name='',
     29                                assess_all_csv_files=True,                           
     30                                create_latex=False,
     31                                verbose=True)
     32
Note: See TracChangeset for help on using the changeset viewer.