Global Index (short | long) | Local contents | Local Index (short | long)
Get ATL CTI
| This script calls | |
|---|---|
clean
ncep
cti = getnc2('skt.mon.mean.nc', 'skt', [340 360 -3 3], 1, 1:648);
cti = squeeze(mean(mean(shiftdim(cti, 1))));
cti = annave(cti);
% Define regions, etc
lims = [0 360 -90 90];
tim = 13:648;
yr = [(1949+1/24):1/12:(2001.99)];
% Get time series
ct = cti(tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
limsatl = [-35.625 5.625 -4.76 4.76];
limspac = [-36 6 -5 5];
% Now, load data for regressions
ncep
[skt, lat, lon] = getnc2('skt.mon.mean.nc', 'skt', lims, 1, tim);
[skt, clim] = annave(skt);
skt = detrend(skt);
vpac = detrend(vatl);
ct = detrend(ct);
% Round 1:
timeseries = vpac;
data = skt;
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(1); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]*0.5);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, SKT', 'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
plot(yr, standardize(timeseries), '-k', yr, standardize(ct), '--r');
axis([1948 2003 -4 4]);
grid on;
legend('VATL1', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_corrmap.ps
% Round 2:
timeseries = lin_remove(vpac, ct);
data = lin_remove(skt, ct);
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(2); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]*0.5);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, SKT; CT\_ATL linearly removed', ...
'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
plot(yr, standardize(timeseries), '-k', yr, standardize(ct), '--r');
axis([1948 2003 -4 4]);
grid on;
legend('VATL1', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_corrmap_noct.ps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Now, do the same for OLR:
clean
% Get ATL CTI
ncep
cti = getnc2('skt.mon.mean.nc', 'skt', [340 360 -3 3], 1, 1:648);
cti = squeeze(mean(mean(shiftdim(cti, 1))));
cti = annave(cti);
% Define regions, etc
lims = [0 360 -90 90];
tim = (12*(1974-1948)+6):(12*(2001-1948)+12);
yr = [(1974+11/24):1/12:(2001.99)];
tim2 = 1:331;
% Get time series
ct = cti(tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
limsatl = [-35.625 5.625 -4.76 4.76];
limspac = [-36 6 -5 5];
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
vatl = vatl(tim-12);
% Now, load data for regressions
cd /home/disk/tao/data/olr/interpolated
[olr, lat, lon] = getnc2('olrint.mon.mean.nc', 'olr', lims, 1, tim2);
[olr, clim] = annave(olr);
olr = detrend_NaN(olr);
vpac = detrend(vatl);
ct = detrend(ct);
kp = find(~isnan(olr(:,1,1)));
olr = olr(kp,:,:);
ct = ct(kp);
vpac = vpac(kp);
% Round 1:
timeseries = vpac;
data = olr;
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(1); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, OLR', 'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
tem = NaN*yr; tem(kp) = yr(kp);
tem1 = NaN*yr; tem1(kp) = standardize(vpac);
tem2 = NaN*yr; tem2(kp) = standardize(ct);
plot(tem, tem1, '-k', tem, tem2, '--r');
axis([1973 2003 -4 4]);
grid on;
legend('VATL1', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_olr_corrmap.ps
% Round 2:
timeseries = lin_remove(vpac, ct);
data = lin_remove(olr, ct);
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(2); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, OLR; CT\_ATL linearly removed', ...
'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
tem = NaN*yr; tem(kp) = yr(kp);
tem1 = NaN*yr; tem1(kp) = standardize(timeseries);
tem2 = NaN*yr; tem2(kp) = standardize(ct);
plot(tem, tem1, '-k', tem, tem2, '--r');
axis([1973 2003 -4 4]);
grid on;
legend('VATL2', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_olr_corrmap_noct.ps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Now, do the same for CMAP:
clean
% Get ATL CTI
ncep
cti = getnc2('skt.mon.mean.nc', 'skt', [340 360 -3 3], 1, 1:648);
cti = squeeze(mean(mean(shiftdim(cti, 1))));
cti = annave(cti);
% Define regions, etc
lims = [0 360 -90 90];
tim = (12*(1979-1948)+1):(12*(1999-1948)+12);
yr = [(1979+1/24):1/12:(1999.99)];
tim2 = 1:252;
% Get time series
ct = cti(tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
limsatl = [-35.625 5.625 -4.76 4.76];
limspac = [-36 6 -5 5];
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
vatl = vatl(tim-12);
% Now, load data for regressions
cd /home/disk/tao/data/merged_precip/cmap
[olr, lat, lon] = ...
getnc2('precipcmap.monthly.7999.nc', 'data', lims, 1, tim2);
[olr, clim] = annave(olr);
olr = detrend(olr);
vpac = detrend(vatl);
ct = detrend(ct);
% Round 1:
timeseries = vpac;
data = olr;
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(1); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, CMAP', 'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
tem = NaN*yr; tem = yr;
tem1 = NaN*yr; tem1 = standardize(timeseries);
tem2 = NaN*yr; tem2 = standardize(ct);
plot(tem, tem1, '-k', tem, tem2, '--r');
axis([1978 2001 -4 4]);
grid on;
legend('VATL1', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_cmap_corrmap.ps
% Round 2:
timeseries = lin_remove(vpac, ct);
data = lin_remove(olr, ct);
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(2); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso', [0 -45]);
map_contour_pn(cpat, 0.2, 'zeroline');
hold on;
hh = map_surface(cpat, -1);
hold off;
caxis([-1 1]);
draw_landmap;
tightmap; gridm on;
setm(gca, 'plinelocation', 30, 'mlinelocation', 45, 'frame', 'off');
colorbar2;
t2 = title2('Correlation: VATL1, CMAP; CT\_ATL linearly removed', ...
'fontsize', 10);
bx = drawboxm(limspac);
set(bx, 'linewidth', 2);
global_axes(6.5, 1.5, 0.5, 0.5, 5.5+.5);
subplot2(1,1);
tem = yr; tem1 = standardize(timeseries); tem2 = standardize(ct);
plot(tem, tem1, '-k', tem, tem2, '--r');
axis([1978 2001 -4 4]);
grid on;
legend('VATL2', 'CT\_ATL');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vatl1_cmap_corrmap_noct.ps