Global Index (short | long) | Local contents | Local Index (short | long)
top = squeeze(mean(vtw500-vta500)); bot = squeeze(mean(vtc500-vta500)); top = squeeze(mean(vtw500)); bot = squeeze(mean(vtc500));
| This script calls | |
|---|---|
cd /cdrom
filin = 'warmenso.nc';
nc = netcdf(filin, 'nowrite');
vtw = nc{'VT'}(2:4,:,:,:);
kew = nc{'VVPUU'}(2:4,:,:,:);
uw = nc{'U'}(2:4,:,:,:);
vw = nc{'V'}(2:4,:,:,:);
airw = nc{'AIRT'}(2:4,:,:,:);
psw = nc{'PS'}(2:4,:,:);
omw = nc{'OMEGA'}(2:4,:,:);
nc = close(nc);
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,:,:);
nc = close(nc);
cd /cdrom
filin = 'coldenso.nc';
nc = netcdf(filin, 'nowrite');
vtc = nc{'VT'}(2:4,:,:,:);
kec = nc{'VVPUU'}(2:4,:,:,:);
uc = nc{'U'}(2:4,:,:,:);
vc = nc{'V'}(2:4,:,:,:);
airc = nc{'AIRT'}(2:4,:,:,:);
psc = nc{'PS'}(2:4,:,:);
omc = nc{'OMEGA'}(2:4,:,:);
nc = close(nc);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data
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,:,:);
lat = nc{'lat'}(:);
lon = nc{'lon'}(:);
hybm = nc{'hyam'}(:);
hyam = nc{'hybm'}(:);
levs = nc{'lev'}(:);
P0 = nc{'P0'}(:);
nc = close(nc);
cd /cdrom
filin = 'clim_dec82-jun98.nc';
nc = netcdf(filin, 'nowrite');
vta = nc{'VT'}(:,:,:,:);
kea = nc{'VVPUU'}(:,:,:,:);
ua = nc{'U'}(:,:,:,:);
va = nc{'V'}(:,:,:,:);
aira = nc{'AIRT'}(:,:,:,:);
psa = nc{'PS'}(:,:,:);
oma = nc{'OMEGA'}(:,:,:,:);
lat = nc{'lat'}(:);
lon = nc{'lon'}(:);
hybm = nc{'hyam'}(:);
hyam = nc{'hybm'}(:);
levs = nc{'lev'}(:);
P0 = nc{'P0'}(:);
nc = close(nc);
vta = vta([12 1 2],:,:,:);
kea = kea([12 1 2],:,:,:);
ua = ua([12 1 2],:,:,:);
va = va([12 1 2],:,:,:);
aira = aira([12 1 2],:,:,:);
oma = oma([12 1 2],:,:,:);
psa = psa([12 1 2],:,:);
lev = 500;
tem = vtw - vw.*airw;
vtw500 = atlev(tem, lev, psw, hyam, hybm, P0);
tem = vtc - vc.*airc;
vtc500 = atlev(tem, lev, psc, hyam, hybm, P0);
tem = vta - va.*aira;
vta500 = atlev(tem, lev, psa, hyam, hybm, P0);
tem = kew - uw.^2 - vw.^2;
ekew500 = atlev(tem, lev, psw, hyam, hybm, P0);
tem = kec - uc.^2 - vc.^2;
ekec500 = atlev(tem, lev, psc, hyam, hybm, P0);
tem = kea - ua.^2 - va.^2;
ekea500 = atlev(tem, lev, psa, hyam, hybm, P0);
lims = [0 360 -90 90];
global XAX YAX FRAME
XAX = lon;
YAX = lat;
FRAME = lims;
figure(1)
range1 = [-40:4:40];
sp(1)
mcont(top, range1);
sp(2)
mcont(bot, range1);
%top = squeeze(mean(ekew500-ekea500));
%bot = squeeze(mean(ekec500-ekea500));
top = squeeze(mean(ekew500));
bot = squeeze(mean(ekec500));
figure(2)
range2 = [-200:50:500];
sp(1);
mcont(top, range2);
sp(2)
mcont(bot, range2);
top = squeeze(mean(vta500));
bot = squeeze(mean(ekea500));
figure(3)
sp(1);
mcont(top, range1);
sp(2)
mcont(bot, range2);
top = squeeze(mean(vtw500-vtc500));
bot = squeeze(mean(ekew500-ekec500));
figure(3)
sp(1);
mcont(top, [-20:2:20]);
sp(2)
mcont(bot, [-100:10:100]);
% Completely different here...
top = squeeze(mean(mean(shiftdim(uw,3))));
bot = squeeze(mean(mean(shiftdim(uc,3))));
range1 = [-30:5:50];
range2 = [-10:1:10];
top = squeeze(mean(mean(shiftdim(omw,3))));
bot = squeeze(mean(mean(shiftdim(omc,3))));
range1 = [-.1:.01:.1];
range2 = [-.1:.0025:.1];
top2 = squeeze(mean(mean(shiftdim(vw,3))));
bot2 = squeeze(mean(mean(shiftdim(vc,3))));
figure(1)
clf
subplot(2,1,1)
pncont(lat, -1*levs, top, range1);
grid on
subplot(2,1,2)
pncont(lat, -1*levs, bot, range1);
grid on
figure(2)
clf
pncont(lat, -1*levs, top-bot, range2)
grid on
figure(1)
clf
subplot(2,1,1)
quiver(lat, -1*levs, top2, -200*top, 15);
grid on
subplot(2,1,2)
quiver(lat, -1*levs, bot2, -200*bot, 15);
grid on
figure(2)
clf
quiver(lat, -1*levs, top2-bot2, -200*(top-bot), 15)
grid on