Documentation of compare_som_fixed

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

Help text

  SLP:

Cross-Reference Information

Listing of script compare_som_fixed

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
filin1 = ['time_wgr_hard.nc'];
filin2 = ['time_clim_hard.nc'];
filin3 = ['time_wgr_5093.nc'];
filin4 = ['time_clim_5093.nc'];
fignum = 1;
filin1 = ['time_cgr_hard.nc'];
filin2 = ['time_clim_hard.nc'];
filin3 = ['time_cgr_5093.nc'];
filin4 = ['time_clim_5093.nc'];
fignum = 2;
filin1 = ['time_wgr_hard.nc'];
filin2 = ['time_cgr_hard.nc'];
filin3 = ['time_wgr_5093.nc'];
filin4 = ['time_cgr_5093.nc'];
[lat, lon, level] = getll('time_cgr_5093.nc');
get_global
default_global
nam_var = 'TS'; lev = 'PSL';
[var1w ts1w] = getnc(filin1, lev, nam_var);
[var1c ts1c] = getnc(filin2, lev, nam_var);
[var2w ts2w] = getnc(filin3, lev, nam_var);
[var2c ts2c] = getnc(filin4, lev, nam_var);
nam_var = 'Z3'; lev = 8;
var1w = getnc(filin1, lev, nam_var);
var1c = getnc(filin2, lev, nam_var);
var2w = getnc(filin3, lev, nam_var);
var2c = getnc(filin4, lev, nam_var);
nam_var = 'PRECL'; lev = 'PRECC';
[a, b] = getnc(filin1, lev, nam_var);
var1w = (a + b) * (3600 * 24 * 1000);
[a, b] = getnc(filin2, lev, nam_var);
var1c = (a + b) * (3600 * 24 * 1000);
[a, b] = getnc(filin3, lev, nam_var);
var2w = (a + b) * (3600 * 24 * 1000);
[a, b] = getnc(filin4, lev, nam_var);
var2c = (a + b) * (3600 * 24 * 1000);
[ntim, nlat, nlon] = size(var1w);
win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);
ns = -1; cint = 0.5; proj = 'stereo';
tit1 = 'SLP:';  tit2 = ' mb';
lims = [ns*90 270];

%  HEIGHT:
ns = 1; cint = 10; proj = 'stereo';
tit1 = '250 mb Height:';  tit2 = ' m';
lims = [ns*90 270];

%  PRECIP:
ns = 1; cint = 1; proj = 'giso';
tit1 = 'PRECIPITATION:';  tit2 = ' mm (day)^-^1';
lims = [ns*0 180]; FRAME = [0 360 -40 40];

[tstat1, dof1] = ttest(var1w(win,:,:), var1c(win,:,:));
score1 = tscore(dof1, 2.5);
[tstat2, dof2] = ttest(var2w(win,:,:), var2c(win,:,:));
score2 = tscore(dof2, 2.5);

top = squeeze(mean(var1w(win,:,:) - var1c(win,:,:))) / 1;
bot = squeeze(mean(var2w(win,:,:) - var2c(win,:,:))) / 1;

maxlev = max([max(max(top)) max(max(bot))]);
minlev = min([min(min(top)) min(min(bot))]);
clev = sort(unique([(-0.5*cint):(-1*cint):minlev (0.5*cint):cint:maxlev]));

figure(fignum)
sp(1)
     clma
     mcont(top, clev, proj, lims);
     greyshd(abs(tstat1), score1);
     title([tit1 '  SOM\_HARD (WARM - CLIM)']);
     xlabel(['Contour Interval:  ' num2str(cint) tit2])

sp(2)
     clma
     mcont(bot, clev, proj, lims);
     greyshd(abs(tstat2), score2);
     title([tit1 '  DOM (WARM - CLIM)']);
     xlabel(['Contour Interval:  ' num2str(cint) tit2])

cd /home/disk/tao/dvimont/matlab/CCM/SOM/SOM_Plots



%  TS:
ns = 1; cint = 0.25; proj = 'giso';
tit1 = 'Surface Temperature:';  tit2 = ' C';
lims = [ns*0 180]; FRAME = [0 360 -90 90];
[tstat1, dof1] = ttest(ts1w(win,:,:), ts1c(win,:,:));
score1 = tscore(dof1, 2.5);
top = squeeze(mean(ts1w(win,:,:) - ts1c(win,:,:)))/2;
maxlev = max([max(max(top)) max(max(bot))]);
minlev = min([min(min(top)) min(min(bot))]);
clev = sort(unique([(-0.5*cint):(-1*cint):minlev (0.5*cint):cint:maxlev]));

figure(2);  orient landscape;
     clma
     mcont(top, clev, proj, lims);
     greyshd(abs(tstat2), score2);
     title([tit1 '  SOM\_HARD (WARM - COLD)/2']);
     xlabel(['Contour Interval:  ' num2str(cint) tit2])

cd /home/disk/tao/dvimont/matlab/CCM/SOM/SOM_Plots
% print -dps TS_DOM_SOMH_win.ps                     



%  Winds:

nam_var = 'U'; nam_var2 = 'V'; lev = 18;
[var1w ts1w] = getnc(filin1, lev, nam_var, nam_var2);
[var1c ts1c] = getnc(filin2, lev, nam_var, nam_var2);
[var2w ts2w] = getnc(filin3, lev, nam_var, nam_var2);
[var2c ts2c] = getnc(filin4, lev, nam_var, nam_var2);

var1w = squeeze(var1w); var2w = squeeze(var2w); 
var1c = squeeze(var1c); var2c = squeeze(var2c); 

ns = 1; cint = 1; proj = 'giso';
tit1 = 'WINDS:';  tit2 = [];
lims = [ns*0 180]; FRAME = [0 360 -40 40];

c1 = var1w + j*ts1w; c1 = reshape(c1, ntim, nlat*nlon);
c2 = var1c + j*ts1c; c2 = reshape(c2, ntim, nlat*nlon);
[fval, dof1, dof2] = vectsig(c1(win,:), c2(win,:));
fval = reshape(fval, nlat, nlon);
top = squeeze(mean(c1(win,:) - c2(win,:))) / 2;
top = reshape(top, nlat, nlon);
figure(1)
sp(1)
     clma;
     mquiv(real(top), imag(top), 1.5, proj, lims);
     greyshd(abs(fval), 3.17);
     title([tit1 '  SOM\_HARD (WARM - COLD)/2']);
%     xlabel(['Contour Interval:  ' num2str(cint) tit2])

c1 = var2w + j*ts2w; c1 = reshape(c1, ntim, nlat*nlon);
c2 = var2c + j*ts2c; c2 = reshape(c2, ntim, nlat*nlon);
[fval, dof1, dof2] = vectsig(c1(win,:), c2(win,:));
fval = reshape(fval, nlat, nlon);
top = squeeze(mean(c1(win,:) - c2(win,:))) / 2;
top = reshape(top, nlat, nlon);
sp(2)
     clma;
     mquiv(real(top), imag(top), 1.5, proj, lims);
     greyshd(abs(fval), 3.17);
     title([tit1 '  DOM (WARM - CLIM)']);
%     xlabel(['Contour Interval:  ' num2str(cint) tit2])