Ignore:
Timestamp:
Jul 25, 2008, 8:53:51 AM (14 years ago)
Author:
kristy
Message:

updated scripts for polyline boundary

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/build_boundary.py

    r5551 r5569  
    3434urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()]
    3535
    36 # Reduce this list slightly
    37 urs_filenames = urs_filenames[:30]
    38 for name in urs_filenames:
    39     print name
     36# AS per David Burbidge email on friday 4th July the mag 9.3 event
     37# has 1m worth of slip on each sub fault therefore mutliple each unit
     38# soucre by the slip (10.4544) and sum the 44 time series together to
     39# get the time series for this event at the points on your boundary.
    4040
    4141weights=10.4544*ones(len(urs_filenames),Float)
    4242
    43 base_name=project.scenario_name
    44 order_filename=os.path.join(project.boundaries_dir, 'boundary_ordering.txt')
     43scenario_name=project.scenario_name
     44order_filename=os.path.join(project.boundaries_dir, 'thinned_bound_order.txt')
    4545
    4646print 'reading', order_filename
    4747# Create ordered sts file
    4848print 'creating sts file'
    49 urs2sts(urs_filenames,basename_out=base_name,
     49
     50urs2sts(urs_filenames,basename_out=scenario_name,
    5051        ordering_filename=order_filename,
    5152        weights=weights,
     
    5455
    5556# Read in boundary from ordered sts file
    56 urs_boundary=create_sts_boundary(base_name)
     57urs_boundary_polygon=create_sts_boundary(scenario_name)
    5758
     59
     60
     61
     62
     63
     64
Note: See TracChangeset for help on using the changeset viewer.