%
% function mapxy_field(X,Y,FIELD,lax,palette,titre,NFIG)
%

     function mapxy_field(X,Y,FIELD,lax,palette,titre,NFIG)

path('/u/u0/czaja/MITGCM/MATLAB',path);
path('/u/u0/czaja/MATLAB',path);

  figure(NFIG), clf, 

% Plot 2D field
  pcolor(X,Y,FIELD'); %shading flat
  set(gca,'Fontsize',14);
  caxis(lax);
  colormap(palette);
  colorbar;
  title(titre,'Fontsize',20);
  xlabel('Longitude','Fontsize',15);
  ylabel('Latitude','Fontsize',15);
  hold on

%Plot continents
  x1 = 105+2.5;
  x2 = 150+2.5;
  x3 = 210+2.5;
  x4 = 255+2.5;
  y1 = 64;
  y2 = -30+1;
  [h]=line([x1 x1],[y1 y2]); set(h,'Linewidth',2,'color','k');
  [h]=line([x1 x2],[y2 y2]); set(h,'Linewidth',2,'color','k');
  [h]=line([x2 x2],[y2 y1]); set(h,'Linewidth',2,'color','k');
  [h]=line([x3 x3],[y1 y2]); set(h,'Linewidth',2,'color','k');
  [h]=line([x3 x4],[y2 y2]); set(h,'Linewidth',2,'color','k');
  [h]=line([x4 x4],[y2 y1]); set(h,'Linewidth',2,'color','k');

  axis([0 360 -64 64]);