Documentation of airt_compare

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

Cross-Reference Information

Listing of script airt_compare

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data
filin = 'wct.nc';
nc = netcdf(filin, 'nowrite');
  airtw = nc{'AIRT'}(1:3,:,:,:);
  tsw = nc{'TS'}(:);
  ccmpslw = nc{'PSL'}(:);
  psw = nc{'PS'}(:);
nc = close(nc);
filin = 'cct.nc';
nc = netcdf(filin, 'nowrite');
  airtc = nc{'AIRT'}(1:3,:,:,:);
  tsc = nc{'TS'}(:);
  ccmpslc = nc{'PSL'}(:);
  psc = nc{'PS'}(:);
  lat  = nc{'lat'}(:);
  lon  = nc{'lon'}(:);
  hybm = nc{'hybm'}(:);
  hyam = nc{'hyam'}(:);
  levs = nc{'lev'}(:);
  P0   = nc{'P0'}(:);
nc = close(nc);
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'air.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  nmclev = nc{'level'}(:);
  nmclat = nc{'lat'}(:);
  nmclon = nc{'lon'}(:);
nc = close(nc);
airtw_nmc = atlev(airtw, nmclev, psw, hyam, hybm, P0);
airtc_nmc = atlev(airtc, nmclev, psc, hyam, hybm, P0);
cd /home/disk/tao/dvimont/matlab/CCM/CT
filin = 'air.mon.mean.nc';
var_nam = 'air';
lev = 1000;
[air1000ct, air1000c] = reg_ct(filin, var_nam, lev);
lev = 925;
[air925ct, air925c] = reg_ct(filin, var_nam, lev);
lev = 500;
[air500ct, air500c] = reg_ct(filin, var_nam, lev);
lev = 250;
[air250ct, air250c] = reg_ct(filin, var_nam, lev);
filin = 'slp.mon.mean.nc';
var_nam = 'slp';
[pslct, pslc] = reg_ct(filin, var_nam, lev);
top = squeeze(mean(airtw(:,2,:,:) - airtc(:,2,:,:)))/(2*1.4);
bot = air925ct;
top = squeeze(mean(tsw - tsc))/(2*1.4);
bot = air1000ct;
int = 0.25;
lev = find(nmclev == 250);
top = squeeze(mean(airtw(:,lev,:,:)-airtc(:,lev,:,:)))/(2*1.4);
bot = air250ct;
int = 0.2;
top = squeeze(mean(ccmpslw - ccmpslc))/(2*1.4*100);
bot = pslct;
int = 0.5;
figure(1)
sp(1)
     XAX = lon;
     YAX = lat;
     gcont(top, [-5:int:5]);
     title(['CT run, DJF ' num2str(nmclev(lev)) 'mb Temperature (WARM - COLD) / (2*1.4)']);
     xlabel(['Contour Interval ' num2str(int) ' K'])
     dc
sp(2)
     XAX = nmclon;
     YAX = nmclat;
     gcont(bot, [-5:int:5]);
     title(['NMC DJF ' num2str(nmclev(lev)) 'mb Temperature regressed on CT*']);
     xlabel(['Contour Interval ' num2str(int) ' K / std'])
     dc
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots