Global Index (short | long) | Local contents | Local Index (short | long)
lims = [89 310 -15 15];
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs2/data
load heat_pcs.mat; pcs = hpcs(:,1);
nfrm = -6:6; lags = nfrm;
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs2/data
load heat_pcs.mat; pcs = lpcs(:,1);
nfrm = -12:12; lags = nfrm;
ct = getct; ct = detrend(ct);
[b, a] = butter(9, 2/18);
[b2, a2] = butter(9, 2/40);
pcs = filtfilt(b, a, ct)-filtfilt(b2, a2, ct);
nfrm = -16:2:16; lags = nfrm;
ct = getct; ct = detrend(ct);
[b, a] = butter(9, 2/20);
[b2, a2] = butter(9, 2/30);
pcs = filtfilt(b, a, ct)-filtfilt(b2, a2, ct);
nfrm = -16:2:16; lags = nfrm;
ct = getct; ct = detrend(ct);
[b, a] = butter(9, 2/4);
[b2, a2] = butter(9, 2/12);
pcs = filtfilt(b, a, ct)-filtfilt(b2, a2, ct);
nfrm = -6:6; lags = nfrm;
[b2, a2] = butter(9, 2/12);
pcs = ct-filtfilt(b2, a2, ct);
nfrm = -8:8; lags = nfrm;
ct = getct; ct = detrend(ct);
[b2, a2] = butter(9, 2/6);
pcs = ct-filtfilt(b2, a2, ct);
nfrm = -4:4; lags = nfrm;
ct = getct; ct = detrend(ct);
[b2, a2] = butter(9, 2/7);
pcs = filtfilt(b2, a2, ct);
nfrm = -8:8; lags = nfrm;
cd ~/matlab/CSIRO/Heat/Old_routines
tim = 101:1000;
lev = 1:3;
lims = [89 310 -82 82];
% Get ubar_dHCprimedx and vbar_dHCprimedy
[ubtp13, vbtp13, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
[uptb13, vptb13, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
% Get the wdtdz terms
cd ~/matlab/CSIRO/Heat/Old_routines
[wbtp13, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
[wptb13, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
cd ~/matlab/CSIRO/Heat
[dhcdt13, lat_out, lon_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);
hflx = getnc('heat', lims, 1, tim);
[lat, lon, depth] = getll('heat', lims);
hflx_reg = regress_eof(hflx, pcs, nfrm);
% Get ubar_dHCprimedx and vbar_dHCprimedy
cd ~/matlab/CSIRO/Heat/Old_routines
lev = 4:7;
[ubtp47, vbtp47, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
[uptb47, vptb47, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
% Get the wdtdz terms
cd ~/matlab/CSIRO/Heat/Old_routines
[wbtp47, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
[wptb47, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
cd ~/matlab/CSIRO/Heat
[dhcdt47, lat_out, lon_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs2/data
%save BP2030_hctend_CT_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
%save HP6_hctend_CT_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
%save BP412_hctend_CT_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
%save LP7_hctend_CT_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
%save HP12_hctend_hpcs_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
save LP10_hctend_hpcs_new.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
ubtp47 vbtp47 wbtp47 uptb47 vptb47 wptb47 ...
dhcdt13 dhcdt47 hflx_reg tim lims nfrm lags pcs
% Get ubar_dHCprimedx and vbar_dHCprimedy
if 0;
cd ~/matlab/CSIRO/Heat/Old_routines
lev = 8:9;
[ubtp89, vbtp89, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
[uptb89, vptb89, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
% Get the wdtdz terms
cd ~/matlab/CSIRO/Heat/Old_routines
[wbtp89, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
[wptb89, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
cd ~/matlab/CSIRO/Heat
[dhcdt89, lat_out, lon_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
save LP9_hctend_89.mat ubtp89 vbtp89 wbtp89 uptb89 vptb89 wptb89 ...
dhcdt89 tim lims nfrm lags
end