How can I superimpose data on a map using Mapping Tool box?

When I superimpose a plot over the world map, the map gets distorted. By distorted, I mean that parts of the boarder disappear. This is the code I use:

coast = load('coast');
axesm('pcarree', 'Origin', 180,'MLabelParallel','south'); 
framem on; gridm on; mlabel on; plabel on; showaxes;
h = geoshow(coast.lat, coast.long,'Color','black');
hold on;

yy = linspace(-pi,pi,NumLong);
xx = linspace(-pi/2, pi/2, NumLat);

[XX, YY] = meshgrid(xx,yy);
surf(YY, XX, data);

[Nancy]

The map is still there and is the right size but when I plot data parts of the boundaries disappear.

Answer 1

[Walter Roberson]

Before you do the surf(), set the axis XLimMode and YLimMode to 'manual'. Otherwise if the surf() wants to expand the xlimit or ylimit, it is not going to take into account the map's existence in figuring out where the bounds should be.

[Nancy]

Thank-you so much for your suggestion. I learned a lot of useful information by looking into this. However, this did not resolve my problem.