Hi,
I'm working with data on the cubed sphere grid
(http://thatsmaths.files.wordpress.com/2012/10/cubed-sphere-grid.jpg).
Even though the data is interpolated on a regular lat-lon grid
for plotting, I would like to plot the edges of the 6 cube faces
for visual reference.
Moreover I would like this to be independent from the map projection type
and map center lon/lat, as you can see in this image:
https://www.dropbox.com/s/79rzhv1buhuk2b8/Desired_cubed_sphere_edges.png
But all I can get using Basemap is this:
https://www.dropbox.com/s/09j62d270znyego/Basemap_cubed_sphere_edges.png
which is definitely not what I would like to do.
This image is produced by the attached script cubed_sphere_edges.py .
Does anyone have any advice on how to solve this issue?
Any help would be really appreciated.
cubed_sphere_edges.py (2.37 KB)
ยทยทยท
--
Pier Giuseppe Fogli
CMCC - Centro Euro-Mediterraneo sui Cambiamenti Climatici (www.cmcc.it)
Viale Aldo Moro, 44
40127 Bologna
ITALY
Phone: +39 051 3782606
FAX: +39 051 3782655
e-mail: piergiuseppe DOT fogli AT cmcc DOT it
skype: beppecmcc