Changeset 578 for inundation/ga/storm_surge/pmesh/load_mesh

Timestamp:

Nov 17, 2004, 3:28:09 PM (20 years ago)

Author:

duncan

Message:

added function to make boundary file for alpha shape

Location:

inundation/ga/storm_surge/pmesh/load_mesh

Files:

• loadASCII.py (modified) (3 diffs)
• loadASCIIHarness.py (modified) (1 diff)

inundation/ga/storm_surge/pmesh/load_mesh/loadASCII.py

@@ -411 +411 @@
         raise IOError, msg 
 
+## used by alpha shapes
+    
+def export_boundary_file( file_name, points, title, delimiter = ','):
+    """
+    export a file, ofile, with the format
+    
+    First line: Title variable
+    Following lines:  [point index][delimiter][point index] 
+
+    file_name - the name of the new file
+    points - List of point index pairs [[p1, p2],[p3, p4]..] 
+    title - info to write in the first line
+    """
+    
+    fd = open(file_name,'w')
+    
+    fd.write(title+"\n")
+    #[point index][delimiter][point index] 
+    for point in points:
+        fd.write( str(point[0]) + delimiter
+                  + str(point[1]) + "\n")
+    fd.close()
+
+
 ###
 #  LOADING XYA FILES
@@ -423 +447 @@
     try:
         fd = open(ofile)
+        #print "ofile",ofile
         xya_dic = read_xya_file(fd, delimiter)
         fd.close()
@@ -510 +535 @@
                   + attlist + "\n")
     fd.close()
-
+     
 if __name__ == "__main__":
     m = import_mesh("tee.txt")

inundation/ga/storm_surge/pmesh/load_mesh/loadASCIIHarness.py

@@ -4 +4 @@
 #FIXME (Ole): I think we need a unit test of write_ASCII_trianglulation 
 #             in order to quash bug in merimbula data
-
+#  - to quash the bug need a raised error in pyvolution for reading xya files
+# when the title attribute information is not consistent with the data in the
+# file - DSG
 
 import sys