MITgcm_contrib/arnaud_matlab/mapxy_field.m

%
% 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]);
  
1	%
2	% function mapxy_field(X,Y,FIELD,lax,palette,titre,NFIG)
3	%
4
5	function mapxy_field(X,Y,FIELD,lax,palette,titre,NFIG)
6
7	path('/u/u0/czaja/MITGCM/MATLAB',path);
8	path('/u/u0/czaja/MATLAB',path);
9
10	figure(NFIG), clf,
11
12	% Plot 2D field
13	pcolor(X,Y,FIELD'); %shading flat
14	set(gca,'Fontsize',14);
15	caxis(lax);
16	colormap(palette);
17	colorbar;
18	title(titre,'Fontsize',20);
19	xlabel('Longitude','Fontsize',15);
20	ylabel('Latitude','Fontsize',15);
21	hold on
22
23	%Plot continents
24	x1 = 105+2.5;
25	x2 = 150+2.5;
26	x3 = 210+2.5;
27	x4 = 255+2.5;
28	y1 = 64;
29	y2 = -30+1;
30	[h]=line([x1 x1],[y1 y2]); set(h,'Linewidth',2,'color','k');
31	[h]=line([x1 x2],[y2 y2]); set(h,'Linewidth',2,'color','k');
32	[h]=line([x2 x2],[y2 y1]); set(h,'Linewidth',2,'color','k');
33	[h]=line([x3 x3],[y1 y2]); set(h,'Linewidth',2,'color','k');
34	[h]=line([x3 x4],[y2 y2]); set(h,'Linewidth',2,'color','k');
35	[h]=line([x4 x4],[y2 y1]); set(h,'Linewidth',2,'color','k');
36
37	axis([0 360 -64 64]);
38