Hi,
Having been a long-term m_map user, I've recently switched to using Matlab's mapping toolbox (v2013a at the moment). On the whole, it's pretty great, and it more-or-less replaces m_map's functionality.
In particular, I've been using it to generate quick, and rather nice, plots of geographical output from a model that I use, NEMO. As aficionados will know, this ocean model has a strange tripolar grid that's regular in the southern hemisphere, but distorts to two poles in the north to avoid the dreaded north pole singularity. When plotting, I use pcolorm as follows:
Where yy is a 2D matrix of latitude locations, xx is a 2D matrix of longitude locations and fld is whatever 2D matrix of model properties I'm interested in at the time. The matrices xx and yy have to be 2D because of this distorted grid, and they are usually the position (in latitude and longitude degrees) of the centre of a given grid cell.
I can't find anything on the Matlab webpages to help me, but I suspect I shouldn't be using the position of grid cell centres for this sort of plot, and should instead be doing something like using the southwest corner of a grid cell. Or, in fact, doing something completely different in which I formally specify the grid cell edges.
Does anyone know of any resources that would explain what the best approach is? The Matlab page is a bit too thin in detail to my mind - especially as a geographical example is used to illustrate it. And though the impact may be less significant, the question of grid cell centres vs. grid cell edges occurs even on a regular grid.
Cheers,
Andrew.
