1 |
% m-file: mit_plotzonalvelocity.m |
2 |
|
3 |
% select timestep |
4 |
k=kmax; |
5 |
|
6 |
up = u(:,:,:,k).*grd.pacific_hfacw; |
7 |
upzm = mit_zonalmean(up,grd.pacific_hfacw,grd.dxg); |
8 |
ua = u(:,:,:,k).*grd.atlantic_hfacw; |
9 |
uazm = mit_zonalmean(ua,grd.atlantic_hfacw,grd.dxg); |
10 |
|
11 |
caxup = [min(upzm(:)) max(upzm(:))]; |
12 |
uplev = -1:.1:1; |
13 |
if max(abs(caxup)) < .1 |
14 |
uplev = .5*uplev; |
15 |
end |
16 |
if max(abs(caxup)) < .2 |
17 |
uplev = .2*uplev; |
18 |
end |
19 |
|
20 |
caxua = [min(uazm(:)) max(uazm(:))]; |
21 |
ualev = -1:.1:1; |
22 |
if max(abs(caxua)) < .1 |
23 |
ualev = .5*ualev; |
24 |
end |
25 |
if max(abs(caxua)) < .2 |
26 |
ualev = .2*ualev; |
27 |
end |
28 |
|
29 |
ixpw = 41; |
30 |
ixpc = 53; |
31 |
ixpe = 65; |
32 |
|
33 |
ixaw = 83; |
34 |
ixac = 85; |
35 |
ixae = 2; |
36 |
|
37 |
zaxis = -grd.zc; |
38 |
zaxis = -grd.zg; |
39 |
yaxis = grd.latg; |
40 |
%zaxis = -[1:grd.nz]; |
41 |
|
42 |
figure('PaperPosition',[0.25 0.621429 8 9.75714]) |
43 |
sh = subplot(4,2,1); |
44 |
pcol(yaxis,zaxis,upzm') |
45 |
title('Pacific Ocean: zonal average') |
46 |
sh(3) = subplot(4,2,3); |
47 |
pcol(yaxis,zaxis,sq(up(ixpw,:,:))') |
48 |
title(['section at ' num2str(grd.long(ixpw)) '^{\circ}E']) |
49 |
sh(5) = subplot(4,2,5); |
50 |
pcol(yaxis,zaxis,sq(up(ixpc,:,:))') |
51 |
title(['section at ' num2str(grd.long(ixpc)) '^{\circ}E']) |
52 |
sh(7) = subplot(4,2,7); |
53 |
pcol(yaxis,zaxis,sq(up(ixpe,:,:))') |
54 |
title(['section at ' num2str(grd.long(ixpe)) '^{\circ}E']) |
55 |
sh(2) = subplot(4,2,2); |
56 |
pcol(yaxis,zaxis,uazm',ualev) |
57 |
title('Atlantic Ocean: zonal average') |
58 |
sh(4) = subplot(4,2,4); |
59 |
pcol(yaxis,zaxis,sq(ua(ixaw,:,:))') |
60 |
title(['section at ' num2str(grd.long(ixaw)) '^{\circ}E']) |
61 |
sh(6) = subplot(4,2,6); |
62 |
pcol(yaxis,zaxis,sq(ua(ixac,:,:))') |
63 |
title(['section at ' num2str(grd.long(ixac)) '^{\circ}E']) |
64 |
sh(8) = subplot(4,2,8); |
65 |
pcol(yaxis,zaxis,sq(ua(ixae,:,:))') |
66 |
title(['section at ' num2str(grd.long(ixae)) '^{\circ}E']) |
67 |
|
68 |
set(sh,'xlim',[-1 1]*30,'ylim',[-400 0]) |
69 |
set(sh,'layer','top') |
70 |
set(sh(1:2:end),'clim',[uplev(1) uplev(end)]) |
71 |
set(sh(2:2:end),'clim',[ualev(1) ualev(end)]) |
72 |
set(gcf,'currentAxes',sh(end-1));colorbar('h') |
73 |
set(gcf,'currentAxes',sh(end));colorbar('h') |
74 |
|
75 |
suptitle(['experiment ' grd.dname ... |
76 |
', timestep = ' num2str(timesteps(k)) ... |
77 |
', ' tuname ' = ' num2str(tim(k)) ', zonal velocity [m/s]']) |