[5413] | 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 | |
---|