Changeset 5569 for anuga_work/production/perth/build_boundary.py
- Timestamp:
- Jul 25, 2008, 8:53:51 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/perth/build_boundary.py
r5551 r5569 34 34 urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()] 35 35 36 # Reduce this list slightly37 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. 40 40 41 41 weights=10.4544*ones(len(urs_filenames),Float) 42 42 43 base_name=project.scenario_name44 order_filename=os.path.join(project.boundaries_dir, ' boundary_ordering.txt')43 scenario_name=project.scenario_name 44 order_filename=os.path.join(project.boundaries_dir, 'thinned_bound_order.txt') 45 45 46 46 print 'reading', order_filename 47 47 # Create ordered sts file 48 48 print 'creating sts file' 49 urs2sts(urs_filenames,basename_out=base_name, 49 50 urs2sts(urs_filenames,basename_out=scenario_name, 50 51 ordering_filename=order_filename, 51 52 weights=weights, … … 54 55 55 56 # Read in boundary from ordered sts file 56 urs_boundary =create_sts_boundary(base_name)57 urs_boundary_polygon=create_sts_boundary(scenario_name) 57 58 59 60 61 62 63 64
Note: See TracChangeset
for help on using the changeset viewer.