source: anuga_work/production/perth/get_timeseries.py @ 6342

Last change on this file since 6342 was 6342, checked in by myall, 15 years ago
File size: 2.0 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 anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
16import project
17
18directory = project.output_dir
19
20##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
21##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
22##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
23time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
24##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
25##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
26
27time_dirs = [time_dir4] #[time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
28
29for time_dir in time_dirs:
30    name = directory+time_dir+sep+project.scenario_name
31    gauge = project.gauges_dir_name
32    sww2csv_gauges(name+'.sww',gauge,
33                   quantities = ['stage','speed','depth','elevation'],
34                   verbose=True)
35
36   
37##    name1 = directory+time_dir+sep+project.scenario_name
38####    name2 = directory+time_dir+sep+project.scenario_name+'_time_43380_0'
39####    name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
40##
41##    names = [name1]#, name2, name3]
42##
43##    for i, name in enumerate(names):
44##        gauge = project.gauges_dir_name
45####        out_name = 'gauge_%d_' % i
46##        print name
47##        sww2csv_gauges(name+'.sww',gauge,#out_name = out_name,
48##                        quantities = ['stage','speed','depth','elevation'],
49##                        verbose=True)
50##
Note: See TracBrowser for help on using the repository browser.