Index: /issm/trunk-jpl/src/m/mesh/roundmesh.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/roundmesh.py	(revision 22319)
+++ /issm/trunk-jpl/src/m/mesh/roundmesh.py	(revision 22320)
@@ -20,12 +20,10 @@
 
 	#Get number of points on the circle
-	pointsonedge=np.floor((2.*np.pi*radius) / resolution)
+	pointsonedge=np.floor((2.*np.pi*radius) / resolution)+1 #+1 to close the outline
 
 	#Calculate the cartesians coordinates of the points
-	x_list=np.ones(pointsonedge)
-	y_list=np.ones(pointsonedge)
-	theta=np.linspace(0.,2.*np.pi,num=pointsonedge,endpoint=False)
-	x_list=roundsigfig(radius*x_list*np.cos(theta),12)
-	y_list=roundsigfig(radius*y_list*np.sin(theta),12)
+	theta=np.linspace(0.,2.*np.pi,pointsonedge)
+	x_list=roundsigfig(radius*np.cos(theta),12)
+	y_list=roundsigfig(radius*np.sin(theta),12)
 	A=OrderedDict()
 	A['x']=[x_list]