Documentation of eddy_time_plots2

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

Help text

 Read Variables from wgr and cgr runs

Cross-Reference Information

Listing of script eddy_time_plots2

clear

  var1 = 'VT'; var2 = 'AIRT'; var3 = 'V';
  var4 = 'UU'; var5 = 'VV'; var6 = 'U'; var7 = 'PS';

  lind = 5:18;

for numind = 1:2;

  if numind == 1;
    cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr2/data
    lab = 'DOM2x'
  elseif numind == 2;
    cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr_som2/data
    lab = 'SOM2x'
  end

vtw = []; kew = []; psw = [];
for i = 3:13;
  if i < 10;
    filin = ['000' num2str(i) '.nc'];
  else
    filin = ['00' num2str(i) '.nc'];
  end
  disp(filin)

  [v1, v2, v3, v4, v5, v6, v7] = getnc(filin, lind, ...
     var1, var2, var3, var4, var5, var6, var7);

  vtw = [vtw; v1 - v2.*v3];
  kew = [kew; v4 + v5 - v6.^2 - v3.^2];
  psw = [psw; v7];
end

  if numind == 1;
    cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr2/data
  elseif numind == 2;
    cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr_som2/data
  end

vtc = []; kec = []; psc = [];
for i = 3:13;
  if i < 10;
    filin = ['000' num2str(i) '.nc'];
  else
    filin = ['00' num2str(i) '.nc'];
  end
  disp(filin)

  [v1, v2, v3, v4, v5, v6, v7] = getnc(filin, lind, ...
     var1, var2, var3, var4, var5, var6, var7);

  vtc = [vtc; v1 - v2.*v3];
  kec = [kec; v4 + v5 - v6.^2 - v3.^2];
  psc = [psc; v7];
end

%  Compress variables, so only yearly data is kept

for i = 1:10;
  ind = (i-1)*12 + [1:6];
  vtw1(i,:,:,:) = mean(vtw(ind,:,:,:));
  kew1(i,:,:,:) = mean(kew(ind,:,:,:));
  psw1(i,:,:,:) = mean(psw(ind,:,:,:));
  vtc1(i,:,:,:) = mean(vtc(ind,:,:,:));
  kec1(i,:,:,:) = mean(kec(ind,:,:,:));
  psc1(i,:,:,:) = mean(psc(ind,:,:,:));
end

clear vtw kew psw vtc kec psc

%  Get variables for interpolation

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(vtw1, newlev, psw1, hyam, hybm, P0);
vtc = atlev(vtc1, newlev, psc1, hyam, hybm, P0);
kew = atlev(kew1, newlev, psw1, hyam, hybm, P0);
kec = atlev(kec1, newlev, psc1, hyam, hybm, P0);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
eval(['save eddy_interped_gr' lab '.mat vtw vtc kew kec lat lon newlev'])

end



%  Plot the data

clear

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data
load eddy_interped_grSOM2x.mat
vtwsg = vtw; kewsg = kew;
vtcsg = vtc; kecsg = kec;
load eddy_interped_grDOM2x.mat
vtwdg = vtw; kewdg = kew;
vtcdg = vtc; kecdg = kec;
load eddy_interped_ctsom.mat
vtwsc = vtw; kewsc = kew;
vtcsc = vtc; kecsc = kec;
load eddy_interped_ctdom.mat
vtwdc = vtw; kewdc = kew;
vtcdc = vtc; kecdc = kec;
get_global; default_global;
FRAME = [0 360 -90 90]; origin = [90 270];

lab = 'dg'; tit1 = ['GR\_DOM']; i = 1;
lab = 'sg'; tit1 = ['GR\_SOM']; i = 1;
lab = 'dc'; tit1 = ['CT\_DOM']; i = 2;
lab = 'sc'; tit1 = ['CT\_SOM']; i = 2;
figure(1);

%  VT
lev = find(newlev == 850); tit = ['V''T'''];
cint = 1; clev = [-10:cint:10];
eval(['v1 = squeeze(vtw' lab '(:,lev,:,:));']); 
eval(['v2 = squeeze(vtc' lab '(:,lev,:,:));']); 

[tstat, dof] = ttest(v1, v2);
score = tscore(dof, 2.5);
top = squeeze(mean(v1 - v2))/2;
sp(i)
     if ismap(gca); clma; end;
     gcont(top, clev); dc;%, 'stereo', origin);
%     greyshd(abs(tstat), score);
     title(['NDJFMA ' tit1 ': ' num2str(newlev(lev)) 'mb ' ...
		      tit ', (WARM - COLD)/2']);
     xlabel(['Contour Interval:  ' num2str(cint) ' K m s^-^1']);

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM 
%  print -dps2 EDDY_VT_DS__W-C.ps

%  EKE
lab = 'dg'; tit1 = ['GR\_DOM']; i = 1;
lab = 'sg'; tit1 = ['GR\_SOM']; i = 1;
lab = 'dc'; tit1 = ['CT\_DOM']; i = 2;
lab = 'sc'; tit1 = ['CT\_SOM']; i = 2;
figure(1);

lev = find(newlev == 300); tit = ['EKE'];
cint = 15; clev = [-90:cint:90];
eval(['v1 = squeeze(kew' lab '(:,lev,:,:));']); 
eval(['v2 = squeeze(kec' lab '(:,lev,:,:));']); 

[tstat, dof] = ttest(v1, v2);
score = tscore(dof, 2.5);
top = squeeze(mean(v1 - v2))/2;
sp(i)
     if ismap(gca); clma; end;
     gcont(top, clev); dc;%, 'stereo', origin);
%     greyshd(abs(tstat), score);
     title(['NDJFMA ' tit1 ': ' num2str(newlev(lev)) 'mb ' ...
		      tit ', (WARM - COLD)/2']);
     xlabel(['Contour Interval:  ' num2str(cint) ' m^2 s^-^2']);

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM 
%  print -dps2 EDDY_EKE_DS__W-C.ps