Documentation of eddy_time_plots

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

Help text

save eddy_interped_grsom.mat vtw vtc kew kec newlev lat lon

Cross-Reference Information

This script calls
atlev default_global frame get_global getll mcont sp tscore ttest

Listing of script eddy_time_plots


clear
if 0;
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/som_5093/data
  var1 = 'VT'; var2 = 'AIRT'; var3 = 'V';
  var4 = 'UU'; var5 = 'VV'; var6 = 'U'; var7 = 'PS'
  ind = 5:18;
  nc = netcdf('EDDY_time_wgr.nc', 'nowrite');
    v1 = squeeze(nc{var1}(:,ind,:,:));
    v2 = squeeze(nc{var2}(:,ind,:,:));
    v3 = squeeze(nc{var3}(:,ind,:,:));
    v4 = squeeze(nc{var4}(:,ind,:,:));
    v5 = squeeze(nc{var5}(:,ind,:,:));
    v6 = squeeze(nc{var6}(:,ind,:,:));
    v7 = squeeze(nc{var7}(:,:,:));
  nc = close(nc);
  psw(:,:,:) = v7(:,:,:);
  e1w = v1-v2.*v3;
  e2w = v4+v5-v6.^2-v3.^2;
  clear v1 v2 v3 v4 v5 v6 v7
  nc = netcdf('EDDY_time_cgr.nc', 'nowrite');
    v1 = squeeze(nc{var1}(:,ind,:,:));
    v2 = squeeze(nc{var2}(:,ind,:,:));
    v3 = squeeze(nc{var3}(:,ind,:,:));
    v4 = squeeze(nc{var4}(:,ind,:,:));
    v5 = squeeze(nc{var5}(:,ind,:,:));
    v6 = squeeze(nc{var6}(:,ind,:,:));
    v7 = squeeze(nc{var7}(:,:,:));
  nc = close(nc);
  e1c = v1-v2.*v3;
  e2c = v4+v5-v6.^2-v3.^2;
  psc(:,:,:) = v7(:,:,:);
  clear v1 v2 v3 v4 v5 v6 v7
  for i = 1:10;
    ind = 12*(i-1)+[1:6];
    vtw(i,:,:,:) = mean(e1w(ind,:,:,:));
    vtc(i,:,:,:) = mean(e1c(ind,:,:,:));
    kew(i,:,:,:) = mean(e2w(ind,:,:,:));
    kec(i,:,:,:) = mean(e2c(ind,:,:,:));
    pw(i,:,:) = mean(psw(ind,:,:));
    pc(i,:,:) = mean(psc(ind,:,:));
  end
  ind = 5:18;
  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];
  vtw = atlev(vtw, newlev, pw, hyam, hybm, P0);
  vtc = atlev(vtc, newlev, pc, hyam, hybm, P0);
  kew = atlev(kew, newlev, pw, hyam, hybm, P0);
  kec = atlev(kec, newlev, pc, hyam, hybm, P0);
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
else
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
  load eddy_interped_grsom.mat
%  load eddy_interped_grdom.mat
end
get_global; default_global;

levt = 5; levb = 3;
[tstatt, dof] = ttest(vtw(:,levt,:,:), vtc(:,levt,:,:));
scoret = tscore(dof, 2.5);
top = squeeze(mean(vtw(:,levt,:,:) - vtc(:,levt,:,:)))/2;
cintt = 1; clevt = [-5:cintt:5];
[tstatb, dof] = ttest(kew(:,levb,:,:), kec(:,levb,:,:));
scoreb = tscore(dof, 2.5);
bot = squeeze(mean(kew(:,levb,:,:) - kec(:,levb,:,:)))/2;
cintb = 15; clevb = [-90:cintb:90];
tit1 = ['NDJFMAM GR Run:  ' num2str(newlev(levt)) 'mb v''T'', (WGR - CGR)/2'];
tit2 = ['NDJFMAM GR Run:  ' num2str(newlev(levb)) 'mb EKE, (WGR - CGR)/2'];
clabt = [' K m s^-^1']; clabb = [' m^2 s^-^2'];

lev = 5;
top = squeeze(mean(vtw(:,lev,:,:)));
bot = squeeze(mean(vtc(:,lev,:,:)));
cintt = 7.5; clevt = [-60:cintt:60]; cintb = cintt; clevb = clevt;
clabt = [' K m s^-^1']; clabb = clabt;
tit1 = ['WGR ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM'];
tit2 = ['CGR ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM'];

lev = 3;
top = squeeze(mean(kew(:,lev,:,:)));
bot = squeeze(mean(kec(:,lev,:,:)));
cintt = 100; clevt = [-600:cintt:600]; cintb = cintt; clevb = clevt;
clabt = [' K m s^-^1']; clabb = clabt;
tit1 = ['WGR ' num2str(newlev(lev)) 'mb EKE, NDJFMAM'];
tit2 = ['CGR ' num2str(newlev(lev)) 'mb EKE, NDJFMAM'];

figure(1); FRAME = [0 360 -90 90];

sp(1)
     if ismap(gca); clma; end;
     mcont(top, clevt);
     title(tit1); 
     xlabel(['Contour Interval:  ' num2str(cintt) clabt])
%     greyshd(abs(tstatt), scoret);
sp(2)
     if ismap(gca); clma; end;
     mcont(bot, clevb);
     title(tit2); 
     xlabel(['Contour Interval:  ' num2str(cintb) clabb])
%     greyshd(abs(tstatb), scoreb);





%  Look at CT runs:

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

  var1 = 'VT'; var2 = 'AIRT'; var3 = 'V';
  var4 = 'UU'; var5 = 'VV'; var6 = 'U'; var7 = '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{var1}(:,ind,:,:));
      v2 = squeeze(nc{var2}(:,ind,:,:));
      v3 = squeeze(nc{var3}(:,ind,:,:));
      v4 = squeeze(nc{var4}(:,ind,:,:));
      v5 = squeeze(nc{var5}(:,ind,:,:));
      v6 = squeeze(nc{var6}(:,ind,:,:));
      v7 = squeeze(nc{var7}(:,:,:));
    nc = close(nc);

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

    e1(tind,:,:,:) = v1-v2.*v3;
    e2(tind,:,:,:) = v4+v5-v6.^2-v3.^2;
    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];
    vtw(i,:,:,:) = mean(e1(indw,:,:,:));
    vtc(i,:,:,:) = mean(e1(indc,:,:,:));
    kew(i,:,:,:) = mean(e2(indw,:,:,:));
    kec(i,:,:,:) = mean(e2(indc,:,:,:));
    pw(i,:,:) = mean(ps(indw,:,:));
    pc(i,:,:) = mean(ps(indc,:,:));
  end

  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];
  vtw = atlev(vtw, newlev, pw, hyam, hybm, P0);
  vtc = atlev(vtc, newlev, pc, hyam, hybm, P0);
  kew = atlev(kew, newlev, pw, hyam, hybm, P0);
  kec = atlev(kec, newlev, pc, hyam, hybm, P0);
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
  %save eddy_interped_ctdom.mat vtw vtc kew kec newlev lat lon

else
  cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
  load eddy_interped_ctsom.mat
%  load eddy_interped_ctdom.mat
end
get_global; default_global;

levt = 5; levb = 3;
[tstatt, dof] = ttest(vtw(:,levt,:,:), vtc(:,levt,:,:));
scoret = tscore(dof, 2.5);
top = squeeze(mean(vtw(:,levt,:,:) - vtc(:,levt,:,:)))/2;
cintt = 1; clevt = [-5:cintt:5];
[tstatb, dof] = ttest(kew(:,levb,:,:), kec(:,levb,:,:));
scoreb = tscore(dof, 2.5);
bot = squeeze(mean(kew(:,levb,:,:) - kec(:,levb,:,:)))/2;
cintb = 15; clevb = [-90:cintb:90];
tit1 = ['NDJFMAM CT Run:  ' num2str(newlev(levt)) 'mb v''T'', (WCT - CCT)/2'];
tit2 = ['NDJFMAM CT Run:  ' num2str(newlev(levb)) 'mb EKE, (WCT - CCT)/2'];
clabt = [' K m s^-^1']; clabb = [' m^2 s^-^2'];

lev = 5;
top = squeeze(mean(vtw(:,lev,:,:)));
bot = squeeze(mean(vtc(:,lev,:,:)));
cintt = 7.5; clevt = [-60:cintt:60]; cintb = cintt; clevb = clevt;
clabt = [' K m s^-^1']; clabb = clabt;
tit1 = ['WCT ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM'];
tit2 = ['CCT ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM'];

lev = 3;
top = squeeze(mean(kew(:,lev,:,:)));
bot = squeeze(mean(kec(:,lev,:,:)));
cintt = 100; clevt = [-600:cintt:600]; cintb = cintt; clevb = clevt;
clabt = [' K m s^-^1']; clabb = clabt;
tit1 = ['WCT ' num2str(newlev(lev)) 'mb EKE, NDJFMAM'];
tit2 = ['CCT ' num2str(newlev(lev)) 'mb EKE, NDJFMAM'];

figure(1); FRAME = [0 360 -90 90];

sp(1)
     if ismap(gca); clma; end;
     mcont(top, clevt);
     title(tit1); 
     xlabel(['Contour Interval:  ' num2str(cintt) clabt])
%     greyshd(abs(tstatt), scoret);
sp(2)
     if ismap(gca); clma; end;
     mcont(bot, clevb);
     title(tit2); 
     xlabel(['Contour Interval:  ' num2str(cintb) clabb])
%     greyshd(abs(tstatb), scoreb);