Documentation of plot_vort

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

Help text

  Try some plots removing the zonal average.

Cross-Reference Information

Listing of script plot_vort

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data
filin = 'wct.nc';
nc = netcdf(filin, 'nowrite');
  vtw  = nc{'VT'}(1:3,:,:,:);
  kew  = nc{'VVPUU'}(1:3,:,:,:);
  uw   = nc{'U'}(1:3,:,:,:);
  vw   = nc{'V'}(1:3,:,:,:);
  airw = nc{'AIRT'}(1:3,:,:,:);
  psw  = nc{'PS'}(1:3,:,:);
  omw  = nc{'OMEGA'}(1:3,:,:);
  zw = nc{'Z3'}(1:3,:,:,:);
nc = close(nc);
filin = 'cct.nc';
nc = netcdf(filin, 'nowrite');
  vtc  = nc{'VT'}(1:3,:,:,:);
  kec  = nc{'VVPUU'}(1:3,:,:,:);
  uc   = nc{'U'}(1:3,:,:,:);
  vc   = nc{'V'}(1:3,:,:,:);
  airc = nc{'AIRT'}(1:3,:,:,:);
  psc  = nc{'PS'}(1:3,:,:);
  omc  = nc{'OMEGA'}(1:3,:,:);
  zc = nc{'Z3'}(1:3,:,:,:);
  lat  = nc{'lat'}(:);
  lon  = nc{'lon'}(:);
  hybm = nc{'hybm'}(:);
  hyam = nc{'hyam'}(:);
  levs = nc{'lev'}(:);
  P0   = nc{'P0'}(:);
nc = close(nc);
lev = 500;
uw500 = atlev(uw, lev, psw, hyam, hybm, P0);
vw500 = atlev(vw, lev, psw, hyam, hybm, P0);
vorw500 = sph_curl(squeeze(mean(uw500)), squeeze(mean(vw500)), lon, lat);
uc500 = atlev(uc, lev, psc, hyam, hybm, P0);
vc500 = atlev(vc, lev, psc, hyam, hybm, P0);
vorc500 = sph_curl(squeeze(mean(uc500)), squeeze(mean(vc500)), lon, lat);
zw500 = atlev(zw, lev, psw, hyam, hybm, P0);
zc500 = atlev(zc, lev, psc, hyam, hybm, P0);
cd /home/disk/tao/dvimont/matlab/CCM/CT/Vars
save 500mb_vars.mat uw500 vw500 zw500 uc500 vc500 zc500
lev = 250;
uw250 = atlev(uw, lev, psc, hyam, hybm, P0);
vw250 = atlev(vw, lev, psc, hyam, hybm, P0);
vorw250 = sph_curl(squeeze(mean(uw250)), squeeze(mean(vw250)), lon, lat);
uc250 = atlev(uc, lev, psc, hyam, hybm, P0);
vc250 = atlev(vc, lev, psc, hyam, hybm, P0);
vorc250 = sph_curl(squeeze(mean(uc250)), squeeze(mean(vc250)), lon, lat);
zw250 = atlev(zw, lev, psw, hyam, hybm, P0);
zc250 = atlev(zc, lev, psc, hyam, hybm, P0);
save 250mb_vars.mat uw250 vw250 zw250 uc250 vc250 zc250
global XAX YAX FRAME
XAX = lon;
YAX = lat;
FRAME = [0 360 -90 90];
top = vorw250;
bot = vorc250;
top = vorw500;
bot = vorc500;
top = divw250;
bot = divc250;
top = divw500;
bot = divc500;
figure(1)
range1 = [-30:1:30]*1e-5;
sp(1)
     mcont(top, range1);
sp(2)
     mcont(bot, range1);
figure(2)
     range2 = [-10:.5:10]*1e-5;
     mcont(top-bot, range2);
divw500 = sph_div(squeeze(mean(uw500)),squeeze(mean(vw500)),lon,lat);
divc500 = sph_div(squeeze(mean(uc500)),squeeze(mean(vc500)),lon,lat);
divw250 = sph_div(squeeze(mean(uw250)),squeeze(mean(vw250)),lon,lat);
divc250 = sph_div(squeeze(mean(uc250)),squeeze(mean(vc250)),lon,lat);

[ntim, nlat, nlon] = size(zw500);

zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zw500, 2))), 1, ntim*nlat);
zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon);
zw5star = zw500 - zon;
zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zc500, 2))), 1, ntim*nlat);
zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon);
zc5star = zc500 - zon;

zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zw250, 2))), 1, ntim*nlat);
zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon);
zw25star = zw250 - zon;
zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zc250, 2))), 1, ntim*nlat);
zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon);
zc25star = zc250 - zon;

top = squeeze(mean(zw5star));
bot = squeeze(mean(zc5star));

top = squeeze(mean(zw25star));
bot = squeeze(mean(zc25star));

figure(1)
     range = [-500:50:500];
sp(1)
     mcont(top, range)
sp(2)
     mcont(bot, range)

figure(2)
     range = [-400:20:400];
     mcont(top-bot, range)