1 |
jmc |
1.1 |
% plot velocity fields |
2 |
|
|
|
3 |
jmc |
1.3 |
% $Header: $ |
4 |
jmc |
1.1 |
% $Name: $ |
5 |
|
|
|
6 |
|
|
figure |
7 |
|
|
colormap(jet0); |
8 |
|
|
for k=pkt; |
9 |
|
|
clear uc vc |
10 |
|
|
uc(1:grd.nx-1,:) = .5*(u(1:grd.nx-1,:,iz,k)+u(2:grd.nx,:,iz,k)); |
11 |
|
|
uc(grd.nx,:) = .5*(u(1,:,iz,k)+u(grd.nx,:,iz,k)); |
12 |
|
|
uc = uc.*grd.cmask(:,:,iz); |
13 |
|
|
vc(:,1:grd.ny-1) = .5*(v(:,1:grd.ny-1,iz,k)+v(:,2:grd.ny,iz,k)); |
14 |
|
|
vc(:,grd.ny) = .5*v(:,grd.ny,iz,k); |
15 |
|
|
vc = vc.*grd.cmask(:,:,iz); |
16 |
|
|
imagesc(grd.lonc,grd.latc,w(:,:,iz,k)'); caxis([-1 1]*max(abs(caxis))); |
17 |
|
|
% hold on,NaNquiver(grd.xc,grd.yc,uc,vc); hold off |
18 |
|
|
hold on, quiver(grd.xc,grd.yc,uc,vc); hold off |
19 |
|
|
axis xy; axis image |
20 |
|
|
set(gca,'clipping','off') |
21 |
|
|
set(gca,'xlim',[0 360]) |
22 |
|
|
% colorbar('h') |
23 |
|
|
title(['velocity field in 3D, depth = ' num2str(grd.zc(iz)) ]) |
24 |
|
|
suptitle(['experiment ' dname ', timestep = ' num2str(timesteps(k)) ... |
25 |
|
|
', ' tuname ' = ' num2str(tim(k))]) |
26 |
|
|
drawnow; pause(delay); |
27 |
|
|
end |
28 |
|
|
|