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