source: anuga_work/development/Hinwood_2008/slope.py @ 5413

Last change on this file since 5413 was 5413, checked in by duncan, 16 years ago

Current Hinwood scenario

File size: 1.4 KB
Line 
1
2"""
3Plot up files from the Hinwood project.
4"""
5from os import sep
6import project
7#from scipy import arange
8from Numeric import arange
9
10from anuga.fit_interpolate.interpolate import interpolate_sww2csv
11
12
13#-------------------------------------------------------------
14if __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
Note: See TracBrowser for help on using the repository browser.