Global Index (short | long) | Local contents | Local Index (short | long)
Get ubar_dHCprimedx and vbar_dHCprimedy
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat
cd ~/matlab/CSIRO/Heat/Old_routines
nfrm = -8:5; lags = nfrm;
tim = 101:1000;
lev = 1:3;
lims = [106 312 1.5 72];
tit = '0-270m LP10';
[ubtp, vbtp, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
[uptb, vptb, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
% Plot the data
[lat_out, lon_out, depth_out] = getll('u', lims);
default_global; XAX = lon_out; YAX = lat_out; FRAME = [110 295 5 65];
lag = 0;
tem1 = -1*ubtp;
tem1 = -1*uptb;
figure(1); figure_orient(1);
cint = 1; clev = [-10:cint:-cint cint:cint:10];
for i = 1:7;
tem = -1 * squeeze(tem1(i,:,:));
subplot(7,2,2*i-1);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
tem = -1 * squeeze(tem1(i+7,:,:));
subplot(7,2,2*i);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i+7))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
for i = 1:2;
subplot(7,2,i)
title([tit ' Ubar dT''dx']);
title([tit ' U'' dTbardx']);
end
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap4/Plots
%print -dps2 ubtp_terms_0-270m_PC1.ps
%print -dps2 uptb_terms_0-270m_PC1.ps
% Plot the vdtdy data
default_global; XAX = lon_out; YAX = lat_out; FRAME = [110 295 5 65];
lag = 0;
figure(2); figure_orient(1);
tem1 = -1*vbtp;
tem1 = -1*vptb;
cint = 1; clev = [-10:cint:-cint cint:cint:10];
for i = 1:7;
tem = -1 * squeeze(tem1(i,:,:));
subplot(7,2,2*i-1);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
tem = -1 * squeeze(tem1(i+7,:,:));
subplot(7,2,2*i);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i+7))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
for i = 1:2;
subplot(7,2,i)
title([tit ' Vbar dT''dy']);
title([tit ' V'' dTbardy']);
end
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap4/Plots
%print -dps2 vbtp_terms_0-270m_PC1.ps
%print -dps2 vptb_terms_0-270m_PC1.ps
% Get the wdtdz terms
cd ~/matlab/CSIRO/Heat/Old_routines
[wbtp, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
[wptb, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
% Plot the wdtdz data
[lat_out, lon_out, depth_out] = getll('wl', lims);
default_global; XAX = lon_out; YAX = lat_out; FRAME = [110 295 5 65];
lag = 0;
figure(3); figure_orient(1);
tem1 = -1*wbtp;
tem1 = -1*wptb;
cint = 1; clev = [-10:cint:-cint cint:cint:10];
for i = 1:7;
tem = -1 * squeeze(tem1(i,:,:));
subplot(7,2,2*i-1);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
tem = -1 * squeeze(tem1(i+7,:,:));
subplot(7,2,2*i);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i+7))]);
if i == 7;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
for i = 1:2;
subplot(7,2,i)
title([tit ' Wbar dT''dz']);
title([tit ' W'' dTbardz']);
end
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap4/Plots
%print -dps2 wbtp_terms_0-270m_PC1.ps
%print -dps2 wptb_terms_0-270m_PC1.ps
% Plot the heat flux and actual heat content terms
hflx = getnc('heat', lims, lev, tim);
[lat, lon, depth] = getll('heat', lims);
cd ~/matlab/CSIRO/Heat
[dhcdt, lat_out, lon_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);
hflx_reg = regress_ceof(hflx, pcs, nfrm);
% Plot the HC and HFLX data
default_global; XAX = lon; YAX = lat; FRAME = [110 295 5 65];
lag = 0;
cint = 1; clev = [-10:cint:-cint cint:cint:10];
figure(4); figure_orient;
for i = 1:6;
tem = squeeze(hflx_reg(i,:,:));
subplot(6,2,2*i);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i))]);
if i == 6;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
subplot(6,2,2)
title([tit ' Heat Flux (Pos Down)']);
for i = 1:6;
tem = squeeze(dhcdt(i,:,:));
subplot(6,2,2*i-1);
gcont(tem, clev);
dc2(tem);
ylabel(['Lag = ' num2str(lags(i))]);
if i == 6;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
subplot(6,2,1)
title([tit ' 0-270m HC Tendency']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 HFLX_dHCdt_0-270m_CPC1.ps
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
%save LP10_hctend_L1-7_yr101-1000_ML.mat uptb ubtp vptb vbtp wptb wbtp lags tim lev nfrm
%save LP10_hctend_L1-3_yr101-1000_ML.mat uptb ubtp vptb vbtp wptb wbtp lags tim lev nfrm
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_hctend_L1-7_yr101-1000_ML.mat
uptb7 = uptb;
ubtp7 = ubtp;
vptb7 = vptb;
vbtp7 = vbtp;
wptb7 = wptb;
wbtp7 = wbtp;
load LP10_hctend_L1-3_yr101-1000_ML.mat
uptb3 = uptb;
ubtp3 = ubtp;
vptb3 = vptb;
vbtp3 = vbtp;
wptb3 = wptb;
wbtp3 = wbtp;
uptb = uptb7 - uptb3;
ubtp = ubtp7 - ubtp3;
vptb = vptb7 - vptb3;
vbtp = vbtp7 - vbtp3;
wptb = wptb7 - wptb3;
wbtp = wbtp7 - wbtp3;
dhcdt7 = dhcdt;
dhcdt3 = dhcdt;
dhcdt = dhcdt7 - dhcdt3;