Changeset 870

Timestamp:

05/17/19 13:41:18 (5 years ago)

Author:

dubos

Message:

devel/Python : bugfix cell bounds, last vertex of pentagons must be repeated

File:

• codes/icosagcm/devel/Python/dynamico/meshes.py (modified) (1 diff)

codes/icosagcm/devel/Python/dynamico/meshes.py

@@ -411 +411 @@
             vertex_lon, vertex_lat = np.zeros((num, maxdeg)), np.zeros((num, maxdeg))
             for cell in range(num):
-                for ivertex in range(deg[cell]):
-                    vtx = vertex[cell,ivertex] - 1
-                    vertex_lon[cell, ivertex] = lon[vtx]
-                    vertex_lat[cell, ivertex] = lat[vtx]
+                cdeg=deg[cell]
+                for ivertex in range(maxdeg):
+                    if ivertex<cdeg:
+                        vtx = vertex[cell,ivertex] - 1
+                        vlon, vlat = lon[vtx], lat[vtx]
+                    vertex_lon[cell, ivertex] = vlon
+                    vertex_lat[cell, ivertex] = vlat
             return vertex_lon, vertex_lat
         mesh.primal_bounds_lon, mesh.primal_bounds_lat = bounds(mesh.primal_deg, mesh.primal_vertex, mesh.lon_v, mesh.lat_v)