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 |
|
|
|