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

"""
Generate time series of nominated "gauges".
Input: sww file from run_perth.py
       gauges project.gauge_filename
Output: csv files stage,speed,depth,elevation over time
Stored in the 'outputs_dir' folder for respective .sww file

Note:
Can run different sww files at the same time
make sure if there is a second sww file that it is placed into a folder called sww2
Can run different gauges at the same time - ie testing boundary index point
"""

from os import sep, rename, listdir, system
from os.path import join
from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
import project

directory = project.output_folder

##time_dir1 = '20090601_133348_run_final_0.8_58280_None_kvanputt'
##time_dir2 = '20090601_133503_run_final_0.8_58292_None_kvanputt'
##time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt'
##time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt'
##time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt'
##time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'
##
##time_dirs = [time_dir3] #3, time_dir4, time_dir5, time_dir6]

time_dir1 = '20090619_165107_run_final_0_58260_None_kvanputt'
time_dir2 = '20090619_165136_run_final_0.8_58260_None_kvanputt'

time_dirs = [time_dir2] #, time_dir2]

print 'time directories', time_dirs

for time_dir in time_dirs:
    name = join(directory,time_dir,project.scenario_name)
    gauge = project.gauges
    sww2csv_gauges(name+'.sww',gauge,
                   quantities = ['stage','speed','depth','elevation'],
                   verbose=True)