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

Timestamp:

Jul 22, 2008, 1:39:46 PM (16 years ago)

Author:

ole

Message:

Work towards using the STS boundary with the Perth production script

File:

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

anuga_work/production/perth/build_boundary.py

@@ -21 +21 @@
 suffix='-z.grd'
 
-urs_filenames=[]
-for i in range(49):
-    if i<10:
-        urs_filenames.append(prefix+'0'+str(i)+suffix)
-    else:
-        urs_filenames.append(prefix+str(i)+suffix)
-print urs_filenames
+
+
+
+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()]
+
+# Reduce this list slightly
+urs_filenames = urs_filenames[:30]
+for name in urs_filenames:
+    print name
 
 weights=10.4544*ones(len(urs_filenames),Float)
 
 base_name=project.scenario_name
-order_filename=os.path.join(project.boundaries_dir,'boundary_ordering.txt')
+order_filename=os.path.join(project.boundaries_dir, 'boundary_ordering.txt')
 
+print 'reading', order_filename
 # Create ordered sts file
 print 'creating sts file'
@@ -40 +51 @@
         weights=weights,
         mean_stage=project.tide,
-        verbose=False)
+        verbose=True)
 
 # Read in boundary from ordered sts file