pyacs.lib.icosahedron.mesh_regional

Build regional icosahedron mesh within bounds.

pyacs.lib.icosahedron.mesh_regional.mesh_regional(num_subdivisions=6, bounds=None)[source]

Build an equilateral triangle mesh on the unit sphere within bounds.

Uses successive subdivisions of an icosahedron; only faces intersecting the given bounds are kept.

Parameters:

  • num_subdivisions (int, optional) – Number of subdivision iterations. Default is 6.

  • bounds (str, optional) – Region as string, e.g. ‘/lon_min/lon_max/lat_min/lat_max’ (degrees).

Returns:

  • verts (list) – Vertices as [x, y, z] in geocentric cartesian on unit sphere.

  • faces (list) – Faces as lists of vertex indices; faces[i][j] is vertex j of face i.