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 -30 75];
sptalk(3,1,1);
  [h, c] = mcont2(-slp1, [-1:.2:-.2 .2:.2:1], 'mollweid');
  set(c, 'linewidth', 2);
  gridm on;
  drawmapm;
  axis_limits;
  hold on;
    mshade(-1*squeeze(sst1).*lm2, .05, [1 .6 .6]);
    mshade(squeeze(sst1).*lm2, .05, [.6 .6 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 -45 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 -30 75];

figure(1); fo; clf;

sptalk(3,1,1);
  [h, c] = mcont2(slp2, [-1.5:.15:-.15 .15:.15:1.5], 'mollweid');
  set(c, 'linewidth', 2);
  gridm on;
  drawmapm;
  axis_limits;
  hold on;
    mshade(1*squeeze(sst2).*lm2, .05, [1 .6 .6]);
    mshade(-1*squeeze(sst2).*lm2, .05, [.6 .6 1]);
    mshade(1*squeeze(sst2).*lm2, .175, [1 .2 .2]);
    mshade(-1*squeeze(sst2).*lm2, .175, [.2 .2 1]);
  hold off;



cd ~/Thesis/Talk
print -dpsc cartoon2.ps