Documentation of make_time_plots

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

Help text

  Look at CT runs:

Cross-Reference Information

Listing of script make_time_plots

clear

if 0
%  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_5093/som_cycle/data
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_5093/cycle/data

  varn1 = 'Z3'; varn2 = 'AIRT'; varn3 = 'V';
  varn4 = 'U';% var5 = 'VV'; var6 = 'U';
  varn7 = 'PS'

  ind = 5:18;

  for i = 1:19;
    i
    if i < 10;
      filin = ['000' num2str(i) '.nc'];
    else
      filin = ['00' num2str(i) '.nc'];
    end

    nc = netcdf(filin, 'nowrite');
      v1 = squeeze(nc{varn1}(:,ind,:,:));
      v2 = squeeze(nc{varn2}(:,ind,:,:));
      v3 = squeeze(nc{varn3}(:,ind,:,:));
      v4 = squeeze(nc{varn4}(:,ind,:,:));
%      v5 = squeeze(nc{var5}(:,ind,:,:));
%      v6 = squeeze(nc{var6}(:,ind,:,:));
      v7 = squeeze(nc{varn7}(:,:,:));
    nc = close(nc);

    if i == 1;
      tind = [1:2];
    else
      tind = 12*(i-2) + 2 + [1:12];
    end

    var1(tind,:,:,:) = v1;
    var2(tind,:,:,:) = v2;
    var3(tind,:,:,:) = v3;
    var4(tind,:,:,:) = v4;
    ps(tind,:,:) = v7;

  end

  clear v1 v2 v3 v4 v5 v6 v7

  for i = 1:9;
    indc = 24*(i-1)+[1:6];
    indw = 24*(i-1)+[13:18];
    var1w(i,:,:,:) = mean(var1(indw,:,:,:));
    var1c(i,:,:,:) = mean(var1(indc,:,:,:));
    var2w(i,:,:,:) = mean(var2(indw,:,:,:));
    var2c(i,:,:,:) = mean(var2(indc,:,:,:));
    var3w(i,:,:,:) = mean(var3(indw,:,:,:));
    var3c(i,:,:,:) = mean(var3(indc,:,:,:));
    var4w(i,:,:,:) = mean(var4(indw,:,:,:));
    var4c(i,:,:,:) = mean(var4(indc,:,:,:));
    pw(i,:,:) = mean(ps(indw,:,:));
    pc(i,:,:) = mean(ps(indc,:,:));
  end

  clear var1 var2 var3 var4 ps

  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
  [lat, lon, lev, hyam, hybm, P0] = getll('wgr_SOM_hard.nc');
  hyam = hyam(ind); hybm = hybm(ind);
  get_global; default_global;

  newlev = [150 250 300 500 850 970 1000];
  z3w = atlev(var1w, newlev, pw, hyam, hybm, P0);
  z3c = atlev(var1c, newlev, pc, hyam, hybm, P0);
  airtw = atlev(var2w, newlev, pw, hyam, hybm, P0);
  airtc = atlev(var2c, newlev, pc, hyam, hybm, P0);
  vw = atlev(var3w, newlev, pw, hyam, hybm, P0);
  vc = atlev(var3c, newlev, pc, hyam, hybm, P0);
  uw = atlev(var4w, newlev, pw, hyam, hybm, P0);
  uc = atlev(var4c, newlev, pc, hyam, hybm, P0);
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
  %save uvzt_ctdom.mat z3w z3c airtw airtc vw vc uw uc newlev lat lon

end

if 1
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
  load uvzt_ctsom.mat
  z3w1 = z3w; z3c1 = z3c;
  airtw1 = airtw; airtc1 = airtc;
  vw1 = vw; vc1 = vc;
  uw1 = uw; uc1 = uc;
  load uvzt_ctdom.mat
end
get_global; default_global;

lev = 4;
varn = 'z3';
tit1 = ['NDJFMA CT\_SOM:  ' num2str(newlev(lev)) 'mb HEIGHT (WCT - CCT)/2'];
tit2 = ['NDJFMA CT\_DOM:  ' num2str(newlev(lev)) 'mb HEIGHT (WCT - CCT)/2'];
tit3 = ['NDJFMA CT:  ' num2str(newlev(lev)) 'mb HEIGHT'];
figure(1); FRAME = [0 360 -90 90];
cint = 10; clev = [-90:cint:90];
xl = ' m';

lev = 4;
varn = 'u';
tit1 = ['NDJFMA CT\_SOM:  ' num2str(newlev(lev)) 'mb U (WCT - CCT)/2'];
tit2 = ['NDJFMA CT\_DOM:  ' num2str(newlev(lev)) 'mb U (WCT - CCT)/2'];
tit3 = ['NDJFMA CT:  ' num2str(newlev(lev)) 'mb U'];
figure(1); FRAME = [0 360 -90 90];
cint = 1; clev = [-9:cint:9];
xl = ' m s^-^1';

eval(['var1 = squeeze(' varn 'w1(:,lev,:,:));']);
eval(['var2 = squeeze(' varn 'c1(:,lev,:,:));']);

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/2;

sp(1)
     if ismap(gca); clma; end;
     mcont(top, clev);
     greyshd(abs(tstat), score);
     title(tit1);
     xlabel(['Contour Interval:  ' num2str(cint) xl]);

eval(['var1 = squeeze(' varn 'w(:,lev,:,:));']);
eval(['var2 = squeeze(' varn 'c(:,lev,:,:));']);

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/2;

sp(2)
     if ismap(gca); clma; end;
     mcont(top, clev);
     greyshd(abs(tstat), score);
     title(tit2);
     xlabel(['Contour Interval:  ' num2str(cint) xl]);



figure(2); figure_orient;

i = 1;
i = 2;

eval(['var1 = squeeze(' varn 'w1(:,lev,:,:)) - squeeze(' varn 'c1(:,lev,:,:));']);
eval(['var2 = squeeze(' varn 'w(:,lev,:,:)) - squeeze(' varn 'c(:,lev,:,:));']);

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/2;

sp(i)
     if ismap(gca); clma; end;
     mcont(top, clev);
     greyshd(abs(tstat), score);
     title([tit3 ':  0.5*((WSOM-CSOM)-(WDOM-CDOM))'])
     xlabel(['Contour Interval:  ' num2str(cint) xl]);