source: anuga_work/production/hobart_2006/get_timeseries.py @ 3852

Last change on this file since 3852 was 3852, checked in by sexton, 18 years ago

updated results based on refined mesh

File size: 1.6 KB
Line 
1"""
2Generate images of "gauges" for production run
3
4Inputs:
5
6production dirs: dictionary of production directories with a
7                 association to that simulation run, eg high tide,
8                 low tide and MSL.
9                   
10Outputs:
11
12* figures used for report stored in the report_figure directory
13
14"""
15
16from os import getcwd, sep, altsep, mkdir, access, F_OK
17import project
18from anuga.abstract_2d_finite_volumes.util import sww2timeseries
19
20# Create report directory
21reportdir = getcwd()+sep+'report'+sep
22if access(reportdir,F_OK) == 0:
23    mkdir (reportdir)
24   
25# 2500 interior resolution for both events
26# use both to compare the events
27production_dirs = {'20061006_062319': 'Mw 8.5',
28                   '20061008_234702': 'Mw 8.7'}
29
30# Generate figures
31swwfiles = {}
32for label_id in production_dirs.keys():
33   
34    file_loc = project.outputdir + label_id + sep
35    swwfile = file_loc + project.basename + '.sww'
36    swwfiles[swwfile] = label_id
37
38texname, elev_output = sww2timeseries(swwfiles,
39                                      project.gauge_filename,
40                                      production_dirs,
41                                      report = False,
42                                      reportname = '',
43                                      plot_quantity = ['stage', 'speed'],
44                                      surface = False,
45                                      time_min = None,
46                                      time_max = None,
47                                      title_on = False,
48                                      verbose = True)
Note: See TracBrowser for help on using the repository browser.