Documentation of ttest_plot

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

Help text

win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);

Cross-Reference Information

This script calls
clf default_global frame get_global getll getnc greyshd mcont mquiv sp tscore ttest vectsig

Listing of script ttest_plot


clear
if 0;
varina = 'U'; varinb = 'V'; filin = 'U_V_5093.nc'; ind = 18;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
  [var1a, var1b] = getnc(filin, ind, varina, varinb);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
  [var2a, var2b] = getnc(filin, ind, varina, varinb);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data
  [var3a, var3b] = getnc(filin, ind, varina, varinb);
  [lat, lon, level] = getll(filin);
end
if 0;
varina = 'FSNS'; varinb = 'FLNS'; varinc = 'LHFLX'; varind = 'SHFLX';
filin = 'FLUX_5093.nc';
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
  [var1a, var1b, var1c, var1d] = ...
     getnc(filin, varina, varinb, varinc, varind);
  ts1 = getnc('wgr_5093.nc', 'TS');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
  [var2a, var2b, var2c, var2d] = ...
     getnc(filin, varina, varinb, varinc, varind);
  ts2 = getnc('cgr_5093.nc', 'TS');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data
  [var3a, var3b, var3c, var3d] = ...
     getnc(filin, varina, varinb, varinc, varind);
  [ts3, oro] = getnc('clim_5093.nc', 'TS', 'ORO');
[lat, lon, level] = getll(filin);
end
if 0;
varina = 'Z3';
varina = 'U';
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
[var1] = getnc([varina '_5093.nc'], varina);
ps1 = getnc('PS_5093.nc', 'PS');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
[var2] = getnc([varina '_5093.nc'], varina);
ps2 = getnc('PS_5093.nc', 'PS');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data
[var3] = getnc([varina '_5093.nc'], varina);
ps3 = getnc('PS_5093.nc', 'PS');
[lat, lon, level] = getll([varina '_5093.nc']);
end;
if 0;
varina = 'PRECC';
varinb = 'PRECL';
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
[var1a, var1b] = getnc([varina '_5093.nc'], varina, varinb);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
[var2a, var2b] = getnc([varina '_5093.nc'], varina, varinb);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data
[var3a, var3b] = getnc([varina '_5093.nc'], varina, varinb);
[lat, lon, lev] = getll([varina '_5093.nc']);
var1 = (var1a + var1b) * (1000 * 3600 * 24);
var2 = (var2a + var2b) * (1000 * 3600 * 24);
var3 = (var3a + var3b) * (1000 * 3600 * 24);
end
[ntim, nlev, nlat, nlon] = size(var1);
get_global
default_global
win = sort([8:12:ntim 9:12:ntim 10:12:ntim]);

clf

ns = -1;
lev = 6;
cint = 10;
var1win = squeeze(var1(win,lev,:,:));
var2win = squeeze(var2(win,lev,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(1);
  clma
  mcont(top, [-105:cint:110], [ns*90 270]);
  greyshd(abs(tstat), score);
  title('250mb Z3:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' m']);

lev = 11;
cint = 10;
var1win = squeeze(var1(win,lev,:,:));
var2win = squeeze(var2(win,lev,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(2);
  clma;
  mcont(top, [-105:cint:110], [ns*90 270]);
  greyshd(abs(tstat), score);
  title('500mb U:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' m']);

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots

%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2D variables %%%%%%%%%%%%%%%%%%%%%%%%%%

get_global
default_global
FRAME = [0 360 -60 60];
figure(1)
[ntim, nlat, nlon] = size(var1);
win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);

var1win = squeeze(var1(win,:,:));
var2win = squeeze(var3(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(1);
  clma;
  mcont(top, [-5.5:1:5.5], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF PRECIPITATION:  (WGR\_5093 - CLIM\_5093) / 2');
  xlabel('Contour Interval:  1 mm / day');

var1win = squeeze(var2(win,:,:));
var2win = squeeze(var3(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(2);
  clma;
  mcont(top, [-5.5:1:5.5], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF PRECIPITATION:  (CGR\_5093 - CLIM\_5093) / 2');
  xlabel('Contour Interval:  1 mm / day');

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots





%%%%%%%%%%%%%%%%%%%%%%%%% FLUX %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%




[ntim, nlev, nlat, nlon] = size(var1a);
get_global
default_global
FRAME = [0 360 -40 40];
%win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);
%win = sort([8:12:ntim 9:12:ntim 10:12:ntim]);
%win = 1:ntim;

var1win = squeeze(var1a(win,:,:));
var2win = squeeze(var2a(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = -1*squeeze(mean(var1win - var2win))/2;
sp(1);
  clma;
  mcont(top, [-52:8:52], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF FSNS:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  8 W / m^2  (Positive Upwards)');

var1win = squeeze(var1b(win,:,:));
var2win = squeeze(var2b(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(2);
  clma;
  mcont(top, [-50:4:50], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF FLNS:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  4 W / m^2  (Positive Upwards)');

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots

%%%%%%%%% LHFLX and SHFLX %%%%%%%%%%%

var1win = squeeze(var1c(win,:,:));
var2win = squeeze(var2c(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = 1*squeeze(mean(var1win - var2win))/2;
sp(1);
  clma;
  mcont(top, [-52:8:52], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF LHFLX:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  8 W / m^2  (Positive Upwards)');

var1win = squeeze(var1d(win,:,:));
var2win = squeeze(var2d(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
sp(2);
  clma;
  mcont(top, [-51:2:51], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF LHFLX:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  2 W / m^2  (Positive Upwards)');

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots

%%%%%%%%% Total Head Flux %%%%%%%%%%%

var1win = -1 * (var1a - var1b - var1c - var1d);
var1win = var1win(win,:,:);
var2win = -1 * (var2a - var2b - var2c - var2d);
var2win = var2win(win,:,:);
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = 1*squeeze(mean(var1win - var2win))/2;
sp(1);
  clma;
  mcont(top, [-55:10:55], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF Total Head Flux:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  10 W / m^2  (Positive Upwards)');

orom = squeeze(mean(oro));
lm = -1 * (orom - 1);
top = squeeze(mean(ts1 - ts2))/2;
top = top .* lm;
sp(2);
  clma;
  mcont(top, [-1:.1:1], 'giso');
  title('SST:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  0.1 K');

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots



%%%%%%%%%%%%%%%%%%%%%% Winds %%%%%%%%%%%%%%%%%%%%%%%%%%%%


get_global
default_global
FRAME = [0 360 -40 40];

for i = 1:3;
  eval(['var' num2str(i) 'a = squeeze(var' num2str(i) 'a);']);
  eval(['var' num2str(i) 'b = squeeze(var' num2str(i) 'b);']);
end;

wspd1 = sqrt(var1a.^2 + var1b.^2);
wspd2 = sqrt(var2a.^2 + var2b.^2);
wspd3 = sqrt(var3a.^2 + var3b.^2);

[ntim, nlat, nlon] = size(var1a);
win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);
%win = sort([8:12:ntim 9:12:ntim 10:12:ntim]);
%win = 1:ntim;

var1win = squeeze(wspd2(win,:,:));
var2win = squeeze(wspd3(win,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = 1*squeeze(mean(var1win - var2win))/1;
sp(1);
  clma;
  mcont(top, [-2.25:.5:2.25], 'giso', [0 180]);
  greyshd(abs(tstat), score);
  title('DJF WIND SPEED:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Contour Interval:  0.5 m / s');

%  Calculate significance of wind vector:

u = var2a(win,:,:);
u = reshape(u, size(u, 1), nlat*nlon);
v = var2b(win,:,:);
v = reshape(v, size(v, 1), nlat*nlon);
c1 = u + j*v;
u = var3a(win,:,:);
u = reshape(u, size(u, 1), nlat*nlon);
v = var3b(win,:,:);
v = reshape(v, size(v, 1), nlat*nlon);
c2 = u + j*v;

[fval, dof1, dof2] = vectsig(c1, c2);
fval = reshape(fval, nlat, nlon);

%  Get f-score from table.
fscore = 3.23;

top = squeeze(mean(c1-c2))/1;
top = reshape(top, nlat, nlon);
sp(2);
  clma;
  mquiv(real(top), imag(top), 2, 'giso', [0 180]);
  greyshd(fval, fscore);
  title('DJF WINDS:  (WGR\_5093 - CGR\_5093) / 2');
  xlabel('Shading indicates vectors that exceed the 95% significance level');

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots