diff --git a/corbass/src/inversion.py b/corbass/src/inversion.py index d47f592a0ea672277f3af35cd69b416c12b06f9e..1b9cc20a379c26f25b89ac358da53e7568309798 100644 --- a/corbass/src/inversion.py +++ b/corbass/src/inversion.py @@ -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