1 | |
---|
2 | """ |
---|
3 | Plot up files from the Hinwood project. |
---|
4 | """ |
---|
5 | from os import sep |
---|
6 | import project |
---|
7 | #from scipy import arange |
---|
8 | from Numeric import arange |
---|
9 | |
---|
10 | from anuga.fit_interpolate.interpolate import interpolate_sww2csv |
---|
11 | |
---|
12 | |
---|
13 | #------------------------------------------------------------- |
---|
14 | if __name__ == "__main__": |
---|
15 | """ |
---|
16 | """ |
---|
17 | from scenarios import scenarios |
---|
18 | |
---|
19 | outputdir_tag = "_good_tri_area_0.01_A" |
---|
20 | scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!! |
---|
21 | for run_data in scenarios: |
---|
22 | point_x = arange(run_data['start_slope_x'], |
---|
23 | run_data['finish_slope_x'], |
---|
24 | 0.01).tolist() |
---|
25 | flume_y_middle = 0.5 |
---|
26 | points = [] |
---|
27 | for gauge_x in point_x: |
---|
28 | points.append([gauge_x, flume_y_middle]) |
---|
29 | id = run_data['scenario_id'] |
---|
30 | |
---|
31 | basename = 'zz_' + run_data['scenario_id'] |
---|
32 | outputdir_name = id + outputdir_tag |
---|
33 | pro_instance = project.Project(['data','flumes','Hinwood_2008'], |
---|
34 | outputdir_name=outputdir_name) |
---|
35 | end = id + ".csv" |
---|
36 | interpolate_sww2csv(pro_instance.outputdir + basename +".sww", |
---|
37 | points, |
---|
38 | pro_instance.outputdir + "depth_" + end, |
---|
39 | pro_instance.outputdir + "velocity_x_" + end, |
---|
40 | pro_instance.outputdir + "velocity_y_" + end) |
---|
41 | |
---|
42 | |
---|
43 | |
---|