Index: /issm/trunk-jpl/src/m/mesh/planet/gmsh/gmshplanet.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/planet/gmsh/gmshplanet.py	(revision 21047)
+++ /issm/trunk-jpl/src/m/mesh/planet/gmsh/gmshplanet.py	(revision 21048)
@@ -155,4 +155,5 @@
 		A=fid.readline().split()
 		mesh.elements[i]=[int(A[5]),int(A[6]),int(A[7])]
+	mesh.elements=mesh.elements.astype(int)
 	A=fid.readline().strip()
 	if not strcmp(A,'$EndElements'):
Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 21047)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 21048)
@@ -23,6 +23,4 @@
 md.slr.sealevel=zeros((md.mesh.numberofvertices,1))
 #antarctica
-print type(md.mesh.elements)
-print md.mesh.x[md.mesh.elements-1]
 late=numpy.sum(md.mesh.lat[md.mesh.elements-1],axis=1)/3
 longe=numpy.sum(md.mesh.long[md.mesh.elements-1],axis=1)/3