Documentation of get_hp_hctend_terms

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  Get ubar_dHCprimedx and vbar_dHCprimedy

Cross-Reference Information

Listing of script get_hp_hctend_terms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_detrend_L1-7_EOF_yr101-1000.mat
cd ~/matlab/CSIRO/Heat/Old_routines
nfrm = -10:10; lags = nfrm;
tim = 101:1000;
lev = 1:3;
lims = [28 355 -75 75];

   [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/Thesis/Data
save HP10_hctend_ALL.mat ubtp13 vbtp13 wbtp13 uptb13 vptb13 wptb13 ...
     ubtp47 vbtp47 wbtp47 uptb47 vptb47 wptb47 ...
     dhcdt13 dhcdt47 hflx_reg tim lims nfrm lags



%  Get ubar_dHCprimedx and vbar_dHCprimedy

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 HP10_hctend_89.mat ubtp89 vbtp89 wbtp89 uptb89 vptb89 wptb89 ...
     dhcdt89 tim lims nfrm lags