source: production/onslow_2006/compare_timeseries_data.py @ 3514

Last change on this file since 3514 was 3514, checked in by duncan, 18 years ago

Hi all,
I'm doing a change in the anuga structure, moving the code to

\anuga_core\source\anuga

After you have done an svn update, the PYTHONPATH has to be changed to;
PYTHONPATH = anuga_core/source/

This is part of changes required to make installation of anuga quicker and reducing the size of our sandpits.

If any imports are broken, try fixing them. With adding anuga. to them for example. If this seems to have really broken things, email/phone me.

Cheers
Duncan

File size: 1.6 KB
Line 
1"""
2Generate latex output comparing time series of given sww files
3
4Inputs:
5
6output dirs: dictionary of output directories with a
7             association to simulation, eg 50m boundary, 100m boundary
8                   
9Outputs:
10
11* latex output to generate figures for report
12* figures used for report stored in the report_figure directory
13NOTE, this directory will need to be committed, as well as
14the latex files.
15
16June 2006                   
17"""
18
19from os import sep
20import project
21from anuga.pyvolution.util import sww2timeseries, get_gauges_from_file
22
23# User defined inputs
24
25# for comparison between DTED and WA DLI data
26production_dirs = {'20060515_001733': 'DTED',
27                   '20060704_063234': 'WA DLI'}
28
29plot_quantity = ['stage', 'speed']
30
31# Create sections and graphs for each designated production directory
32latex_output = []
33swwfiles = {}
34for label_id in production_dirs.keys():
35   
36    file_loc = project.outputdir + label_id + sep
37    swwfile = file_loc + project.basename + '.sww'
38    swwfiles[swwfile] = label_id
39
40texname, vec = sww2timeseries(swwfiles,
41                              project.gauge_filename,
42                              production_dirs,
43                              report = True,
44                              reportname = 'compare_output_datasets',
45                              surface = False,
46                              plot_quantity = plot_quantity,
47                              time_min = None,
48                              time_max = None,
49                              title_on = False,
50                              verbose = True)
51
52latex_output.append(texname)
Note: See TracBrowser for help on using the repository browser.