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

Timestamp:

Aug 7, 2008, 8:54:24 AM (16 years ago)

Author:

sexton

Message:

updated build_boundary based on 44 unit sources as opposed to all mux2 files that were sitting in the mux directory (i.e. 49)

File:

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

anuga_work/production/perth/build_boundary.py

@@ -10 +10 @@
 import project
 
-
 ######################################
 # Create urs boundary from mux2files #
 ######################################
 
-#print project.muxhome
 dir=os.path.join(project.muxhome,'mux')
-#print dir 
-prefix=os.path.join(dir,'Java-00')
-suffix='-z.grd'
 
-urs_filenames={}
-for filename in os.listdir(dir):
-    basename = filename[:-7]
+# taken from David's event list
+urs_filenames = [
+    os.path.join(dir,'Java-0000-z.grd'),
+    os.path.join(dir,'Java-0001-z.grd'),
+    os.path.join(dir,'Java-0002-z.grd'),
+    os.path.join(dir,'Java-0003-z.grd'),
+    os.path.join(dir,'Java-0005-z.grd'), 
+    os.path.join(dir,'Java-0006-z.grd'), 
+    os.path.join(dir,'Java-0007-z.grd'), 
+    os.path.join(dir,'Java-0008-z.grd'), 
+    os.path.join(dir,'Java-0010-z.grd'), 
+    os.path.join(dir,'Java-0011-z.grd'), 
+    os.path.join(dir,'Java-0012-z.grd'), 
+    os.path.join(dir,'Java-0013-z.grd'),  
+    os.path.join(dir,'Java-0015-z.grd'),  
+    os.path.join(dir,'Java-0016-z.grd'),  
+    os.path.join(dir,'Java-0017-z.grd'),  
+    os.path.join(dir,'Java-0018-z.grd'),  
+    os.path.join(dir,'Java-0020-z.grd'),  
+    os.path.join(dir,'Java-0021-z.grd'),  
+    os.path.join(dir,'Java-0022-z.grd'),  
+    os.path.join(dir,'Java-0023-z.grd'),  
+    os.path.join(dir,'Java-0025-z.grd'),  
+    os.path.join(dir,'Java-0026-z.grd'),  
+    os.path.join(dir,'Java-0027-z.grd'),  
+    os.path.join(dir,'Java-0028-z.grd'),  
+    os.path.join(dir,'Java-0029-z.grd'),  
+    os.path.join(dir,'Java-0030-z.grd'),  
+    os.path.join(dir,'Java-0031-z.grd'),  
+    os.path.join(dir,'Java-0032-z.grd'),  
+    os.path.join(dir,'Java-0033-z.grd'),  
+    os.path.join(dir,'Java-0034-z.grd'),  
+    os.path.join(dir,'Java-0035-z.grd'),  
+    os.path.join(dir,'Java-0036-z.grd'),  
+    os.path.join(dir,'Java-0037-z.grd'),  
+    os.path.join(dir,'Java-0038-z.grd'),  
+    os.path.join(dir,'Java-0039-z.grd'),  
+    os.path.join(dir,'Java-0040-z.grd'),  
+    os.path.join(dir,'Java-0041-z.grd'),  
+    os.path.join(dir,'Java-0042-z.grd'),  
+    os.path.join(dir,'Java-0043-z.grd'),  
+    os.path.join(dir,'Java-0044-z.grd'),  
+    os.path.join(dir,'Java-0045-z.grd'),  
+    os.path.join(dir,'Java-0046-z.grd'),  
+    os.path.join(dir,'Java-0047-z.grd'),
+    os.path.join(dir,'Java-0048-z.grd')]  
+ 
+print 'number of sources', len(urs_filenames)
 
-    urs_filenames[basename] = True
+#from os import sys
+#sys.exit()
+# for all mux files
+#urs_filenames={}
+##for filename in os.listdir(dir)
+##    basename = filename[:-7]
+##
+##    urs_filenames[basename] = True
+    
+# Make list of basenames  
+#urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()]
 
-# Make list of basenames  
-urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()]
-
-# AS per David Burbidge email on friday 4th July the mag 9.3 event
+# 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
+# source 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)
+weight_factor = 10.4544
+weights=weight_factor*ones(len(urs_filenames),Float)
 
-scenario_name=project.scenario_name
+scenario_name='perth' #project.scenario_name
 order_filename=os.path.join(project.boundaries_dir, 'thinned_bound_order.txt')
 
@@ -50 +98 @@
         mean_stage=project.tide,
         verbose=True)
-
-
-
-
-