Global Index (short | long) | Local contents | Local Index (short | long)
xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
This script calls | |
---|---|
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