Documentation of regress_EQ_SST_indices

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

Help text

  For fun, let's see how much mass the atmosphere loses.

Cross-Reference Information

Listing of script regress_EQ_SST_indices

clean
lims = [0 360 -90 90];
tim = 1:636;
cd /home/disk/tao/data/nmc.reanalysis/monthly
[pres, lat, lon] = getnc2('pres.mon.mean.nc', 'pres', lims, 1, tim);
pres = cosweight(pres, lat);
figure_landscape(1); clf;
dmass = squeeze(mean(mean(shiftdim(pres, 1))));
plot(1948:(1/12):2000.99, dmass);


cd /home/disk/tao/data/nmc.reanalysis
lm2 = getnc2('landmask.nc', 'lsmask', lims, 1, 1);
lm2(lm2==-1) = NaN;

cd /home/disk/tao/data/nmc.reanalysis/monthly
lm1 = getnc2('land.sfc.gauss.nc', 'land', lims, 1, 1);
lm1(lm1==1) = NaN; lm1(~isnan(lm1))=1;

[sst, lat, lon] = getnc2('skt.mon.mean.nc', 'skt', lims, 1, tim);
sst = annave(sst);
[ntim, nlat, nlon] = size(sst);

[slp, lat2, lon2] = getnc2('slp.mon.mean.nc', 'slp', lims, 1, tim);
slp = annave(slp);
[ntim, nlat2, nlon2] = size(slp);


%  Now, get indices
times = repmat(NaN, [ntim 11]);

glob = [0 360 -90 90];

ct = [180 270 -6 6];
nino3 = [210 270 -5 5];
nino4 = [160 210 -5 5];
nino3p4 = [190 240 -5 5];
nino1p2 = [270 280 -10 0];

lims2 = [glob; ct; nino3; nino4; nino3p4; nino1p2];

for i = 1:5;
  [xk, yk] = keep_var(lims2(i,:), lon, lat);
  tem = find(~isnan(lm1(yk,xk)));
  tem2 = sst(:,yk,xk);
  
  times(:,i) = mean(sst(:,tem)')';
end

tahiti = [209 211 -18 -16];
darwin = [129 131 -13 -12];
np = [160 220 30 65];

soidat = [0 360 -20 20];
soistar1 = [180 285 -20 20];