Documentation of regress_velocity

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

Cross-Reference Information

Listing of script regress_velocity

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 4:7;
u = getheat(lims, lev, tim, 'u');
v = getheat(lims, lev, tim, 'v');
[lat, lon, depth] = getll('u', lims);
[tem, tem1, depthw] = getll('wl', lims);
depthw
ureg = regress_eof(u, pcs, 0);
vreg = regress_eof(v, pcs, 0);
default_global;  FRAME = [105 290 -60 60];
figure(1); figure_landscape;
subplot(2,2,1);
  gquiv(ureg, vreg, -0.25, 1, 'cm sec^-^1');
  dc2(ureg);
  title('LP10:  80-270m OCEAN VELOCITY');
subplot(2,2,3);
  gquiv(mean(u), mean(v), -7.5, 1, 'cm sec^-^1');
  dc2(ureg);
  title('CLIM:  80-270m OCEAN VELOCITY');
lev = 8:10;
u = getheat(lims, lev, tim, 'u');
v = getheat(lims, lev, tim, 'v');
[lat, lon, depth] = getll('u', lims);
ureg = regress_eof(u, pcs, 0);
vreg = regress_eof(v, pcs, 0);
subplot(2,2,2);
  gquiv(ureg, vreg, -0.1, 1, 'cm sec^-^1');
  dc2(ureg);
  title('LP10:  270-470m OCEAN VELOCITY');
subplot(2,2,4);
  gquiv(mean(u), mean(v), -3, 1, 'cm sec^-^1');
  dc2(ureg);
  title('CLIM:  270-470m OCEAN VELOCITY');
subplot(2,2,1);
  ylabel('Lag = 0');
subplot(2,2,2);
  ylabel('Lag = 0');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 0 65];
tim = 101:1000;
lev = 1:3;
v = getheat(lims, lev, tim, 'v');
[lat, lon] = getll('v', lims);
[strf, lat_out, lon_out] = vel_to_strfnc(squeeze(mean(v)), lat, lon);
default_global; FRAME = [105 290 2.5 60];
XAX = lon_out; YAX = lat_out;
figure(1); figure_orient;
subplot(2,1,1);
  gcont(1e-6*strf, [-50:5:50]);
  dc2(strf);
  title('0-80m Climatological Stream Function');
lev = 1:7;
v = getheat(lims, lev, tim, 'v');
[lat, lon] = getll('v', lims);
[strf, lat_out, lon_out] = vel_to_strfnc(squeeze(mean(v)), lat, lon);
subplot(2,1,2);
  gcont(1e-6*strf, [-50:1:50]);
  dc2(strf);
  title('0-270m Climatological Stream Function');