| 1 |
cnh |
1.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]']) |
Parent Directory
| 
Revision Log
| 
Revision Graph