Index: /issm/trunk-jpl/src/m/classes/model/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14017)
+++ /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14018)
@@ -363,7 +363,5 @@
 		if md1.mesh.dimension==2:
 			[md2.mesh.vertexconnectivity]=NodeConnectivity(md2.mesh.elements,md2.mesh.numberofvertices)
-			md2.mesh.vertexconnectivity=md2.mesh.vertexconnectivity.astype(int)
 			[md2.mesh.elementconnectivity]=ElementConnectivity(md2.mesh.elements,md2.mesh.vertexconnectivity)
-			md2.mesh.elementconnectivity=md2.mesh.elementconnectivity.astype(int)
 			md2.mesh.segments=contourenvelope(md2)
 			md2.mesh.vertexonboundary=numpy.zeros(numberofvertices2,bool)
@@ -372,7 +370,5 @@
 			#First do the connectivity for the contourenvelope in 2d
 			[md2.mesh.vertexconnectivity]=NodeConnectivity(md2.mesh.elements2d,md2.mesh.numberofvertices2d)
-			md2.mesh.vertexconnectivity=md2.mesh.vertexconnectivity.astype(int)
 			[md2.mesh.elementconnectivity]=ElementConnectivity(md2.mesh.elements2d,md2.mesh.vertexconnectivity)
-			md2.mesh.elementconnectivity=md2.mesh.elementconnectivity.astype(int)
 			md2.mesh.segments=contourenvelope(md2)
 			md2.mesh.vertexonboundary=numpy.zeros(numberofvertices2/md2.mesh.numberoflayers,bool)
@@ -381,7 +377,5 @@
 			#Then do it for 3d as usual
 			[md2.mesh.vertexconnectivity]=NodeConnectivity(md2.mesh.elements,md2.mesh.numberofvertices)
-			md2.mesh.vertexconnectivity=md2.mesh.vertexconnectivity.astype(int)
 			[md2.mesh.elementconnectivity]=ElementConnectivity(md2.mesh.elements,md2.mesh.vertexconnectivity)
-			md2.mesh.elementconnectivity=md2.mesh.elementconnectivity.astype(int)
 
 		#Boundary conditions: Dirichlets on new boundary
Index: /issm/trunk-jpl/src/m/mesh/triangle.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/triangle.py	(revision 14017)
+++ /issm/trunk-jpl/src/m/mesh/triangle.py	(revision 14018)
@@ -60,8 +60,6 @@
 
 	#Now, build the connectivity tables for this mesh.
-	[md.mesh.vertexconnectivity]= NodeConnectivity(md.mesh.elements, md.mesh.numberofvertices)
-	md.mesh.vertexconnectivity=md.mesh.vertexconnectivity.astype(int)
+	[md.mesh.vertexconnectivity] = NodeConnectivity(md.mesh.elements, md.mesh.numberofvertices)
 	[md.mesh.elementconnectivity] = ElementConnectivity(md.mesh.elements, md.mesh.vertexconnectivity)
-	md.mesh.elementconnectivity=md.mesh.elementconnectivity.astype(int)
 
 	#type of model
Index: /issm/trunk-jpl/src/m/parameterization/contourenvelope.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/contourenvelope.py	(revision 14017)
+++ /issm/trunk-jpl/src/m/parameterization/contourenvelope.py	(revision 14018)
@@ -41,8 +41,6 @@
 	if numpy.size(md.mesh.vertexconnectivity,axis=0)!=md.mesh.numberofvertices and numpy.size(md.mesh.vertexconnectivity,axis=0)!=md.mesh.numberofvertices2d:
 		[md.mesh.vertexconnectivity]=NodeConnectivity(md.mesh.elements,md.mesh.numberofvertices)
-		md.mesh.vertexconnectivity=md.mesh.vertexconnectivity.astype(int)
 	if numpy.size(md.mesh.elementconnectivity,axis=0)!=md.mesh.numberofelements and numpy.size(md.mesh.elementconnectivity,axis=0)!=md.mesh.numberofelements2d:
 		[md.mesh.elementconnectivity]=ElementConnectivity(md.mesh.elements,md.mesh.vertexconnectivity)
-		md.mesh.elementconnectivity=md.mesh.elementconnectivity.astype(int)
 
 	#get nodes inside profile