[6770] | 1 | """ |
---|
| 2 | Generate time series of nominated "gauges". |
---|
| 3 | Input: sww file from run_perth.py |
---|
| 4 | gauges project.gauge_filename |
---|
| 5 | Output: csv files stage,speed,depth,elevation over time |
---|
| 6 | Stored in the 'outputs_dir' folder for respective .sww file |
---|
| 7 | |
---|
| 8 | Note: |
---|
| 9 | Can run different sww files at the same time |
---|
| 10 | make sure if there is a second sww file that it is placed into a folder called sww2 |
---|
| 11 | Can run different gauges at the same time - ie testing boundary index point |
---|
| 12 | """ |
---|
| 13 | |
---|
| 14 | from os import sep, rename, listdir, system |
---|
| 15 | from os.path import join |
---|
| 16 | from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs |
---|
| 17 | import project |
---|
| 18 | |
---|
| 19 | directory = project.output_folder |
---|
| 20 | |
---|
| 21 | time_dir1 = '20090324_135427_run_final_0.6_27255_Bt_kvanputt' |
---|
| 22 | |
---|
| 23 | time_dirs = [time_dir1] #1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] |
---|
| 24 | |
---|
| 25 | print 'time directories', time_dirs |
---|
| 26 | |
---|
| 27 | for time_dir in time_dirs: |
---|
| 28 | name = join(directory,time_dir,project.scenario_name) |
---|
| 29 | gauge = project.gauges |
---|
| 30 | sww2csv_gauges(name+'.sww',gauge, |
---|
| 31 | quantities = ['stage','speed','depth','elevation'], |
---|
| 32 | verbose=True) |
---|
| 33 | |
---|
| 34 | |
---|