def get_timeseries():
"""Get time series data."""
name = os.path.join(project.output_folder, project.scenario_name+'.sww')
log.debug('get_timeseries: input SWW file=%s' % name)
log.debug('get_timeseries: gauge file=%s' % project.gauges)
anuga.sww2csv_gauges(name, project.gauges, quantities=project.layers_list,
verbose=False)
def get_timeseries():
"""Get time series data"
Returns a list of generated files.
"""
# generate the result files
name = os.path.join(project.output_folder, project.scenario+'.sww')
log.debug('get_timeseries: input SWW file=%s' % name)
log.debug('get_timeseries: gauge file=%s' % project.gauge_file)
anuga.sww2csv_gauges(name, project.gauge_file, quantities=project.layers_list,
verbose=False)
# since ANUGA code doesn't return a list of generated files,
# look in output directory for 'gauge_*.csv' files.
glob_mask = os.path.join(project.output_folder, 'gauge_*.csv')
return glob.glob(glob_mask)
def get_timeseries():
"""Get time series data"
Returns a list of generated files.
"""
# check contents of project.gauge_file, if empty, do nothing
with open(project.gauge_file, 'r') as fd:
lines = fd.readlines()
if len(lines) < 2:
log.info('get_timeseries: Nothing to do')
return []
# generate the result files
name = os.path.join(project.output_folder, project.scenario+'.sww')
log.debug('get_timeseries: input SWW file=%s' % name)
log.debug('get_timeseries: gauge file=%s' % project.gauge_file)
anuga.sww2csv_gauges(name, project.gauge_file,
quantities=['stage'], verbose=False)
# since ANUGA code doesn't return a list of generated files,
# look in output directory for 'gauge_*.csv' files.
glob_mask = os.path.join(project.output_folder, 'gauge_*.csv')
return glob.glob(glob_mask)
"""
Generate time series of nominated "gauges" read from project.gauge_filename. This
is done by first running sww2csv_gauges on two different directories to make
'csv' files. Then running csv2timeseries_graphs detailing the two directories
containing the csv file and produces one set of graphs in the 'output_dir' containing
the details at the gauges for both these sww files.
Note, this script will only work if pylab is installed on the platform
"""
from os import sep
import project
import anuga
anuga.sww2csv_gauges('cairns_slide.sww',
project.gauge_filename,
quantities=['stage','speed','depth','elevation'],
verbose=True)
anuga.sww2csv_gauges('cairns_fixed_wave.sww',
project.gauge_filename,
quantities=['stage', 'speed','depth','elevation'],
verbose=True)
try:
import pylab
anuga.csv2timeseries_graphs(directories_dic={'slide'+sep: ['Slide',0, 0],
'fixed_wave'+sep: ['Fixed Wave',0,0]},
output_dir='fixed_wave'+sep,
base_name='gauge_',
plot_numbers='',
quantities=['stage','speed','depth'],
gauge_writer.writerow([xi + X0 + 0.4 * sigmaG, yi * 1.4, 'p3', 0.])
gauge_writer.writerow([xi + X0 + 0.4 * sigmaG, yi * 1.6, 'p4', 0.])
gauge_writer.writerow([xi + X0 + 0.4 * sigmaG, yi * 1.8, 'p5', 0.])
gauge_writer.writerow([xi + X0 + 0.4 * sigmaG, yi * 1.9, 'p6', 0.])
gauge_writer.writerow([xi + X0 + 0.4 * sigmaG, yi * 1.98, 'p7', 0.])
gaugefile.close()
# ------------------------------------------------------------------------------
###STEP-5 Generate Gauge Files
# ------------------------------------------------------------------------------
if p_stepGaugeChart:
os.chdir(prj_output_dir)
anuga.sww2csv_gauges(prj_output_dir + os.sep + sww_file,
prj_output_dir + os.sep + gauge_file,
domain_name + '_gauge_',
quantities=['stage', 'speed', 'elevation'],
verbose=True)
import pylab
anuga.csv2timeseries_graphs(directories_dic={prj_output_dir + os.sep: [r'', 0, 0]},
output_dir=prj_output_dir + os.sep + 'gauges' + os.sep,
base_name=domain_name + '_gauge_',
plot_numbers='',
quantities=['stage', 'speed', 'elevation'],
extra_plot_name='',
assess_all_csv_files=True,
create_latex=False,
verbose=True)
