Documentation of regress_gr_eddy

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

Help text

  Look at 850mb VT:

Cross-Reference Information

Listing of script regress_gr_eddy

clear
cd /home/disk/hayes2/dvimont/data
load gr_ct_29March1999.mat

cd /home/disk/hayes2/dvimont/nmc/vgrd.tmp
ntim = 480; [nlat, nlon] = size(grpat);
%  Note that 35=November, 460 = April
tkp = 35:460;
gr = gr(tkp);
ct = ctstar(tkp);

lev = 850;
  filin = ['vgrd.tmp.' num2str(lev) '.mon.bin']
  fid = fopen(filin, 'r', 'b');
  tmp = fread(fid, inf, 'float32');
  tmp = reshape(tmp, nlon, nlat, ntim);
  xdim = size(tmp);
  [tmp, xdim] = shift_dim(tmp, xdim, [2 3]);
  tmp = reshape(tmp, ntim, nlat, nlon);
  tmp = reshape(shiftdim(tmp, 1), nlat, nlon*ntim);
  tmp = flipud(tmp);
  tmp = shiftdim(reshape(tmp, nlat, nlon, ntim), 2);
  tmp = tmp(tkp,:,:);
  [tmp, clim] = annave(tmp);

[ntim, nlat, nlon] = size(tmp);
tmp = reshape(tmp, ntim, nlat*nlon);

for i = 1:36;
  ind = 12*(i-1)+[1:6];
  tmtim(i) = mean(gr(ind));
  tempat(i,:) = mean(tmp(ind,:));
end;
[ntim, npt] = size(tempat);
[tempat, tem] = remove_mean(tempat);

for i = 1:nlat*nlon;
  variance(i) = tempat(:,i)'*tempat(:,i)/ntim;
  rcoef(i) = corr(tmtim, tempat(:,i), 0);
end
dofx = get_dof(tmtim);
dofy = get_dof(tempat);
doftot = dofx + dofy - 2;

score = rcoef .* sqrt(doftot) ./ sqrt(1 - rcoef.^2);
tlim = tscore(min(doftot), 2.5);


FRAME = [110 290 15 90];
figure(1);  

tem = tmtim * tempat ./ ntim;
shpat = reshape(score, nlat, nlon);
tem = reshape(tem, nlat, nlon);
cint = 0.75; clev = [-6:cint:6];
sp(1)
     if ismap(gca); clma; end;
     mcont(tem, clev, 'stereo', [90 mean(FRAME(1:2))]);
     greyshd(abs(shpat), tlim);
     title(['NMC:  ' num2str(lev) 'mb v''T'' Regressed on GR, 1960-1995']);
     xlabel(['Contour Interval:  ' num2str(cint) ' K m s^-^1 std^-^1'])

for i = 1:36;
  ind = 12*(i-1)+[1:6];
  tmtim(i) = mean(ct(ind));
  tempat(i,:) = mean(tmp(ind,:));
end;
[ntim, npt] = size(tempat);
[tempat, tem] = remove_mean(tempat);

for i = 1:nlat*nlon;
  variance(i) = tempat(:,i)'*tempat(:,i)/ntim;
  rcoef(i) = corr(tmtim, tempat(:,i), 0);
end
dofx = get_dof(tmtim);
dofy = get_dof(tempat);
doftot = dofx + dofy - 2;

score = rcoef .* sqrt(doftot) ./ sqrt(1 - rcoef.^2);
tlim = tscore(min(doftot), 2.5);

tem = tmtim * tempat ./ ntim;
shpat = reshape(score, nlat, nlon);
tem = reshape(tem, nlat, nlon);
cint = 0.75; clev = [-6:cint:6];
sp(2)
     if ismap(gca); clma; end;
     mcont(tem, clev, 'stereo', [90 mean(FRAME(1:2))]);
     greyshd(abs(shpat), tlim);
     title(['NMC:  ' num2str(lev) 'mb v''T'' Regressed on CT*, 1960-1995']);
     xlabel(['Contour Interval:  ' num2str(cint) ' K m s^-^1 std^-^1'])




cd /home/disk/hayes2/dvimont/nmc/ugrd.vgrd

lev = 300; ntim = 480;
  filin = ['ugrd.vgrd.' num2str(lev) '.mon.bin']
  fid = fopen(filin, 'r', 'b');
  tmp = fread(fid, inf, 'float32');
  tmp = reshape(tmp, nlon, nlat, ntim);
  xdim = size(tmp);
  [tmp, xdim] = shift_dim(tmp, xdim, [2 3]);
  tmp = reshape(tmp, ntim, nlat, nlon);
  tmp = reshape(shiftdim(tmp, 1), nlat, nlon*ntim);
  tmp = flipud(tmp);
  tmp = shiftdim(reshape(tmp, nlat, nlon, ntim), 2);
  tmp = tmp(tkp,:,:);

[ntim, nlat, nlon] = size(tmp);
tmp = reshape(tmp, ntim, nlat*nlon);

FRAME = [110 290 15 90];
figure(1);  

tem = gr1' * tmp(win,:) ./ length(win);
tem = reshape(tem, nlat, nlon);
cint = 2.5; clev = [-60:cint:60];
sp(1)
     if ismap(gca); clma; end;
     mcont(tem, clev, 'stereo', [90 mean(FRAME(1:2))]);
     title(['NMC:  ' num2str(lev) 'mb u''v'' Regressed on GR, 1960-1995']);
     xlabel(['Contour Interval:  ' num2str(cint) ' m^2 s^-^2 std^-^1'])

tem = ct1' * tmp(win,:) ./ length(win);
tem = reshape(tem, nlat, nlon);
cint = 2.5; clev = [-60:cint:60];
sp(2)
     if ismap(gca); clma; end;
     mcont(tem, clev, 'stereo', [90 mean(FRAME(1:2))]);
     title(['NMC:  ' num2str(lev) 'mb u''v'' Regressed on CT*, 1960-1995']);
     xlabel(['Contour Interval:  ' num2str(cint) ' K m s^-^1 std^-^1'])