Changeset 5415 for anuga_work/production/perth

Timestamp:

Jun 24, 2008, 10:24:21 AM (17 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/perth

Files:

• GetTimeseries.py (modified) (1 diff)
• export_results.py (modified) (1 diff)
• export_results_all.py (modified) (4 diffs)
• plot_timeseries.py (added)
• project.py (modified) (2 diffs)
• run_perth.py (modified) (1 diff)

anuga_work/production/perth/GetTimeseries.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 images of "gauges" for production run
 
-Note, this script will only work if pylab is installed on the platform
+Inputs:
+
+production dirs: dictionary of production directories with a
+                 association to that simulation run, eg high tide,
+                 low tide and MSL.
+                    
+Outputs:
+
+* figures used for report stored in the report_figure directory
 
 """
+from os import getcwd, sep, altsep, mkdir, access, F_OK
+import project
+from anuga.abstract_2d_finite_volumes.util import sww2timeseries
 
-from os import sep
-import project
-from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
+production_dirs = {'20080619_115643_run_trial_0.6_exmouth_kvanputt': ''}
 
-sww2csv_gauges('slide'+sep+'slidesource.sww',
-                project.gauge_filename,
-                quantities = ['stage','speed','depth','elevation'],
-                verbose=True)
-                
-sww2csv_gauges('fixed_wave'+sep+'fixed_wavesource.sww',
-               project.gauge_filename,
-               quantities = ['stage', 'speed','depth','elevation'],
-               verbose=True)
+is_parallel = False
 
-try: 
-    import pylab
-    csv2timeseries_graphs(directories_dic={'slide'+sep:['Slide',0, 0],
-                                           'fixed_wave'+sep:['Fixed Wave',0,0]},
-                            output_dir='fixed_wave'+sep,
-                            base_name='gauge_',
-                            plot_numbers='',
-                            quantities=['stage','speed','depth'],
-                            extra_plot_name='',
-                            assess_all_csv_files=True,                            
-                            create_latex=False,
-                            verbose=True)
-except ImportError:
-    #ANUGA does not rely on pylab to work 
-    print 'must have pylab installed to generate plots'
+if is_parallel is True:
+    no_nodes = 4
     
+if is_parallel is False:
+    
+    swwfiles = {}    
+    for label_id in production_dirs.keys():
+        
+        file_loc = project.output_dir + label_id + sep
+        swwfile = file_loc + project.scenario_name + '.sww'
+        swwfiles[swwfile] = label_id
+            
+        texname, elev_output = sww2timeseries(swwfiles,
+                                              project.gauges_dir_name,
+                                              production_dirs,
+                                              report = False,
+                                              plot_quantity = ['stage', 'speed', 'depth'],
+                                              generate_fig = False,
+                                              surface = False,
+                                              time_min = None,
+                                              time_max = None,
+                                              time_unit = 'mins',
+                                              title_on = False,
+                                              use_cache = True,
+                                              verbose = True)
 
+if is_parallel is True:                                      
+
+    for count in range(no_nodes):
+        swwfiles = {} 
+        for label_id in production_dirs.keys():
+
+            file_loc = project.output_dir + label_id + sep
+            extra = '_P%i_%i' %(count,no_nodes)
+            swwfile = file_loc + project.scenario_name + extra + '.sww'
+            swwfiles[swwfile] = label_id
+
+            texname, elev_output = sww2timeseries(swwfiles,
+                                                  project.gauges_dir_name,
+                                                  production_dirs,
+                                                  report = False,
+                                                  plot_quantity = ['stage', 'speed'],
+                                                  generate_fig = False,
+                                                  surface = False,
+                                                  time_min = None,
+                                                  time_max = None,
+                                                  time_unit = 'mins',
+                                                  title_on = False,
+                                                  use_cache = True,
+                                                  verbose = True)
+                                      

anuga_work/production/perth/export_results.py

@@ -6 +6 @@
 
 #time_dir = '20080526_104946_run_final_0.6_test_kvanputt'
-time_dir = '20080530_170833_run_final_0.6_exmouth_kvanputt'
-
+#time_dir = '20080530_170833_run_final_0.6_exmouth_kvanputt'
+time_dir = '20080619_115643_run_trial_0.6_exmouth_kvanputt'
 
 cellsize = 25

anuga_work/production/perth/export_results_all.py

@@ -7 +7 @@
 #time_dir = '20080526_104946_run_final_0.6_test_kvanputt'
 #time_dir = '20080530_170833_run_final_0.6_exmouth_kvanputt'
-time_dir = '20080612_102742_run_final_0.6_exmouth_kvanputt'
+time_dir = '20080619_115643_run_trial_0.6_exmouth_kvanputt'
 
 
 
-cellsize = 25
+cellsize = 10
 #cellsize = 150
-#timestep = 0
+timestep = 0
 directory = project.output_dir
 name = directory+time_dir+sep+project.scenario_name
@@ -53 +53 @@
         northing_max = project.ymaxRockingham
 
-    var = [2,3,4] # depth and speed
+   # var = [2,3,4] # depth and speed
 #var = [2] # depth
-#var = [4]
+    var = [0,4]
 
     for which_var in var:
@@ -86 +86 @@
                 sww2dem(namei, basename_out = outnamei,
                             quantity = quantityname,
-                            #timestep = timestep,
+                            timestep = timestep,
                             cellsize = cellsize,      
                             easting_min = project_grad.e_min_area,
@@ -99 +99 @@
             sww2dem(name, basename_out = outname,
                         quantity = quantityname,
-                        #timestep = timestep,
+                        timestep = timestep,
                         cellsize = cellsize,      
                         easting_min = easting_min,

anuga_work/production/perth/project.py

@@ -37 +37 @@
 alpha = 0.1
 friction=0.01
-starttime=10000
+starttime=0
 finaltime=80000
-export_cellsize=50
-setup='final'
+export_cellsize=25
+setup='trial'
 source='exmouth'
 
@@ -143 +143 @@
 #gauges
 gauge_name = 'perth.csv'
-gauges_dir = anuga_dir+'gauges'+sep
+
+gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
+beach_gauges = gauges_dir + 'beach_gauges.csv'
 gauges_dir_name = gauges_dir + gauge_name
 

anuga_work/production/perth/run_perth.py

@@ -193 +193 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']
-                       ,skip_initial_step = True): 
+    for t in domain.evolve(yieldstep = project.yieldstep, finaltime = kwargs['finaltime']
+                       ,skip_initial_step = False): 
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')