source: anuga_work/production/busselton/get_timeseries.py @ 5714

Last change on this file since 5714 was 5702, checked in by kristy, 17 years ago
File size: 1.5 KB
Line 
1"""
2Generate time series of nominated "gauges" read from project.gauge_filename. This
3is done by first running sww2csv_gauges on two different directories to make
4'csv' files. Then running csv2timeseries_graphs detailing the two directories
5containing the csv file and produces one set of graphs in the 'output_dir' containing
6the details at the gauges for both these sww files.
7
8Note, this script will only work if pylab is installed on the platform
9
10"""
11
12from os import sep
13from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
14
15import project
16
17
18#timestamp='20080724_121200_run_trial_0.6_polyline_alpha0.1_kvanputt'
19#timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
20timestamp1='20080815_103818_run_final_0_polyline_newExtent_kvanputt'
21timestamp2='20080815_103708_run_final_0.6_polyline_newExtent_kvanputt'
22
23timestamps = [timestamp1, timestamp2]
24for timestamp in timestamps:
25
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'
28
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                   
38
39   
40
Note: See TracBrowser for help on using the repository browser.