source: anuga_work/production/dampier_2006/get_timeseries.py @ 4308

"""
Generate images of "gauges" for production run

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

# Create report directory
reportdir = getcwd()+sep+'report'+sep
if access(reportdir,F_OK) == 0:
    mkdir (reportdir)
    
# sww file created from URS boundary
production_dirs = {'boundaries': '10000 yr wave height'
                   #'20070307_224017_run': 'URS boundary normal'
                   }

    
is_parallel = False
if is_parallel == True:
    nodes = 2
    
# Generate figures
swwfiles = {}

#J:\inundation\data\western_australia\onslow_tsunami_scenario_2006\anuga\boundaries SU-AU_high_res1.sww

#boundary_dir_filename= os.path.join(home,data,state,'onslow_tsunami_scenario_2006',an,bo,'SU-AU_high_res1.sww')
#gauges_dir_name_simple = 
if is_parallel == True:
    print 'buggur'
   

else:
    
    for label_id in production_dirs.keys():
        
        if label_id == 'boundaries':
                print 'boundaries'
                file_loc = project.boundaries_in_dir
                swwfile = project.boundaries_dir_name3 + '.sww'
#                swwfile = boundary_dir_filename
        else:
            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_simple,
                                          production_dirs,
                                          report = False,
                                          #reportname = report_name,
#                                          plot_quantity = ['stage', 'momentum'],
                                          plot_quantity = ['stage', 'xmomentum', 'ymomentum'],
                                          surface = False,
                                          time_min = None,
                                          time_max = None,
                                          title_on = False,
                                          verbose = True)