Global Index (short | long) | Local contents | Local Index (short | long)
Load PCS:
| This script calls | |
|---|---|
clear
biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
load LP10_L1-7_CEOF.mat; tit = 'LP10';
elseif biff2 == 2;
load HP8_L1-7_CEOF.mat; tit = 'HP8';
elseif biff2 == 3;
load RAW_L1-7_CEOF.mat; tit = 'RAW';
end
cd ~/matlab/CSIRO/Heat/Temp_tend
nfrm = 6;
tim = 101:550;
lev = 5;
lims = [106 312 -63 63];
% Get ubar_dHCprimedx and vbar_dHCprimedy
[ubtp, vbtp, lat_out, lon_out, depth_out] = ...
get_ubtp(pcs, lims, nfrm, tim, lev);
[uptb, vptb, lat_out, lon_out, depth_out] = ...
get_uptb(pcs, lims, nfrm, tim, lev);
% Plot the data
default_global; XAX = lon_out; YAX = lat_out; FRAME = [110 299 -60 60];
lag = 0;
figure(1); figure_orient;
cint = 1; clev = [-20:cint:-cint cint:cint:20];
for i = 1:6;
tem = -1e9 * squeeze(ubtp(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm Ubar dT''/dx regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_ubtp_L5.ps
figure(1); figure_orient;
for i = 1:6;
tem = -1e9 * squeeze(uptb(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm U'' dTbar/dx regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_uptb_L5.ps
% Plot the vdtdy data
figure(1); figure_orient;
for i = 1:6;
tem = -1e9 * squeeze(vbtp(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm Vbar dT''/dy regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_vbtp_L5.ps
figure(1); figure_orient;
for i = 1:6;
tem = -1e9 * squeeze(vptb(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm V'' dTbar/dy regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_vptb_L5.ps
% Get the wdtdz terms
lev = lev;
cd ~/matlab/CSIRO/Heat/Temp_tend
[wbtp, lat_out, lon_out, depth_out] = ...
get_wbtp(pcs, lims, nfrm, tim, lev);
[wptb, lat_out, lon_out] = ...
get_wptb(pcs, lims, nfrm, tim, lev);
% Plot the wdtdz data
default_global; XAX = lon_out; YAX = lat_out; FRAME = [110 299 -60 60];
lag = 0;
figure(1); figure_orient;
for i = 1:6;
tem = -1e9 * squeeze(wbtp(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm Wbar dT''/dz regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_wbtp_L5.ps
figure(1); figure_orient;
for i = 1:6;
tem = -1e9 * squeeze(wptb(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm W'' dTbar/dz regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_wptb_L5.ps
% Plot the heat flux and actual heat content terms
hflx = getnc('heat', lims, lev, tim);
[lat, lon, depth] = getll('heat', lims);
heat = getnc('temp', lims, lev, tim);
hflx_reg = regress_ceof(hflx, pcs, nfrm);
heat_reg = regress_ceof(heat, pcs, nfrm);
% Plot the HC and HFLX data
default_global; XAX = lon; YAX = lat; FRAME = [110 299 -60 60];
lag = 0;
figure(1); figure_orient;
cint = 1; clev = [-20:cint:-cint cint:cint:20];
for i = 1:6;
tem = squeeze(hflx_reg(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2']);
end
end
subplot(3,2,3)
ylabel([tit ': Downward Heat Flux regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_HFLX_down.ps
figure(1); figure_orient;
cint = 0.05; clev = [-10:cint:-cint 0 cint:cint:10];
for i = 1:6;
tem = squeeze(heat_reg(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' K std^-^1']);
end
end
subplot(3,2,3)
ylabel([tit ': ' num2str(depth_out) ...
'm Temperature regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_Temp_L5.ps
% Plot sfc temp and zonal wind
taux = getnc('taux', lims, lev, tim);
[latx, lonx, depthx] = getll('taux', lims);
xreg = regress_ceof(taux, pcs, nfrm);
% Plot the HC and HFLX data
default_global; XAX = lonx; YAX = latx; FRAME = [110 299 -60 60];
lag = 0;
figure(1); figure_orient;
cint = 1; clev = [-10:cint:-cint 0 cint:cint:10];
for i = 1:6;
tem = 100*squeeze(xreg(i,:,:));
subplot(3,2,i);
gcont(tem, clev);
dc2(tem);
title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
if i >= 5;
xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^2 dynes cm^-^2']);
end
end
subplot(3,2,3)
ylabel([tit ': Zonal Wind Stress regressed on CPC1 of 270m Heat Content']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Temp_tend/Plot_adv/LP10
print -dps2 LP10_Taux.ps