-fs=open(f1, 'r')
-ft=open(f3, 'w')
-ft.write('MeshVersionFormatted 1\n')
-ft.write('\nDimension\n3\n')
-tmp=fs.readline()
-nb=int(tmp)
-ft.write('\nVertices\n')
-ft.write(tmp)
-for i in xrange(0,nb):
- lig=fs.readline()
- ft.write(lig)
-fs.close()
+with codecs.open(f3, 'w') as ft:
+ with codecs.open(f1, 'r') as fs:
+ ft.write('MeshVersionFormatted 1\n')
+ ft.write('\nDimension\n3\n')
+ tmp = fs.readline()
+ nb = int(tmp)
+ ft.write('\nVertices\n')
+ ft.write(tmp)
+ for i in range(nb):
+ lig = fs.readline()
+ ft.write(lig)
+ with codecs.open(f2, 'r') as fs:
+ lig = fs.readline()
+ nb = int(lig.split()[0])
+ ft.write('\nTriangles\n' + lig.split()[0] + '\n')
+ for i in range(nb):
+ lig = fs.readline()
+ lig = lig.split()
+ ft.write(lig[1] + ' ' + lig[2] + ' ' + lig[3] + ' ' + lig[4] + '\n')
+ ft.write('\nEnd\n')