Documentation of get_dhcdt

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

Function Synopsis

[dhcdt, lon_out, lat_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);

Help text


  [dhcdt, lon_out, lat_out, tim] = get_dhcdt(pcs, lims, nfrm, tim, lev);

Cross-Reference Information

This function calls	This function is called by
getheat regress_ceof regress_eof	get_hctend_terms get_hp_hctend_terms get_hctend_terms plot_30s30n_terms get_hctend_terms plot_30s30n_terms get_hctend_terms get_hp_hctend_terms plot_15n65n_terms plot_30s30n_EOF_terms

Listing of function get_dhcdt

function [dhcdt, lon_out, lat_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);

cdtem = ['cd ' eval('pwd')];

if size(tim, 2) == 1; tim = tim'; end;

if max(tim) ~= 1000;
  maxout = 0;
  tim2 = [min(tim)-1 tim max(tim)+1];
else
  maxout = 1;
  tim2 = [min(tim)-1 tim];
end

[hc, lat_out, lon_out] = getheat(lims, lev, tim2);

[ntim, nlat, nlon] = size(hc);

dhcdt = (hc(3:ntim,:,:) - hc(1:(ntim-2),:,:)) ./ (3600*24*365*2);

tim2 = tim2(2:(ntim-1));
timind = tim2 - min(tim2) + 1;

if isreal(pcs);
  dhcdt = regress_eof(dhcdt(timind,:,:), pcs(timind,:), nfrm);
else
  dhcdt = regress_ceof(dhcdt(timind,:,:), pcs(timind,:), nfrm);
end