Skip to content
Snippets Groups Projects

Smooth global plots

Merged Stefan Mauerberger requested to merge central_longitue into master
1 file
+ 5
0
Compare changes
  • Side-by-side
  • Inline
+ 5
0
@@ -1047,6 +1047,11 @@ class Inversion:
Hannes' library.
"""
theta, phi = zip(*pyfield.equi_sph(n))
# Phi's range is [0,2*pi]. To achieve smooth global plots with cartopy
# the interval depends on the central longitude (cl) chosen. AFAIK it
# shall be [-180+ct,180+ct] to have points that hit the projection's
# edges. If we agree on cl=0, subtracting pi shall do the job.
phi = np.subtract(phi, np.pi)
colat, lon = np.rad2deg((theta, phi))
n_desi = len(colat) # Actual number of design points
# Allocate array and store locations as expected by Hannes libraries
Loading