Documentation of ocean_velocity_vert_ave

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])

Cross-Reference Information

Listing of script ocean_velocity_vert_ave

clear
lev = 1:3;
tim = 101:550;
[u, v] = getnc('u', 'v', [-0.1 360 -90 90], lev, tim);
[lat, lon, depth] = getll('u', [-0.1 360 -90 90]);
u = squeeze(mean(u));
v = squeeze(mean(v));
[nlev, nlat, nlon] = size(u);
d1 = depth(1);
depthw = [];
for i = 1:nlev;
  depthw = [depthw; depth(i)+d1];
  d1 = depth(i+1) - depth(i) - d1;
end
dz = diff([0; depthw]);
u = reshape(u, nlev, nlat*nlon);
v = reshape(v, nlev, nlat*nlon);
u = dz'*u ./ sum(dz);
v = dz'*v ./ sum(dz);
u = reshape(u, nlat, nlon);
v = reshape(v, nlat, nlon);
u = u - u2;
v = v - v2;
default_global;
figure(1);
sp(1)
     [h, hh] = gquiv(u, v, -20, 1, 'cm s^-^1');
     dc2(u);
     title('0:80m Vertically Averaged Ocean Current, years 101-550');

cd /home/disk/tao/dvimont/matlab/CSIRO/Plots_Clim
print -dps2 oce_vel_L1-3_L4-7_yr101-550.ps