Documentation of zonal_u

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

Help text

  Just a side note here.  Recall that the SSTA pattern
  used to force this run is weighted by some variable
  moave which is the mean of the warm event and cold event
  time progressions from SEP to APR.  This is contained
  in ct_sst.nc.  The average value for moave (in units of
  standard deviations of the CT* pattern) is 1.446.  So,
  I think we want to divide the linear difference between
  warm-cold by 2*1.446 in order to get reasonable values
  for the 500mb response.

Cross-Reference Information

Listing of script zonal_u

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data

filin = 'wct.nc';
nc = netcdf(filin, 'nowrite');
  pslw = nc{'PSL'}(:);
  hgtw = nc{'Z3'}(:);
  airtw = nc{'AIRT'}(:);
  psw = nc{'PS'}(:);
  latct = nc{'lat'}(:);
  lonct = nc{'lon'}(:);
  hyam = nc{'hyam'}(:);
  hybm = nc{'hybm'}(:);
  P0 = nc{'P0'}(:);
nc = close(nc);

filin = 'cct.nc';
nc = netcdf(filin, 'nowrite');
  pslc = nc{'PSL'}(:);
  hgtc = nc{'Z3'}(:);
  airtc = nc{'AIRT'}(:);
  psc = nc{'PS'}(:);
nc = close(nc);

pslw = squeeze(mean(pslw));
psw = squeeze(mean(psw));
hgtw = squeeze(mean(hgtw));
airtw = squeeze(mean(airtw));
pslc = squeeze(mean(pslc));
psc = squeeze(mean(psc));
hgtc = squeeze(mean(hgtc));
airtc = squeeze(mean(airtc));

weight = 2*1.446;

global XAX YAX FRAME
XAX = lonct;
YAX = latct;
FRAME = [0 360 -90 90];

%  Compare with NMC data