Documentation of hgt_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 hgt_gr_compare

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/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{'Z3'}(1:3,:,:,:);
  psw = nc{'PS'}(1:3,:,:);
  tsw = nc{'TS'}(1:3,:,:);
nc = close(nc);

filin = 'cgr_TropPac.nc';
nc = netcdf(filin, 'nowrite');
  hgtc = nc{'Z3'}(1:3,:,:,:);
  psc = nc{'PS'}(1:3,:,:);
  tsc = nc{'TS'}(1:3,:,:);
  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 = 'hgt.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/dvimont/matlab/CCM/GR

filin = 'hgt.mon.mean.nc';
var_nam = 'hgt';
lev = 1000;
[hgt1000gr, hgt1000cl] = reg_gr(filin, var_nam, lev);
lev = 500;
[hgt500gr, hgt500cl] = reg_gr(filin, var_nam, lev);
lev = 250;
[hgt250gr, hgt250cl] = 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
%save hgtlp.mat hgt1000gr hgt1000cl hgt500gr hgt500cl hgt250gr hgt250cl hgtw_nmc hgtc_nmc
%load hgt.mat
load hgtlp.mat

cd /home/disk/tao/dvimont/matlab/CCM/GR
[lat, lon] = getll('wgr.nc');
[nmclat, nmclon] = getll('hgt.mon.mean.nc');
[nmclev] = getnc('hgt.mon.mean.nc', [1:17], 'level');

define_global

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

top = hgt250gr;
lev = 250;
cint = 5;

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

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


XAX = lon;
YAX = lat;

lev = 250;
cint = 5;

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

figure(2)
sp(1)
mcont(top, [-500:cint:500], [90 270]);
title(['GR Run:  NH ' num2str(lev) 'mb HGT (WARM - COLD)/2']);
xlabel(['Contour Interval:  ' num2str(cint) 'm'])
sp(2)
mcont(1 * top, [-500:cint:500], [-90 270]);
title(['GR Run:  SH ' num2str(lev) 'mb HGT (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 individual maps of hgtw_nmc and hgtc_nmc

lev = 500;
cint = 50;

lind = find(nmclev == lev);
top = squeeze(mean(hgtw_nmc(:,lind,:,:)));
top = top - (mean2(top'))' * ones(1, size(top, 2));
bot = squeeze(mean(hgtc_nmc(:,lind,:,:)));
bot = bot - (mean2(bot'))' * ones(1, size(bot, 2));

figure(1)
sp(1)
mcont(top, [-500:cint:500], [90 270]);
title(['GR Run:  NH ' num2str(lev) 'mb HGT (WARM - ZONAL\_MEAN)']);
xlabel(['Contour Interval:  ' num2str(2*cint) 'm'])
sp(2)
mcont(bot, [-500:cint:500], [90 270]);
title(['GR Run:  NH ' num2str(lev) 'mb HGT (COLD - ZONAL\_MEAN)']);
xlabel(['Contour Interval:  ' num2str(2*cint) 'm'])

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