Ignore:
Timestamp:
Mar 21, 2007, 2:17:41 PM (18 years ago)
Author:
nick
Message:

update to dampier

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/dampier_2006/get_timeseries.py

    r4282 r4308  
    2323   
    2424# sww file created from URS boundary
    25 production_dirs = {#'20070206_005106_run': '10000 yr wave height'}
    26                    '20070223_031527_run': 'URS boundary condition'}
     25production_dirs = {'boundaries': '10000 yr wave height'
     26                   #'20070307_224017_run': 'URS boundary normal'
     27                   }
    2728
     29   
    2830is_parallel = False
    2931if is_parallel == True:
     
    3335swwfiles = {}
    3436
     37#J:\inundation\data\western_australia\onslow_tsunami_scenario_2006\anuga\boundaries SU-AU_high_res1.sww
     38
     39#boundary_dir_filename= os.path.join(home,data,state,'onslow_tsunami_scenario_2006',an,bo,'SU-AU_high_res1.sww')
     40#gauges_dir_name_simple =
    3541if is_parallel == True:
     42    print 'buggur'
     43   
    3644
    37     for i in range(nodes):
    38         print 'Sending node %d of %d' %(i,nodes)
    39         swwfiles = {}
    40         #reportname = report_name + '_%s' %(i)
    41         for label_id in production_dirs.keys():
    42             if label_id == 'boundaries':
    43                 file_loc = project.boundaries_in_dir
    44                 swwfile = file_loc + sep + '20070223_031527_run' + sep + 'dampier.sww'
    45                 print 'hi'
    46             #else:
    47             #file_loc = project.output_dir + label_id + sep
    48             #sww_extra = '_P%s_%s' %(i,nodes)
    49             #swwfile = file_loc + project.scenario_name + sww_extra + '.sww'
    50             swwfiles[swwfile] = label_id
    51 
    52         texname, elev_output = sww2timeseries(swwfiles,
    53                                               project.gauges_dir_name_test,
    54                                               production_dirs,
    55                                               report = False,
    56                                               #reportname = reportname,
    57                                               plot_quantity = ['stage', 'momentum'],
    58                                               generate_fig = False,
    59                                               surface = False,
    60                                               time_min = None,
    61                                               time_max = None,
    62                                               title_on = False,
    63                                               verbose = True)
    64            
    65             #latex_output.append(texname)
    66    
    6745else:
    6846   
     
    7048       
    7149        if label_id == 'boundaries':
     50                print 'boundaries'
    7251                file_loc = project.boundaries_in_dir
    73                 swwfile = file_loc + sep + '20070223_031527_run' + sep + 'dampier.sww'
     52                swwfile = project.boundaries_dir_name3 + '.sww'
     53#                swwfile = boundary_dir_filename
    7454        else:
    7555            file_loc = project.output_dir + label_id + sep
     
    7858       
    7959    texname, elev_output = sww2timeseries(swwfiles,
    80                                           project.gauges_dir_name_test,
     60                                          project.gauges_dir_name_simple,
    8161                                          production_dirs,
    8262                                          report = False,
    8363                                          #reportname = report_name,
    84                                           plot_quantity = ['stage', 'momentum'],
     64#                                          plot_quantity = ['stage', 'momentum'],
     65                                          plot_quantity = ['stage', 'xmomentum', 'ymomentum'],
    8566                                          surface = False,
    8667                                          time_min = None,
     
    8869                                          title_on = False,
    8970                                          verbose = True)
     71                                         
Note: See TracChangeset for help on using the changeset viewer.