Changeset 5569 for anuga_work/production/perth/build_boundary.py

Timestamp:

Jul 25, 2008, 8:53:51 AM (15 years ago)

Author:

kristy

Message:

updated scripts for polyline boundary

File:

• anuga_work/production/perth/build_boundary.py (modified) (2 diffs)

anuga_work/production/perth/build_boundary.py

@@ -34 +34 @@
 urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()]
 
-# Reduce this list slightly
-urs_filenames = urs_filenames[:30]
-for name in urs_filenames:
-    print name
+# AS per David Burbidge email on friday 4th July the mag 9.3 event
+# has 1m worth of slip on each sub fault therefore mutliple each unit
+# soucre by the slip (10.4544) and sum the 44 time series together to
+# get the time series for this event at the points on your boundary. 
 
 weights=10.4544*ones(len(urs_filenames),Float)
 
-base_name=project.scenario_name
-order_filename=os.path.join(project.boundaries_dir, 'boundary_ordering.txt')
+scenario_name=project.scenario_name
+order_filename=os.path.join(project.boundaries_dir, 'thinned_bound_order.txt')
 
 print 'reading', order_filename
 # Create ordered sts file
 print 'creating sts file'
-urs2sts(urs_filenames,basename_out=base_name,
+
+urs2sts(urs_filenames,basename_out=scenario_name,
         ordering_filename=order_filename,
         weights=weights,
@@ -54 +55 @@
 
 # Read in boundary from ordered sts file
-urs_boundary=create_sts_boundary(base_name)
+urs_boundary_polygon=create_sts_boundary(scenario_name)
 
+
+
+
+
+
+