Documentation of airt_gr_compare

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

Help text

  Load GR pattern and time series

Cross-Reference Information

Listing of script airt_gr_compare

cd /home/disk/hayes2/dvimont/data
load ct_gr.mat;

%  Get CCM data:

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
filin = 'wgr_TropPac.nc';
nc = netcdf(filin, 'nowrite');
  hgtw = nc{'AIRT'}(1:3,:,:,:);
  psw = nc{'PS'}(1:3,:,:);
  prw = (nc{'PRECL'}(:) + nc{'PRECC'}(:));
  tsw = nc{'TS'}(1:3,:,:);
nc = close(nc);

filin = 'cgr_TropPac.nc';
nc = netcdf(filin, 'nowrite');
  hgtc = nc{'AIRT'}(1:3,:,:,:);
  psc = nc{'PS'}(1:3,:,:);
  tsc = nc{'TS'}(1:3,:,:);
  prc = (nc{'PRECL'}(:) + nc{'PRECC'}(:));
  lat  = nc{'lat'}(:);
  lon  = nc{'lon'}(:);
  hybm = nc{'hybm'}(:);
  hyam = nc{'hyam'}(:);
  levs = nc{'lev'}(:);
  P0   = nc{'P0'}(:);
nc = close(nc);

%  Get NMC data:

cd /home/disk/tao/data/nmc.reanalysis/monthly

filin = 'prate.mon.mean.nc';
%filin = 'air.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  nmclev = nc{'level'}(:);
  nmclat = nc{'lat'}(:);
  nmclon = nc{'lon'}(:);
nc = close(nc);

hgtw_nmc = atlev(hgtw, nmclev, psw, hyam, hybm, P0);
hgtc_nmc = atlev(hgtc, nmclev, psc, hyam, hybm, P0);

cd /home/disk/tao/data/nmc.reanalysis/monthly
cd /home/disk/tao/dvimont/matlab/CCM/GR

filin = 'air.mon.mean.srfc.nc';
var_nam = 'air';
[sfctgr, sfctcl] = reg_gr(filin, var_nam, 1);

filin = 'prate.mon.mean.nc';
var_nam = 'prate';
[prategr, pratecl] = reg_gr(filin, var_nam, 1);
prategrt = prategr * (3600 * 24 * 1e3 / 1e3);

filin = 'air.mon.mean.nc';
var_nam = 'air';
lev = 1000;
[hgt1000gr, hgt1000cl] = reg_gr(filin, var_nam, lev);
lev = 500;
[hgt500gr, hgt500cl] = reg_gr(filin, var_nam, lev);
lev = 200;
[hgt200gr, hgt200cl] = reg_gr(filin, var_nam, lev);
lev = 850;
[hgt850gr, hgt850cl] = reg_gr(filin, var_nam, lev);


%cd /home/disk/tao/dvimont/matlab/CCM/GR/DATA
%load hgt.mat

define_global

XAX = nmclon;
YAX = nmclat;
FRAME = [0 360 -90 90];

top = hgt200gr;
lev = 200;
cint = 0.2;

figure(1)
sp(1)
mcont(-1 * top, [-10:cint:10], [90 270]);
title(['NMC:  NH ' num2str(lev) 'mb TEMP Regressed on GR']);
xlabel(['Contour Interval:  ' num2str(cint) ' K (std)^-^1'])
sp(2)
mcont(-1 * top, [-10:cint:10], [-90 270]);
title(['NMC:  SH ' num2str(lev) 'mb TEMP Regressed on GR']);
xlabel(['Contour Interval:  ' num2str(cint) ' K (std)^-^1'])

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


XAX = lon;
YAX = lat;

lev = 200;
cint = 5;

lind = find(nmclev == lev);
top = squeeze(mean(hgtw_nmc(:,lind,:,:) - hgtc_nmc(:,lind,:,:)))/2;

figure(1)
sp(1)
mcont(1 * top, [-50:cint:50], [90 270]);
title(['GR Run:  NH ' num2str(lev) 'mb TEMP (WARM - COLD)/2']);
xlabel(['Contour Interval:  ' num2str(cint) 'm'])
sp(2)
mcont(1 * top, [-50:cint:50], [-90 270]);
title(['GR Run:  SH ' num2str(lev) 'mb TEMP (WARM - COLD)/2']);
xlabel(['Contour Interval:  ' num2str(cint) 'm'])

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

%  Check TS out:

figure(2)
orient landscape
gcont(squeeze(mean(tsw - tsc)/2), [-2:.25:2]);
dc

%  Look at global maps of airt

top = sfctgr;
bot = squeeze(mean(tsw(:,:,:) - tsc(:,:,:)))/2;
cint = 0.25

top = prategrt;
bot = squeeze(mean(prw - prc)) * (3600 * 24 * 1e3 / 2);
cint = 0.5;
yticks = 0; xticks = 0;

figure(1)

FRAME = [0 360 -90 90];
sp(1)
XAX = nmclon;
YAX = nmclat;
mcont(1*top, [-5:cint:5], 'gortho', [0 180]);
manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
title(['NMC:  PRECIP Regressed on GR']);
xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);
sp(2)
XAX = lon;
YAX = lat;
mcont(bot, [-5:cint:5], 'gortho', [0 180]);
manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
title(['PRECIP (WARM - COLD)/2:  GR\_TPAC run']);
xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);

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