Documentation of pos_vs_neg_LPCT

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

Help text

  The following files are in the directory below:

  filin = 'temp_A_L1_1000_years.nc'    % SST
  filin = 'temp_A_L5_1000_years.nc'    % Layer 5 Temp.
  filin = 'temp_A_L1-5_1000_years.nc'  % Averaged Pot. Temp., Top 160m

Cross-Reference Information

Listing of script pos_vs_neg_LPCT

clear
  filin = 'temp_A_L1-10_1000_years.nc' % Averaged Pot. Temp., Top 620m

cd /home/disk/hayes2/dvimont/csiro/data

nc = netcdf(filin, 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);

%  ctlim = [130 285 -50 65];
%  ctlim = [180 270 -6 6];
%  ctlim = [110 275 -75 65];
  ctlim = [-.1 360 -90 90];
  [xk, yk] = keep_var(ctlim, lon, lat);

%  temp = nc{'temp'}(:,1,yk,xk);
  temp = nc{'temp'}(:,yk,xk);
  mv = nc{'temp'}.missing_value(:);
nc = close(nc);
temp = squeeze(temp);
[ntim, nlat, nlon] = size(temp);

temp = reshape(temp, ntim, nlat*nlon);
for i = 1:ntim
  temp(i,find(temp(i,:) == mv)) = NaN * ones(size(find(temp(i,:) == mv)));
end;
clim = mean(temp);
kp = find(~isnan(clim));
temp = temp(:,kp);
temp = detrend(temp);

lat = lat(yk); lon = lon(xk);
get_global
default_global; FRAME = [0 360 -90 90];

cd /home/disk/tao/dvimont/matlab/CSIRO/Data
load butter_4.5_ctstar.mat
ctlow = detrend(ctlow);
ctlow = (ctlow - mean(ctlow)) / std(ctlow);

pos = find(ctlow > 0);
neg = find(ctlow < 0);

pospat = NaN * ones(1, nlat*nlon);
negpat = NaN * ones(1, nlat*nlon);

pospat(kp) = ctlow(pos)' * temp(pos, :) ./ sum(ctlow(pos));
negpat(kp) = ctlow(neg)' * temp(neg, :) ./ sum(ctlow(neg));

pospat = reshape(pospat, nlat, nlon);
negpat = reshape(negpat, nlat, nlon);
pospat(1,1) = -.4; pospat(1,2) = .4;
negpat(1,1) = -.4; negpat(1,2) = .4;

figure(1); figure_orient;
sp(1)
     gshade(pospat, [-.4:.05:.4]);
     dc;
     colorbar
sp(2)
     gshade(negpat, [-.4:.05:.4]);
     dc;
     colorbar