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	edhill	1.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