Documentation of ceof_regress_velocities

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

Help text

  Plot regressions on CPC1 of T20'

Cross-Reference Information

Listing of script ceof_regress_velocities

for lev = [1 5 7];
cd /home/disk/hayes2/dvimont/csiro/data
ctlim = [110 300 -65 60];
tim = 101:550;
nc = netcdf('ul_L1-10.nc', 'nowrite');
  depth = nc{'depth'}(:);
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  [xk, yk] = keep_var(ctlim, lon, lat);
  u = nc{'u'}(tim, lev, yk, xk);
  mv = nc{'u'}.missing_value(:);
nc = close(nc);
u(u == mv) = NaN;
nc = netcdf('vl_L1-10.nc', 'nowrite');
  v = nc{'v'}(tim, lev, yk, xk);
nc = close(nc);
v(v == mv) = NaN;
lat = lat(yk); lon = lon(xk); 
u = squeeze(u); v = squeeze(v);
[ntim, nlat, nlon] = size(u);
[u, uclim] = remove_mean(u);
[v, vclim] = remove_mean(v);
u = reshape(u, ntim, nlat*nlon);
v = reshape(v, ntim, nlat*nlon);
vel_spd = sqrt(u.^2 + v.^2);

for biff = [4];
cd /home/disk/hayes2/dvimont/csiro/matlab_data
if biff == 1;
  load 20deg_CEOF_HP10_yr101-550.mat
  tit = ['HP ( < 10 yr ) '];
  [b, a] = butter(6, 2/10);
  oce_vel = vel_spd - filtfilt(b, a, vel_spd);
  ptit = 'HP10';
elseif biff == 2;
  load 20deg_CEOF_LP10_yr101-550.mat;
  tit = ['LP ( > 10 yr ) '];
  [b, a] = butter(6, 2/10);
  oce_vel = filtfilt(b, a, vel_spd);
  ptit = 'LP10';
elseif biff == 3;
  load 20deg_CEOF_HP8_yr101-550.mat
  tit = ['HP ( < 8 yr ) '];
  [b, a] = butter(6, 2/8);
  oce_vel = vel_spd - filtfilt(b, a, vel_spd);
  ptit = 'HP8';
elseif biff == 4;
  load 20deg_CEOF_RAW_yr101-550.mat
  tit = ['Unfiltered Data'];
  ptit = 'RAW';
  oce_vel = vel_spd;
end
tit2 = [num2str(depth(lev)/100) 'm Ocean Velocity'];

cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('ul_L1-10.nc', 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  [xk, yk] = keep_var(ctlim, lon, lat);
nc = close(nc);
lat = lat(yk); lon = lon(xk);
default_global

timeseries = sqrt(2)*pcs(:,1)./std(pcs(:,1));

lag = 0; lg = lag*pi/180;
num = 1; nfrm = 6;

  for ind = 1:nfrm;
    wgt = conj(exp(j * ((ind-1) * pi/nfrm + lg) ));
    temtim = real(wgt .* timeseries);
    temu = temtim'*u./ntim;
    temv = temtim'*v./ntim;
    ccoef = corr_sig(temtim, oce_vel);
    temu = reshape(temu, nlat, nlon); 
    temv = reshape(temv, nlat, nlon); 
    tstat = reshape(abs(ccoef), nlat, nlon);
    subplot(3,2,ind);
      cla;
      gshade2(abs(tstat), 0.25); hold on;
      [h, a, c] = gquiv(temu, temv, 2, 2); hold off;
      dc;
      title(['Phase = ' num2str((ind-1)*180/nfrm + lag)]);
      xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
   end
   subplot(3,2,3);
     ylabel([tit ':  ' tit2 ' regressed on CPC1 of 20deg T'', years 101-550;'...
             ' (r( CPC1, |VEL| ) > 0.25 shaded)' ]);

cd /home/disk/tao/dvimont/matlab/CSIRO/Velocity/Plots1
pcomm = ['print -dps2 ' ptit '_CPC1_OCE_VEL_L' num2str(lev) '_yr101-550.ps'];
eval(pcomm)

end;  %  End biff loop
end;  %  End lev loop