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 | """ |
---|
2 | Generate images of "gauges" for production run |
---|
3 | |
---|
4 | Inputs: |
---|
5 | |
---|
6 | production dirs: dictionary of production directories with a |
---|
7 | association to that simulation run, eg high tide, |
---|
8 | low tide and MSL. |
---|
9 | |
---|
10 | Outputs: |
---|
11 | |
---|
12 | * figures used for report stored in the report_figure directory |
---|
13 | |
---|
14 | """ |
---|
15 | |
---|
16 | from os import getcwd, sep, altsep, mkdir, access, F_OK |
---|
17 | import project |
---|
18 | from anuga.abstract_2d_finite_volumes.util import sww2timeseries |
---|
19 | |
---|
20 | # Create report directory |
---|
21 | reportdir = getcwd()+sep+'report'+sep |
---|
22 | if access(reportdir,F_OK) == 0: |
---|
23 | mkdir (reportdir) |
---|
24 | |
---|
25 | # 2500 interior resolution for both events |
---|
26 | # use both to compare the events |
---|
27 | production_dirs = {'20061006_062319': 'Mw 8.5', |
---|
28 | '20061008_234702': 'Mw 8.7'} |
---|
29 | |
---|
30 | # Generate figures |
---|
31 | swwfiles = {} |
---|
32 | for 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 | |
---|
38 | texname, 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.