Documentation of zonal_eddy

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

Cross-Reference Information

Listing of script zonal_eddy

cd /home/disk/tao/dvimont/matlab/CCM/CT/Vars
load 500mb_vars.mat
load 250mb_vars.mat
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data
filin = 'wct.nc';
nc = netcdf(filin, 'nowrite');
  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/dvimont/matlab/CCM/CT
vorw500 = sph_curl(squeeze(mean(uw500)), squeeze(mean(vw500)), lon, lat);
vorc500 = sph_curl(squeeze(mean(uc500)), squeeze(mean(vc500)), lon, lat);
vorw250 = sph_curl(squeeze(mean(uw250)), squeeze(mean(vw250)), lon, lat);
vorc250 = sph_curl(squeeze(mean(uc250)), squeeze(mean(vc250)), lon, lat);
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(vw250);
global XAX YAX FRAME
XAX = lon;
YAX = lat;
FRAME = [0 360 -90 90];
zon = ones(nlat, 1) * mean(vorw250);
top = vorw250 - zon;
zon = ones(nlat, 1) * mean(vorc250);
bot = vorc250 - zon;
top = vorw250;
bot = vorc250;
range1 = [-30:1:30]*1e-5;
range2 = [-10:.5:10]*1e-5;
top = divw250;
bot = divc250;
range1 = [-30:1:30]*1e-6;
range2 = [-10:1:10]*1e-6;
figure(1)
sp(1)
     mcont(top, range1);
sp(2)
     mcont(top, range1, [-90 270]);
figure(2)
sp(1)
     mcont(bot, range1);
sp(2)
     mcont(bot, range1, [-90 270]);
figure(3)
sp(1)
     mcont(top-bot, range2, 'miller');
sp(2)
     mcont(top-bot, range2, [-90 270]);