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