Global Index (short | long) | Local contents | Local Index (short | long)
[dhcdt, lon_out, lat_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);
[dhcdt, lon_out, lat_out, tim] = get_dhcdt(pcs, lims, nfrm, tim, lev);
This function calls | This function is called by |
---|---|
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