Changeset 1395 for inundation/ga/storm_surge/pmesh/mesh.py

Timestamp:

May 16, 2005, 11:22:40 AM (19 years ago)

Author:

duncan

Message:

load_mesh refactoring

File:

• inundation/ga/storm_surge/pmesh/mesh.py (modified) (5 diffs)

inundation/ga/storm_surge/pmesh/mesh.py

@@ -1731 +1731 @@
                             self.regions)
         
-        # exportASCIImeshfile   
+        # exportASCIImeshfile   - this function is used
     def export_mesh_file(self,ofile):
         """
@@ -1773 +1773 @@
             raise RuntimeError
         numVertAttrib = str(len(Vertices[0].attributes))
-        title = numVert + " " + numVertAttrib + " # <vertex #> <x> <y> [attributes]"
+        #title = numVert + " " + numVertAttrib + " # <vertex #> <x> <y> [attributes]"
 
         #Convert the Vertices to pointlist and pointattributelist
@@ -1788 +1788 @@
         xya_dict['geo_reference'] = self.geo_reference
 
-        load_mesh.loadASCII.export_xya_file(ofile, xya_dict, title, delimiter = " ")
+        load_mesh.loadASCII.export_points_file(ofile, xya_dict)
 
         
@@ -2732 +2732 @@
 def importMeshFromFile(ofile):
     """returns a mesh object, made from a .xya/.pts or .tsh/.msh file
-    Often raises SyntaxError, IOError
+    Often raises IOError,RuntimeError
     """
     newmesh = None
-    if ofile[-4:]== ".xya":
-        #print "loading " + ofile
-        try:
-            dict = load_mesh.loadASCII.load_points_file(ofile, delimiter = ',')
-        except SyntaxError:
-            dict = load_mesh.loadASCII.load_points_file(ofile, delimiter = ' ')
-        #print "dict",dict   outline_     
-        dict['segmentlist'] = []
-        dict['segmenttaglist'] = []
-        dict['regionlist'] = []
-        dict['regionattributelist'] = []
-        dict['regionmaxarealist'] = []
-        dict['holelist'] = []
-        newmesh= Mesh()
-        newmesh.setMesh(dict) #FIXME use IOOutline2Mesh
-        counter = newmesh.removeDuplicatedUserVertices()
-        if (counter >0):
-            print "%i duplicate vertices removed from dataset" % (counter)
-    elif ofile[-4:]== ".pts":
-        #print "loading " + ofile
-        dict = load_mesh.loadASCII.load_points_file(ofile)
-        #print "dict",dict
+    if (ofile[-4:]== ".xya" or ofile[-4:]== ".pts"):
+        dict = load_mesh.loadASCII.import_points_file(ofile)
         dict['points'] = dict['pointlist']
         dict['outline_segments'] = []
@@ -2765 +2745 @@
         dict['holes'] = [] 
         newmesh= Mesh(geo_reference = dict['geo_reference'])
-        newmesh.IOOutline2Mesh(dict) #FIXME use IOOutline2Mesh
+        newmesh.IOOutline2Mesh(dict)
         counter = newmesh.removeDuplicatedUserVertices()
         if (counter >0):
-            print "%i duplicate vertices removed from dataset" % (counter)  
+            print "%i duplicate vertices removed from dataset" % (counter)
     elif (ofile[-4:]== ".tsh" or ofile[-4:]== ".msh"):
-        dict = load_mesh.loadASCII.import_triangulation(ofile)
+        dict = load_mesh.loadASCII.import_mesh_file(ofile)
         #print "********"
         #print "zq mesh.dict",dict