Documentation of cartoon2

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

Help text

tx1 = regress_eof(txann, lpcs, 0);
ty1 = regress_eof(tyann, lpcs, 0);
tau = sqrt(tx1.^2 + ty1.^2);

Cross-Reference Information

Listing of script cartoon2

clear
csirod
load ML_ANN_slp.mat;
load ML_ANN_sst.mat;
load taux_ML_annave.mat;
load tauy_ML_annave.mat;
back
data
load ML_SLP_eof_nhem.mat
back
slp1 = regress_eof(slp, lpcs, 0);
sst1 = regress_eof(sst2, lpcs, 0);

lims = [-0.1 360 -90 90];
[lat1, lon1, depth1, lm] = getll('temp', lims);
lm2 = ~isnan(lm);

figure(1); fo; clf;

default_global; FRAME = [100 300 -45 75];
sptalk(3,1,1); cla;
  [h, c] = mcont2(-slp1, .16, 'mollweid');
  set(c, 'linewidth', 2);
  gridm on;
  drawmapm;
  axis_limits;
  hold on;
    mshade(-1*squeeze(sst1).*lm2, .08, .39*[1 1 1]);%[1 .4 .4]); %
    mshade(1*squeeze(sst1).*lm2, .08, .8*[1 1 1]);%[.4 .4 1]); %
%    mshade(-1*squeeze(sst1).*lm2, .2, [1 .2 .2]);
%    mshade(squeeze(sst1).*lm2, .2, [.2 .2 1]);
 hold off;

%cd ~/Thesis/Talk
%print -dpsc cartoon1.ps


loadpcs;

lims = [95 305 -90 90];
tim = 101:1000;

sst = getnc('temp', lims, 1, tim);
slp = getflx('psl', lims, tim);

sst2 = regress_eof(sst, lpcs, 0);
slp2 = regress_eof(slp, lpcs, 0);

[lat, lon, depth1, lm] = getll('temp', lims);
lm2 = ~isnan(lm);
default_global; FRAME = [100 300 -45 75];

%figure(1); fo; clf;

sptalk(3,1,2); cla
  [h, c] = mcont2(slp2, .16, 'mollweid');
  set(c, 'linewidth', 2);
  gridm on;
  drawmapm;
  axis_limits;
  hold on;
    mshade(1*squeeze(sst2).*lm2, .08, .39*[1 1 1]);%[1 .4 .4]);%
    mshade(-1*squeeze(sst2).*lm2, .08, .8*[1 1 1]);%[.4 .4 1]);%
%    mshade(1*squeeze(sst2).*lm2, .175, [1 .2 .2]);
%    mshade(-1*squeeze(sst2).*lm2, .175, [.2 .2 1]);
  hold off;



cd ~/Thesis/Chap5
print -dps2 cartoon2bw.ps