Documentation of storm_1pt_corr_map

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

Help text

  Look at correlation maps

Cross-Reference Information

Listing of script storm_1pt_corr_map

clean
cd /home/disk/tao/jyin/ccm3/matlab
load storms.mat
[b, a] = butter(9, 2/10);
hpsl31 = psl_inst_T31_5 - filtfilt(b, a, psl_inst_T31_5);
hpsl42 = psl_inst_T42_5 - filtfilt(b, a, psl_inst_T42_5);
hz50031 = z500_inst_T31_5 - filtfilt(b, a, z500_inst_T31_5);
hz50042 = z500_inst_T42_5 - filtfilt(b, a, z500_inst_T42_5);
x31 = find(XAX_T31==180.);
x42 = find(XAX_T42==180.);
y31 = find(YAX_T31 >= 46. & YAX_T31 <= 47);
y42 = find(YAX_T42 >= 46. & YAX_T42 <= 47);
tim1 = squeeze(hpsl31(:,y31,x31));
tim2 = squeeze(hpsl42(:,y42,x42));
lags = -2:2;
[reg1, c1] = regress_eof(hpsl31, tim1, lags);
[reg2, c2] = regress_eof(hz50031, tim1, lags);
[reg1a, c1a] = regress_eof(hpsl42, tim2, lags);
[reg2a, c2a] = regress_eof(hz50042, tim2, lags);
get_global; FRAME = [90 270 10 90];

XAX = XAX_T31; YAX = YAX_T31;
figure(1); fo(1); clf;
for i = 1:5;
  subplot(5,2,2*i-1);
    mcont3(c1(i,:,:), 0.2, 'stereo', [90 180]);
    gridm on;
    axis_limits;
  subplot(5,2,2*i);
    mcont3(c2(i,:,:), 0.2, 'stereo', [90 180]);
    gridm on;
    axis_limits;
end
subplot(5,2,1);
  title('T31 corr(SLP, SLP(45^\circ,180^\circ)')
subplot(5,2,2);
  title('T31 corr(Z500, SLP(45^\circ,180^\circ)')
subplot(5,2,9);
  xlabel('Contour:  0.2')
subplot(5,2,10);
  xlabel('Contour:  0.2')

cd ~/matlab/CCM/Jeff/Figs
print -dps2 T31_1pt_corr_maps.ps


XAX = XAX_T42; YAX = YAX_T42;
figure(2); fo(1); clf;
for i = 1:5;
  subplot(5,2,2*i-1);
    mcont3(c1a(i,:,:), 0.2, 'stereo', [90 180]);
    gridm on;
    axis_limits;
  subplot(5,2,2*i);
    mcont3(c2a(i,:,:), 0.2, 'stereo', [90 180]);
    gridm on;
    axis_limits;
end
subplot(5,2,1);
  title('T42 corr(SLP, SLP(45^\circ,180^\circ)')
subplot(5,2,2);
  title('T42 corr(Z500, SLP(45^\circ,180^\circ)')
subplot(5,2,9);
  xlabel('Contour:  0.2')
subplot(5,2,10);
  xlabel('Contour:  0.2')

cd ~/matlab/CCM/Jeff/Figs
print -dps2 T42_1pt_corr_maps.ps


%  Look at regression maps

XAX = XAX_T31; YAX = YAX_T31;
figure(1); fo(1); clf;
for i = 1:5;
  subplot(5,2,2*i-1);
    mcont3(0.01*reg1(i,:,:), 1, 'stereo', [90 180]);
    gridm on;
    axis_limits;
  subplot(5,2,2*i);
    mcont3(reg2(i,:,:), 7.5, 'stereo', [90 180]);
    gridm on;
    axis_limits;
end
subplot(5,2,1);
  title('T31 reg(SLP, SLP(45^\circ,180^\circ)')
subplot(5,2,2);
  title('T31 reg(Z500, SLP(45^\circ,180^\circ)')
subplot(5,2,9);
  xlabel('Contour:  1mb')
subplot(5,2,10);
  xlabel('Contour:  7.5m')

cd ~/matlab/CCM/Jeff/Figs
print -dps2 T31_1pt_reg_maps.ps


XAX = XAX_T42; YAX = YAX_T42;
figure(2); fo(1); clf;
for i = 1:5;
  subplot(5,2,2*i-1);
    mcont3(0.01*reg1a(i,:,:), 1, 'stereo', [90 180]);
    gridm on;
    axis_limits;
  subplot(5,2,2*i);
    mcont3(reg2a(i,:,:), 7.5, 'stereo', [90 180]);
    gridm on;
    axis_limits;
end
subplot(5,2,1);
  title('T42 reg(SLP, SLP(45^\circ,180^\circ)')
subplot(5,2,2);
  title('T42 reg(Z500, SLP(45^\circ,180^\circ)')
subplot(5,2,9);
  xlabel('Contour:  1mb')
subplot(5,2,10);
  xlabel('Contour:  7.5m')

cd ~/matlab/CCM/Jeff/Figs
print -dps2 T42_1pt_reg_maps.ps