Documentation of make_time_plots

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

Help text




  Load variables directly from netcdf files

Cross-Reference Information

This script calls
dc dc2 dcm2 default_global frame get_global getll getnc greyshd mcont mquiv sp tscore ttest vectsig

Listing of script make_time_plots


clear
lab = 'U';
lab = 'V';
lab = 'Z3';
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
eval(['load ' lab '_som_dom.mat']);
get_global; default_global;
eval(['[ntim, nlev, nlat, nlon] = size(' lab 'dw);']);
tim = sort([1:12:120 2:12:120 3:12:120 4:12:120 5:12:120 6:12:120]);
lab = 'U';
Udw = sqrt(Udw.^2 + Vdw.^2);
Udc = sqrt(Udc.^2 + Vdc.^2);
Usw = sqrt(Usw.^2 + Vsw.^2);
Usc = sqrt(Usc.^2 + Vsc.^2);
hgt = 850;
lind = find(lev2 == hgt);
i = 1;  % i = 1 ==> W-C, i=2 ==> S-D
if i == 1;
  eval(['var1 = squeeze(' lab 'dw(tim,lind,:,:));']);
  eval(['var2 = squeeze(' lab 'dc(tim,lind,:,:));']);
elseif i == 2;
  eval(['var1 = squeeze(' lab 'sw(tim,lind,:,:));']);
  eval(['var2 = squeeze(' lab 'dw(tim,lind,:,:));']);
end
[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
cint = 0.25; clev = [-10:cint:10];
sp(1);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' num2str(hgt) 'mb ' lab ':  (WGR\_DOM - CGR\_DOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' num2str(hgt) 'mb ' lab ':  (WGR\_SOM - WGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' m s^-^1']);
if i == 1;
  eval(['var1 = squeeze(' lab 'sw(tim,lind,:,:));']);
  eval(['var2 = squeeze(' lab 'sc(tim,lind,:,:));']);
elseif i == 2;
  eval(['var1 = squeeze(' lab 'sc(tim,lind,:,:));']);
  eval(['var2 = squeeze(' lab 'dc(tim,lind,:,:));']);
end
[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
sp(2);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' num2str(hgt) 'mb ' lab ':  (WGR\_SOM - CGR\_SOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' num2str(hgt) 'mb ' lab ':  (CGR\_SOM - CGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' m s^-^1']);
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM

clear

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data

varid1 = 'U'; varid2 = 'V';
varid1 = 'TS'; varid2 = 'PSL';
varid1 = 'PRECC'; varid2 = 'PRECL';
[dw1, dw2] = getnc('time_wgr_5093.nc', 18, varid1, varid2);
[dc1, dc2] = getnc('time_cgr_5093.nc', 18, varid1, varid2);
[sw1, sw2] = getnc('time_wgr_hard.nc', 18, varid1, varid2);
[sc1, sc2] = getnc('time_cgr_hard.nc', 18, varid1, varid2);
[lat, lon, lev, hyam, hybm, P0] = getll('wgr_SOM_hard.nc');
[ntim, nlev, nlat, nlon] = size(dw1);
get_global; default_global;

%  PRECIP
dw = (dw1 + dw2) * (3600*24*1000);
dc = (dc1 + dc2) * (3600*24*1000);
sw = (sw1 + sw2) * (3600*24*1000);
sc = (sc1 + sc2) * (3600*24*1000);
FRAME = [0 360 -50 50];
FRAME = [30 120 -45 30];

% WINDSPEED
dw = squeeze(sqrt(dw1.^2 + dw2.^2));
dc = squeeze(sqrt(dc1.^2 + dc2.^2));
sw = squeeze(sqrt(sw1.^2 + sw2.^2));
sc = squeeze(sqrt(sc1.^2 + sc2.^2));

dw = squeeze(dw1);% / 100;
dc = squeeze(dc1);% / 100;
sw = squeeze(sw1);% / 100;
sc = squeeze(sc1);% / 100;

tim = sort([1:12:120 2:12:120 3:12:120 4:12:120 5:12:120 6:12:120]);

cint = 0.5; clev = [-5:cint:5]; lab = 'PRECIP';
cint = 0.1; clev = [-5:cint:5]; lab = 'TS';
cint = 0.5; clev = [-5:cint:5]; lab = 'SLP';
cint = 0.5; clev = [-5:cint:5]; lab = '970MB U';
cint = 2.5; clev = [-20:cint:20]; lab = 'Sigma * TS^4';
cint = 0.25; clev = [-5:cint:5]; lab = '970MB WINDSPEED';
i = 1;

if i == 1;
  var1 = squeeze(dw(tim,:,:));
  var2 = squeeze(dc(tim,:,:));
elseif i == 2;
  var1 = squeeze(sw(tim,:,:));
  var2 = squeeze(dw(tim,:,:));
end

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
sp(1);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_DOM - CGR\_DOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (WGR\_SOM - WGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' C']);% mm (day)^-^1']);

if i == 1;
  var1 = squeeze(sw(tim,:,:));
  var2 = squeeze(sc(tim,:,:));
elseif i == 2;
  var1 = squeeze(sc(tim,:,:));
  var2 = squeeze(dc(tim,:,:));
end

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
sp(2);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_SOM - CGR\_SOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (CGR\_SOM - CGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' C']);% mm (day)^-^1']);

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM
% print -dps2 






%%%%%%%%%%%%%%%%% For winds, do multivariate t-test %%%%%%%%%%%%%%%%%%%

lab = '990 mb WIND';

i = 1;

if i == 1;
  dw = squeeze(dw1(tim,:,:) + j*dw2(tim,:,:));
  dac = squeeze(dc1(tim,:,:) + j*dc2(tim,:,:));
  sw = squeeze(sw1(tim,:,:) + j*sw2(tim,:,:));
  sc = squeeze(sc1(tim,:,:) + j*sc2(tim,:,:));
elseif i == 2;
  dw = squeeze(sw1(tim,:,:) + j*sw2(tim,:,:));
  dac = squeeze(dw1(tim,:,:) + j*dw2(tim,:,:));
  sw = squeeze(sc1(tim,:,:) + j*sc2(tim,:,:));
  sc = squeeze(dc1(tim,:,:) + j*dc2(tim,:,:));
end

dw = reshape(dw, size(dw, 1), nlat*nlon);
dac = reshape(dac, size(dac, 1), nlat*nlon);
sw = reshape(sw, size(sw, 1), nlat*nlon);
sc = reshape(sc, size(sc, 1), nlat*nlon);
FRAME = [0 360 -60 60];

[fval, dof1, dof2] = vectsig(dw, dac);
fval = reshape(fval, nlat, nlon);
fscore = 3.07;
top = squeeze(mean(dw - dac)) / (3-i);
top = reshape(top, nlat, nlon);
sp(1);
  clma; 
  [h,a,c] = mquiv(real(top), imag(top), 1.2, 'giso', [0 180]);
  dcm2
%  gshade(fval, fscore);
  greyshd(fval, fscore);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_DOM - CGR\_DOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (WGR\_SOM - WGR\_DOM)']);
  end
  xlabel(['Maximum Windspeed:  ' num2str(c(2)) ' m s^-^1']);
%  xlabel(['Contour Interval:  ' num2str(cint) ' m s^-^1']);

[fval, dof1, dof2] = vectsig(sw, sc);
fval = reshape(fval, nlat, nlon);
fscore = 3.07;
top = squeeze(mean(sw - sc)) / (3-i);
top = reshape(top, nlat, nlon);
sp(2);
  clma; 
  [h,a,c] = mquiv(real(top), imag(top), 1.2, 'giso', [0 180]);
  dcm2
  greyshd(fval, fscore);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_SOM - CGR\_SOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (CGR\_SOM - CGR\_DOM)']);
  end
  xlabel(['Maximum Windspeed:  ' num2str(c(2)) ' m s^-^1']);
%  xlabel(['Contour Interval:  ' num2str(cint) ' m s^-^1']);

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