Global Index (short | long) | Local contents | Local Index (short | long)
Define regions, etc
| This script calls | |
|---|---|
clean
cd /home/disk/tao/dvimont/Papers/obs/matlab/Data
load cti2.mat;
lims = [0 360 -90 90];
tim = 13:648;
yr = [(1949+1/24):1/12:(2001.99)];
% Get time series
ct = ct2(tim);
ncep
limsatl = [360-35.625 5.625 -4.76 4.76];
%limspac = [180 138.75 -4.76 4.76];
limspac = [180 222 -5 5];
[vwnd1, lat1, lon1] = ...
getnc2('vwnd10m.mon.mean.nc', 'vwnd', limspac, 1, tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
[ntim, nlat, nlon] = size(vwnd1);
tem = ones(2,nlat,nlon);
tem = cosweight(tem, lat1);
tem = cosweight(tem, lat1);
tem = sum(sum(squeeze(tem(1,:,:))));
vwnd1 = cosweight(vwnd1, lat1);
vwnd1 = cosweight(vwnd1, lat1);
vpac2 = squeeze(sum(sum(shiftdim(vwnd1, 1))))./tem;
[vpac2, vclim] = annave(vpac2);
% 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(vpac);
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');
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: VPAC1, SKT', 'fontsize', 10);
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([1945 2005 -4 4]);
grid on;
legend('VPAC1', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_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');
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: VPAC1, SKT; CTI 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([1945 2005 -4 4]);
grid on;
legend('VPAC1', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_corrmap_noct.ps
% Round 3: Get rid of outliers
kp = find(abs(standardize(vpac)) < 3);
timeseries = lin_remove(vpac(kp), ct(kp));
data = lin_remove(skt(kp,:,:), ct(kp));
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(3); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso');
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: VPAC1, SKT; CTI linearly removed; ' ...
'Outliers 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);
tem1 = NaN*yr; tem1(kp) = yr(kp);
tem2 = NaN*vpac; tem2(kp) = standardize(vpac(kp));
tem3 = NaN*ct; tem3(kp) = standardize(ct(kp));
plot(tem1, tem2, '-k', ...
tem1, tem3, '--r');
axis([1945 2005 -4 4]);
grid on;
legend('VPAC3', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_corrmap_noct2.ps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Now, do the same for OLR:
clean
cd /home/disk/tao/dvimont/Papers/obs/matlab/Data
load cti2.mat;
% 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;
limspac = [180 222 -5 5];
% Get time series
ct = ct2(tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
vpac = vpac(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(vpac);
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');
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: VPAC1, OLR', 'fontsize', 10);
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([1970 2005 -4 4]);
grid on;
legend('VPAC1', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_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');
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: VPAC1, OLR; CTI 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([1970 2005 -4 4]);
grid on;
legend('VPAC2', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_olr_corrmap_noct.ps
% Round 3: Get rid of outliers
kp2 = find(abs(standardize(vpac)) < 3);
timeseries = lin_remove(vpac(kp2), ct(kp2));
data = lin_remove(olr(kp2,:,:), ct(kp2));
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(3); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso');
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: VPAC1, OLR; CTI linearly removed; ' ...
'Outliers 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);
kp3 = kp(kp2);
tem1 = NaN*yr; tem1(kp3) = yr(kp3);
tem2 = NaN*yr; tem2(kp3) = standardize(vpac(kp2));
tem3 = NaN*yr; tem3(kp3) = standardize(ct(kp2));
plot(tem1, tem2, '-k', ...
tem1, tem3, '--r');
axis([1970 2005 -4 4]);
grid on;
legend('VPAC3', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_olr_corrmap_noct2.ps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Now, do the same for CMAP:
clean
cd /home/disk/tao/dvimont/Papers/obs/matlab/Data
load cti2.mat;
% 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;
limspac = [180 222 -5 5];
% Get time series
ct = ct2(tim);
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Data
load V_timeseries.mat
vpac = vpac(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(vpac);
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');
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: VPAC1, CMAP', 'fontsize', 10);
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('VPAC1', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_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');
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: VPAC1, CMAP; CTI 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('VPAC2', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_cmap_corrmap_noct.ps
% Round 3: Get rid of outliers
kp = find(abs(standardize(vpac)) < 3);
timeseries = lin_remove(vpac(kp), ct(kp));
data = lin_remove(olr(kp,:,:), ct(kp));
[rpat, cpat] = regress_eof(data, timeseries);
figure_tall(3); clf;
global_axes(6.5, 4, 0.5, 0.5, 1.5);
global_latlon(lat, lon, lims);
subplot2(1,1);
map_axis('giso');
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: VPAC1, CMAP; CTI linearly removed; ' ...
'Outliers 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);
tem1 = NaN*yr; tem1(kp) = yr(kp);
tem2 = NaN*vpac; tem2(kp) = standardize(vpac(kp));
tem3 = NaN*ct; tem3(kp) = standardize(ct(kp));
plot(tem1, tem2, '-k', ...
tem1, tem3, '--r');
axis([1978 2001 -4 4]);
grid on;
legend('VPAC3', 'CTI');
cd /home/disk/tao/dvimont/matlab/NMC/Atlantic/Plots
print -painters -dpsc2 vpac1_cmap_corrmap_noct2.ps