Documentation of regress_thermocline

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script regress_thermocline

clear

%for biff2 = 1:1;
biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
  load LP10_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
elseif biff2 == 2;
  load HP8_L1-7_CEOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
elseif biff2 == 3;
  load RAW_L1-7_CEOF.mat; tit = 'Unfiltered Data';
end

nfrm = 6; 
tim = 101:550;
lev = 1:10;
%lims = [105 310 -5 5];
%lims = [105 310 10 20];
%lims = [150 180 -60 60];
%lims = [120 150 -60 60];
lims = [180 210 -60 60];
%lims = [210 240 -60 60];

temp = getnc('temp', lims, lev, tim);
[lat, lon, depth] = getll('temp', lims);
[ntim, nlev, nlat, nlon] = size(temp);

[temp, clim] = remove_mean(temp);
temp = shiftdim(temp, 3);
temp2 = squeeze(mean2(temp));
%temp2 = shiftdim(temp2, 1);
[treg, temtim] = regress_ceof(temp2, pcs, nfrm);

%  Plot the data

figure(biff2); figure_orient(1);
cint = 0.025; clev = [-5:cint:-cint cint:cint:5];
lag = 0;
lon = lat;

for i = 1:nfrm;
  tem = squeeze(treg(i,:,:));
  subplot(6,1,i);
    pncont(lon, -1*depth, tem, clev, 0, 'k');
    dcv(lon, -1*depth, tem);
    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' cm s^-^1 std^-^1']);
  end
%  axis([105 299 -500 0]);
%  set(gca, 'XTick', 105:15:300, ...
%           'YTick', -500:125:0, 'YTickLabel', ['500'; '   '; '250'; '   '; '  0']);
   axis([-60 60 -500 0]);
   set(gca, 'XTick', -60:15:60, ...
           'YTick', -500:125:0, 'YTickLabel', ['500'; '   '; '250'; '   '; '  0']);
  box on
end
subplot(6,1,1)
  title([tit ':  Meridional Velocity Anomalies (' num2str(lims(1)) ...
         'E to ' num2str(lims(2)) 'E)']);

if 1;
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Plot_Ceof
if biff2 == 1;
  eval(['print -dps2 LP10_w_zonal_' num2str(lims(1)) 'to' num2str(lims(2)) '.ps']);
elseif biff2 == 2;
  print -dps2 HP8_T_zonal_150to180.ps
elseif biff2 == 3; 
  print -dps2 RAW_T_zonal_150to180.ps
end
end

end