mdpa working with exponential notation

Author: e-dub

pygem/mdpahandler.py

@@ -44,6 +44,7 @@ def parse(self, filename):
                     if line.startswith('End Nodes'):
                         break
                     else:
+                        line = line.replace('D', 'E')
                         li = []
                         for t in line.split()[1:]:
                             try:
@@ -79,11 +80,18 @@ def write(self, mesh_points, filename):
                     if num == index + 1:
                         if line.startswith('End Nodes'):
                             index = -9
+                            output_file.write(line)
                         else:
+                            output_file.write(3 * ' ' + '{:5d}'.format(i+1))
                             for j in range(0, 3):
-                                line = (line[:5 + 15 * (j)] +
-                                        '{:15.10f}'.format(mesh_points[i][j]) +
-                                        line[5 + 15 * (j + 1):])
+                                output_file.write(3 * ' ' + '{:.16E}'.format(
+                                    mesh_points[i][j]))
+                            output_file.write('\n')
+#
+#                                line = (line[:5 + 15 * (j)] +
+#                                        '{:15.10f}'.format(mesh_points[i][j]) +
+#                                        line[5 + 15 * (j + 1):])
                             i += 1
                             index = num
-                    output_file.write(line)
+                    else:
+                        output_file.write(line)