Changeset 6017

Timestamp:

Nov 28, 2008, 11:32:55 AM (16 years ago)

Author:

kristy

Message:

restructured get_timeseries.py to fit style of all other scripts

File:

• anuga_work/production/busselton/get_timeseries.py (modified) (1 diff)

anuga_work/production/busselton/get_timeseries.py

@@ -1 +1 @@
 """
-Generate time series of nominated "gauges" read from project.gauge_filename. This 
-is done by first running sww2csv_gauges on two different directories to make 
-'csv' files. Then running csv2timeseries_graphs detailing the two directories 
-containing the csv file and produces one set of graphs in the 'output_dir' containing
-the details at the gauges for both these sww files.
+Generate time series of nominated "gauges".
+Input: sww file from run_perth.py
+       gauges project.gauge_filename
+Output: csv files stage,speed,depth,elevation over time
+Stored in the 'outputs_dir' folder for respective .sww file
 
-Note, this script will only work if pylab is installed on the platform
-
+Note:
+Can run different sww files at the same time
+make sure if there is a second sww file that it is placed into a folder called sww2
+Can run different gauges at the same time - ie testing boundary index point
 """
 
-from os import sep
+from os import sep, rename, listdir, system
 from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
-
 import project
 
+directory = project.output_dir
 
-#timestamp='20080724_121200_run_trial_0.6_polyline_alpha0.1_kvanputt'
-#timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
-timestamp1='20080815_103818_run_final_0_polyline_newExtent_kvanputt'
-timestamp2='20080815_103708_run_final_0.6_polyline_newExtent_kvanputt'
+##time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt'
+##time_dir2 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt'
+time_dir3 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'
+##time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'
+##time_dir5 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt'
+##time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'
 
-timestamps = [timestamp1, timestamp2]
-for timestamp in timestamps:
+time_dirs = [time_dir3] #1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
-    filename1=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
-    filename2=project.output_dir+timestamp+sep+'busselton_time_38340'+sep+project.scenario_name+'_time_38340_0.sww'
+for time_dir in time_dirs:
+    
+    name1 = directory+time_dir+sep+project.scenario_name
+    #name2 = directory+time_dir+sep+project.scenario_name+'_time_39600_0'
+    #name3 = directory+time_dir+sep+project.scenario_name+'_time_79200_0'
 
-    filenames = [filename1, filename2]
-    for filename in filenames:
+    names = [name1] #, name2, name3]
 
-        sww2csv_gauges(filename,
-                        project.gauges_dir_name,
-                        #project.gauges_dir_name2,
+    for i, name in enumerate(names):
+        gauge = project.gauges_dir_name
+                  
+        sww2csv_gauges(name+'.sww',gauge,
                         quantities = ['stage','speed','depth','elevation'],
                         verbose=True)
-                    
+##        for filename in listdir(time_dir):
+##            total_filename = 'All_' + filename
+##            if i == 0:
+##                if filename.startswith('gauge_'):
+##                    rename(filename,total_filename)
+##            else:
+##                system('cat %s %s > %s' %(total_filename,filename,total_filename))
+##
+       
 
-    
-