source: anuga_work/production/hobart_2009/combine_gauges.py @ 7355

Last change on this file since 7355 was 7355, checked in by kristy, 15 years ago
File size: 1.5 KB
Line 
1# Standard modules
2from os import sep
3import os, glob
4from os.path import join
5
6# Related major packages
7from anuga.utilities.csv_tools import merge_csv_key_values
8
9# Application specific imports
10from setup_model import project  # Definition of file names and polygons
11import csv
12directory = project.output_folder
13
14time_dir1 = 'Event1_MSL'
15time_dir2 = 'Event2_MSL'
16time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
17time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
18time_dir5 = '20090619_165107_run_final_0_58260_None_kvanputt'
19time_dir6 = '20090619_165136_run_final_0.8_58260_None_kvanputt'
20
21time_dirs = [time_dir1, time_dir2] #, time_dir3, time_dir4] 
22
23data = 'stage'
24output_filename = 'gauges_combined_' + data + '.csv' 
25for time_dir in time_dirs:
26    print time_dir
27    folder = join(directory, time_dir)
28    output_file = join(folder, output_filename) 
29    #Looking for 'gauge_<alphabetic>.csv' filename
30    input_csv = glob.glob(folder + sep + 'gauge_[A-Za-z]*.csv')
31   
32    csv_files = []
33    for infile in input_csv:
34        title = os.path.basename(infile)[:-4]
35        title = title.replace('gauge_', '')
36        if title in csv_files:
37            print 'Title (%s) already in use' % title
38            sys.exit(10)
39        csv_files.append([infile, title])
40
41   # print 'csv_files=%s' % str(csv_files)
42    merge_csv_key_values(csv_files, output_file,
43                         key_col='hours', data_col=data)
Note: See TracBrowser for help on using the repository browser.