Documentation of gr_ice_noice

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

Help text

  Plot out patterns regressed on LP filtered GR time series:

Cross-Reference Information

Listing of script gr_ice_noice

clear
if 0;
ind = 11;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
var1 = getnc('Z3_wgr.nc', ind, 'Z3');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
var2 = getnc('Z3_cgr.nc', ind, 'Z3');
[ntim, nlev, nlat, nlon] = size(var1);
[lat, lon] = getll('Z3_cgr.nc');
end
if 0;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data
[a, b] = getnc('PREC_wgr.nc', 'PRECL', 'PRECC');
var1 = (a+b)*(1200*24*3600);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data
[a, b] = getnc('PREC_cgr.nc', 'PRECL', 'PRECC');
var2 = (a+b)*(1200*24*3600);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data
[a, b] = getnc('PREC_clim.nc', 'PRECL', 'PRECC');
var3 = (a+b)*(1200*24*3600);
[ntim, nlat, nlon] = size(var1);
[lat, lon] = getll('Z3_cgr.nc');
end
get_global
default_global
win = sort([2:12:ntim 3:12:ntim 4:12:ntim]);
cint = 7.5;
var1win = squeeze(var1(win,:,:,:));
var2win = squeeze(var2(win,:,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/2;
figure(2)
sp(1);
  clma
  mcont(top, [-90:cint:110], [90 270]);
  greyshd(abs(tstat), score);
  title('500mb Z3:  (WGR\_NOICE - CGR\_NOICE) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' m']);
sp(2);
  clma
  mcont(top, [-90:cint:110], [-90 270]);
  greyshd(abs(tstat), score);
  title('500mb Z3:  (WGR\_NOICE - CGR\_NOICE) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' m']);
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots
figure(2)
FRAME = [40 300 -45 45];
cint = 1;
var1win = squeeze(var1(win,:,:,:));
var2win = squeeze(var3(win,:,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/1;
sp(1);
  clma
  mcont(top, [-10:cint:10], 'giso');
  greyshd(abs(tstat), score);
  title('PRECIP:  (WGR\_NOICE - CGR\_NOICE) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);
var1win = squeeze(var2(win,:,:,:));
var2win = squeeze(var3(win,:,:,:));
[tstat, dof] = ttest(var1win, var2win);
score = tscore(dof, 2.5);
top = squeeze(mean(var1win - var2win))/1;
sp(2);
  clma
  mcont(top, [-10:cint:10], 'giso');
  greyshd(abs(tstat), score);
  title('PRECIP:  (WGR\_NOICE - CGR\_NOICE) / 2');
  xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);

get_global
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'hgt.mon.mean.nc';
[lat, lon] = getll(filin);
default_global

cd /home/disk/tao/dvimont/matlab/CCM/GR
filin = 'hgt.mon.mean.nc';
var = 'hgt';
[pat500, clim] = reg_gr(filin, var, 250);

figure(1)
sp(1)
     cla
     mcont(pat500, [-100:10:100]);
     title('NMC 250mb HEIGHT REGRESSED ON GR')
     xlabel('Contour Interval:  10 m (std)^-^1')
sp(2)
     cla
     mcont(pat500, [-100:10:100], [-90 270]);
     title('NMC 250mb HEIGHT REGRESSED ON GR')
     xlabel('Contour Interval:  10 m (std)^-^1')

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