Documentation of prec_compare

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

Cross-Reference Information

Listing of script prec_compare

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data
filin = 'wct.nc';
nc = netcdf(filin, 'nowrite');
  pw = (nc{'PRECL'}(:) + nc{'PRECC'}) * 24 * 3600 * 1000;
  lat  = nc{'lat'}(:);
  lon  = nc{'lon'}(:);
nc = close(nc);
filin = 'cct.nc';
nc = netcdf(filin, 'nowrite');
  pc = (nc{'PRECL'}(:) + nc{'PRECC'}) * 24 * 3600 * 1000;
nc = close(nc);
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'prate.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  nmclev = nc{'level'}(:);
  nmclat = nc{'lat'}(:);
  nmclon = nc{'lon'}(:);
nc = close(nc);
cd /home/disk/tao/dvimont/matlab/CCM/CT
filin = 'prate.mon.mean.nc';
var_nam = 'prate';
[prategr, pratecl] = reg_ct(filin, var_nam, 1);
prategrt = prategr * (3600 * 24 * 1e3 / 1e3);
get_global
define_global
top = prategrt;
bot = squeeze(mean(pw - pc))/2;
cint = 1;
yticks = 0; xticks = 0;
figure(1)
sp(1)
XAX = nmclon;
YAX = nmclat;
mcont(top, [-10:cint:10], 'gortho', [0 180]);
manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
title(['NMC:  PRECIP Regressed on CT']);
xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);
sp(2)
XAX = lon;
YAX = lat;
mcont(bot, [-10:cint:10], 'gortho', [0 180]);
manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
title(['CT Run:  PRECIP (WARM - COLD)/2']);
xlabel(['Contour Interval:  ' num2str(cint) ' mm (day)^-^1']);
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots