Changeset 6784


Ignore:
Timestamp:
Apr 14, 2009, 10:15:49 AM (15 years ago)
Author:
kristy
Message:

Reverted back to previous version to pass through 3 arguments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/melbourne/build_urs_boundary.py

    r6760 r6784  
    9292# @param output_dir Directory to write STS data to.
    9393# @note 'event_file' is produced by EventSelection.
    94 def build_urs_boundary(event_file, output_dir):
     94def build_urs_boundary(event_file, output_dir,urs_file):
    9595    '''Build a boundary STS file from a set of MUX files.'''
    96 
     96    print 'reading', urs_file
    9797    # if we are using an EventSelection multi-mux file
    9898    if project.multi_mux:
     
    137137        urs2sts(mux_filenames,
    138138                basename_out=output_dir,
    139                 ordering_filename=project.urs_order,
     139                ordering_filename=urs_file,
    140140                weights=mux_weights,
    141141                central_meridian=project.central_meridian,
     
    149149        weight_factor = 1.0
    150150        mux_weights = weight_factor*num.ones(len(mux_filenames), num.Float)
    151            
    152         order_filename = project.urs_order
    153151
    154         print 'reading', order_filename
     152
    155153        # Create ordered sts file
    156154        urs2sts(mux_filenames,
    157155                basename_out=output_dir,
    158                 ordering_filename=order_filename,
     156                ordering_filename=urs_file,
    159157                weights=mux_weights,
    160158                central_meridian=project.central_meridian,
     
    163161
    164162    # report on progress so far
    165     sts_file = os.path.join(project.event_folder, project.scenario_name)
    166     quantities, elevation, time = get_sts_gauge_data(sts_file, verbose=False)
     163##    sts_file = os.path.join(project.event_folder, project.scenario_name)
     164    quantities, elevation, time = get_sts_gauge_data(output_dir, verbose=False)
    167165    print len(elevation), len(quantities['stage'][0,:])
    168166
Note: See TracChangeset for help on using the changeset viewer.