Documentation of getheat

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

Function Synopsis

[heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn);

Help text


 [heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn);

 lev = level indices over which to get the heat content 
 tim = time period
 lims = lat-lon frame
 varn = variable to be averaged

Cross-Reference Information

This function calls	This function is called by
getll getnc	heat_ceof plot_60s60n_terms heat_ceof eq_tend_terms eq_tend_terms2 lagged_correlation lagged_regressions regress_velocity wbar_dtpdz_hc0-80 plot_lp_pac_regressions plot_regressions heat_ceof eq_tend_terms eq_tend_terms2 tendency_terms calc_eq_tend_terms eq_tend_terms eq_tend_terms2 heat_ceof plot_30s30n_EOF_regressions plot_30s30n_regressions eq_tend_terms calc_clim calc_eq_tend_terms complex_eof_heat_content composite corr_ct_mlcurl corr_indices cospectrum eq_ind_corr eq_tend_terms eq_tend_terms2 eq_tend_terms3 get_dhcdt gu_philander gu_philander_thesis HC_clim hc_eof hc_eofs hc_sst_lagged_regression hc_sst_lagged_regression_color heat_ceof heat_content_clim hoffmoeller_corr_plot hoffmoeller_plot HP10_eq_tend_terms HP10_LP9_wreg HP_HC_evolution_30s30n lag_regress_on_ct lagged_1000yr_regression lagged_correlation lagged_regression lagged_regressions look_hc_terms LP9_eq_tend_terms LP_HC_evolution_30s30n make_regressions movie_hc_25yrlp ocean_current_clim ocean_vel plot_15n65n_regressions plot_30s30n_EOF_regressions plot_30s30n_EOF_terms plot_30s30n_regressions plot_30s30n_terms plot_60s60n_terms plot_EOF1_EOF2_0-270mHC plot_eq_currents plot_lag_regress2 plot_lag_regressions plot_lag_regressions_color plot_lp_pac_regressions plot_pac_regressions plot_reg_HEAT_SFC plot_regressions plot_regressions_HP_LP plot_regressions_HP_LP_thesis plot_variance_explained regress_sv regress_taux_hanom sst_complex_eof sst_tau_hc_color_fig sst_tau_hc_fig tendency_terms tendency_terms_lreg time_freq_anal wbar_dtpdz_hc0-80 write_netcdf

This function calls

This function is called by

getll
getnc

Listing of function getheat

function [heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn);

cdback = eval('cd');
cd /home/disk/hayes2/dvimont/csiro/data

if nargin < 4;  varn = 'temp';  end;
varn
temp = getnc(varn, lims, lev, tim);
[lat, lon, depth] = getll(varn, lims);

[ntim, nlev, nlat, nlon] = size(temp);

d1 = depth(1);
middepth = [0];
for i = 1:(length(depth)-1);
  middepth = [middepth; depth(i)+d1];
  d1 = depth(i+1) - depth(i) - d1;
end
middepth = [middepth; depth(length(depth))+d1];
wgt = diff(middepth);

temp = shiftdim(temp, 1);
temp = reshape(temp, nlev, nlat*nlon*ntim);

heat = wgt(lev)' * temp;
heat = shiftdim(reshape(heat, nlat, nlon, ntim), 2);

%  If necessary, Convert heat to units J / m^2
if strcmp(varn, 'temp');
  heat = heat * 1000 * 4.218e3;
else
  heat = heat ./ sum(wgt(lev));
end

eval(['cd ' cdback]);