Changeset 22320

Timestamp:

01/03/18 11:12:58 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixed python version of roundmesh

File:

• issm/trunk-jpl/src/m/mesh/roundmesh.py (modified) (1 diff)

issm/trunk-jpl/src/m/mesh/roundmesh.py

@@ -20 +20 @@
 
         #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]