Documentation of plot_regressions_HP_LP

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

Help text

load RAW_detrend_L1-7_EOF_yr101-1000.mat; rpc = -1*pcs;
load LP9_detrend_L1-7_EOF_yr101-1000.mat; lppc = -1*pcs;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; hppc = 1*pcs;

Cross-Reference Information

This script calls
axis_limits color_shade dc2 default_global detrend drawmapm figure_orient fl fo frame gcont getflx getheat getll getnc globll gquiv keep_var mcont3 mean2 mshade regress_eof sph_curl1 subplot

Listing of script plot_regressions_HP_LP


clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data

%load PDO_detrend_L1_EOF_yr101-1000.mat; rpc = -1*pcs;
%load PDO_LP9_detrend_L1_EOF_yr101-1000.mat; lppc = -1*pcs;
%load PDO_HP10_detrend_L1_EOF_yr101-1000.mat; hppc = 1*pcs;

load GR9_detrend_L1_EOF_yr101-1000.mat; lppc = -1*pcs;
load GR10_detrend_L1_EOF_yr101-1000.mat; hppc = -1*pcs;

lims = [-0.1 360 -90 90];
tim = 101:1000;
lev = 1;
[lat1, lon1, depth1, lm] = getll('temp', lims);
colr = 0.55;

%varn = 'temp'; tit = 'SST'; units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25;
%varn = 'psl'; tit = 'SLP'; units = 'hPa'; cint = .1; pn = -1; cint2 = 0.25;
%varn = 'rnd'; tit = 'Precipitation'; units = 'mm day^-^1'; cint = .075; pn = -1;
   cint2 = 0.1;
%varn = 'z250'; tit = '250mb Height'; units = 'm'; cint = 2; pn = -1; cint2 = 0.25;
%varn = 'ico'; tit = 'Ice Concentration'; units = '%'; 
%       cint = 1; pn = 1; cint2 = 0.25;
%varn = 'taux'; tit = 'Wind Stress'; units = 'dynes cm^-^2'; pn = 1; cint2 = 0.1;
%varn2 = 'tauy';

varn = 'u'; tit = '80-270m Averaged Ocean Current'; 
       units = 'cm s^-^1'; pn = 1; cint2 = 0.1;
varn2 = 'v';

varn = 'temp'; tit = '0-270m HC'; units = 'x 10^6 J m^-^2'; cint = 50;
       pn = 1; cint2 = 0.25;

var = getheat(lims, 1:7, tim, varn);
var = 1e-6*detrend(var);
var2 = getheat(lims, 4:7, tim, varn2);
var2 = 1e0*detrend(var2);

var = getflx(varn, lims, tim);
var = detrend(var);

var = getnc(varn, lims, 1, tim);
var = detrend(var);

var2 = getnc(varn2, lims, 1, tim);
var3 = sqrt(var2.^2 + var.^2);
var = detrend(var);
var2 = detrend(var2);
var3 = detrend(var3); 

[hreg, hcoef] = regress_eof(var, hppc(:,2), 0);
[lreg, lcoef] = regress_eof(var, lppc, -1);
[lat, lon] = getll(varn, lims);

%  Plot the data

figure(1); figure_orient(1);
default_global;  FRAME = [0 360 -90 90];

subplot(2,1,1);
  gcont(hreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(hcoef.^2), cint2, 0.87*[1 1 1]);
  title(['HP10:  ' tit ' Regressed on HP10 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1;  ' ...
          'Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);

subplot(2,1,2);
  gcont(lreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(lcoef.^2), cint2, 0.87*[1 1 1]);
  title(['LP9:  ' tit ' Regressed on LP9 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1;  ' ...
          'Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);

cd ~/Thesis/Chap2
%  print -dps2 HC270_HP_LP_global.ps




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Polar Plots
%  Plot the data

figure(2); figure_orient(1);
default_global;  

subplot(2,2,1);  FRAME = [0 360 40 90];
  mcont3(100*hreg, cint, 'stereo', [90 270]); gridm on;
  drawmapm; 
%  [c, h] = mshade3(squeeze(hcoef.^2), cint2, 0.87);
%  hold on; 
%    mcont3(100*hreg, cint, 'stereo', [90 270]); gridm on;
%  hold off;
  axis_limits;
  title(['NH:  ' tit ' Regressed on HP10 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);
subplot(2,2,2);  FRAME = [0 360 -90 -40];
  mcont3(100*hreg, 1*cint, 'stereo', [-90 270]); gridm on;
  drawmapm; 
%  [c, h] = mshade(squeeze(hcoef.^2), cint2, 0.87);
%  hold on; 
%    mcont3(1*hreg, cint, 'stereo', [-90 270]); gridm on;
%  hold off;
  axis_limits;
  title(['SH:  ' tit ' Regressed on HP10 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);

subplot(2,2,3);  FRAME = [0 360 40 90];
  mcont3(100*lreg, 1*cint, 'stereo', [90 270]); gridm on;
  drawmapm; 
%  [c, h] = mshade(squeeze(lcoef.^2), cint2, 0.87);
%  hold on; 
%    mcont3(1*lreg, cint, 'stereo', [90 270]); gridm on;
%  hold off;
  axis_limits;
  title(['NH:  ' tit ' Regressed on LP9 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);
subplot(2,2,4);  FRAME = [0 360 -90 -40];
  mcont3(100*lreg, 1*cint, 'stereo', [-90 270]); gridm on;
  drawmapm; 
  [c, h] = mshade(squeeze(lcoef.^2), cint2, 0.87);
  hold on; 
    mcont3(1*lreg, cint, 'stereo', [-90 270]); gridm on;
  hold off;
  axis_limits;
  title(['SH:  ' tit ' Regressed on LP9 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);


cd ~/Thesis/Chap2
%  print -dps2 SICE_HP_LP_global.ps


%  Plot wind stress 333333333333333333333333333333333

[hreg2, hcoef2] = regress_eof(var2, hppc, 0);
[lreg2, lcoef2] = regress_eof(var2, lppc, 0);
[hreg3, hcoef3] = regress_eof(var3, hppc, 0);
[lreg3, lcoef3] = regress_eof(var3, lppc, 0);

figure(1); figure_orient(1);
default_global;  FRAME = [0 360 -90 90];

subplot(2,1,1);
  gquiv(hreg, hreg2, -.3, 2, units);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(hcoef3.^2), cint2, 0.87*[1 1 1]);
  title(['HP10:  ' tit ' Regressed on HP10 PC1']);
  xl = get(gca, 'xlabel'); xl = get(xl, 'string');
  xlabel([xl ' std^-^1;  Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);
subplot(2,1,2);
  gquiv(lreg, lreg2, -.3, 2, units);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(lcoef3.^2), cint2, 0.87*[1 1 1]);
  title(['LP9:  ' tit ' Regressed on LP9 PC1']);
  xl = get(gca, 'xlabel'); xl = get(xl, 'string');
  xlabel([xl ' std^-^1;  Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);

cd ~/Thesis/Chap2
%  print -dps2 UV_80-270m_HP_LP_global.ps


ct = getnc('temp', [180 270 -6 6], 1, 101:1000);
ct = squeeze(mean(mean(ct, 2), 3));
ct = detrend(ct) ./ std(detrend(ct));






%%%%%%%  Get seasonal data

cd /home/disk/hayes2/dvimont/csiro/annual_data
varn = 'temp';

filin = 'djf_L1_temp.nc'; tit = 'DJF SST'; tit3 = 'CT\_DJF';
     units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25;

filin = 'jja_L1_temp.nc';tit = 'JJA SST'; tit3 = 'CT\_JJA';
     units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25;


filin = 'ndjfma_L1_temp.nc';tit = 'NDJFMA SST'; tit3 = 'CT\_NDJFMA';
     units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25;


filin = 'mjjaso_L1_temp.nc';tit = 'MJJASO SST'; tit3 = 'CT\_MJJASO';
     units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25;


nc = netcdf(filin, 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  [xk, yk] = keep_var(lims, lon, lat);

  var = nc{'temp'}(:, 1, yk, xk);
  mv = nc{'temp'}.missing_value(:);
nc = close(nc);
var(var == mv) = NaN;
var = detrend(var);
var = squeeze(var);
lat = lat(yk); lon = lon(xk);

[ntim, nlat, nlon] = size(var);
var = [var(2:ntim,:,:); var(1,:,:)];

[xk, yk] = keep_var([180 270 -6 6], lon, lat);
ct = squeeze(mean2(mean2(shiftdim(var(:,yk,xk), 1))));

[b1, a1] = butter(9, 2/10);
[b2, a2] = butter(9, 2/9);
hppc = ct - filtfilt(b1, a1, ct); stdh = std(hppc);
lppc = filtfilt(b2, a2, ct); stdl = std(lppc);
hppc = detrend(hppc) ./ (std(detrend(hppc)));
lppc = detrend(lppc) ./ (std(detrend(lppc)));


[hreg, hcoef] = regress_eof(var, hppc, 0);
[lreg, lcoef] = regress_eof(var, lppc, 0);
[lat, lon] = getll(varn, lims);

%  Plot the data

figure(1); figure_orient(1);
default_global;  FRAME = [0 360 -90 90];

subplot(2,1,1);
  gcont(hreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(hcoef.^2), cint2, 0.87*[1 1 1]);
  title(['HP10:  ' tit ' Regressed on HP10 ' tit3]);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1;  ' ...
          'Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);
  ylabel(['STD HPPC = ' num2str(round(100*stdh)/100) ' K std^-^1']);

subplot(2,1,2);
  gcont(lreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
  color_shade(squeeze(lcoef.^2), cint2, 0.87*[1 1 1]);
  title(['LP9:  ' tit ' Regressed on LP9 ' tit3]);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1;  ' ...
          'Light Shading indicates ' num2str(100*cint2) ...
          '% Unfiltered Variance Explained']);
  ylabel(['STD LPPC = ' num2str(round(100*stdl)/100) ' K std^-^1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Seasonal
%  print -dps2 SST_MJJASO_HP_LP_global.ps











%%%%%%%%%%%%%%%  Test to make sure everything's OK

nc = netcdf(filin, 'nowrite');
  tim = nc{'month'}(:);
  lreg = nc{'temp'}(1,1,:,:);
  mv = nc{varn}.missing_value(:);
nc = close(nc);
lreg(lreg == mv) = NaN;

nc = netcdf('temp_M_L1_1000_years_new.nc', 'nowrite');
  tim = nc{'month'}(:);
  tind = find(tim == 11988)
  hreg = nc{'temp'}(tind+[0:2], 1, :, :);
nc = close(nc);
hreg(hreg == mv) = NaN;

figure(1); fo
subplot(2,1,1)
     gcont(mean2(hreg), [-2:2:40]);
     dc2(mean2(hreg));
subplot(2,1,2)
     gcont(lreg, [-2:2:40]);
     dc2(lreg);

figure(2); fl(1)
     gcont(squeeze(mean(hreg))-squeeze(lreg), [-5:.2:5]);
     dc2(lreg)


%  Plot wind stress curl

lreg1 = lreg;
hreg1 = hreg;

[lreg, lat2, lon2] = sph_curl1(lreg1, lreg2, lat, lon, 1);
[hreg, lat2, lon2] = sph_curl1(hreg1, hreg2, lat, lon, 1);

lat = lat2; lon = lon2;
lreg = 1e8 * lreg; hreg = 1e8*hreg;
units = 'x 10^-^8 dynes cm^-^2 m^-^1'; tit = 'Wind Stress Curl'; cint = 1.5; pn = 1;

%  Plot the data

figure(1); figure_orient(1);
default_global;  FRAME = [0 360 -90 90];

subplot(2,1,1);
  gcont(hreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
%  color_shade(squeeze(hcoef.^2), cint2, 0.87*[1 1 1]);
  title(['HP10:  ' tit ' Regressed on HP10 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);%  ;  ' ...
%          'Light Shading indicates ' num2str(100*cint2) ...
%          '% Unfiltered Variance Explained']);

subplot(2,1,2);
  gcont(lreg, cint);
  XAX = lon1; YAX = lat1;
  dc2(lm, colr, pn); globll;
%  color_shade(squeeze(lcoef.^2), cint2, 0.87*[1 1 1]);
  title(['LP9:  ' tit ' Regressed on LP9 PC1']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units ' std^-^1']);%  ;  ' ...
%          'Light Shading indicates ' num2str(100*cint2) ...
%          '% Unfiltered Variance Explained']);

cd ~/Thesis/Chap2
%  print -dps2 HC270_HP_LP_global.ps