source: anuga_work/production/hobart_2009/get_timeseries.py @ 7083

Last change on this file since 7083 was 7082, checked in by kristy, 16 years ago
File size: 1.4 KB
Line 
1"""
2Generate time series of nominated "gauges".
3Input: sww file from run_perth.py
4       gauges project.gauge_filename
5Output: csv files stage,speed,depth,elevation over time
6Stored in the 'outputs_dir' folder for respective .sww file
7
8Note:
9Can run different sww files at the same time
10make sure if there is a second sww file that it is placed into a folder called sww2
11Can run different gauges at the same time - ie testing boundary index point
12"""
13
14from os import sep, rename, listdir, system
15from os.path import join
16from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
17import project
18
19directory = project.output_folder
20
21time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
22time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'
23time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt'
24time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt'
25time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt'
26time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'
27
28time_dirs = [time_dir1, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6]
29 
30print 'time directories', time_dirs
31
32for time_dir in time_dirs:
33    name = join(directory,time_dir,project.scenario_name)
34    gauge = project.gauges
35    sww2csv_gauges(name+'.sww',gauge,
36                   quantities = ['stage','speed','depth','elevation'],
37                   verbose=True)
38 
39
Note: See TracBrowser for help on using the repository browser.