Global Index (short | long) | Local contents | Local Index (short | long)
Now get tprime
This script calls | |
---|---|
clear tim = 101:550; ctlim = [110 300 -60 60]; lev = 1; cd /home/disk/hayes2/dvimont/csiro/annual_data filin = 'temp_half_L1_yr101-550.nc'; nc = netcdf(filin, 'nowrite'); depth = nc{'depth'}(:); latt = nc{'latitude'}(:); lont = nc{'longitude'}(:); [xk, yk] = keep_var(ctlim, lont, latt); temp = nc{'temp'}(:, 1, yk, xk); mv = nc{'temp'}.missing_value(:); nc = close(nc); temp(temp == mv) = NaN; latt = latt(yk); lont = lont(xk); temp = squeeze(temp); [ntim, nlat, nlon] = size(temp); temp = reshape(temp, ntim, nlat*nlon); dtdt = diff(temp) ./ (365 * 24 * 3600); ntim = ntim-1; % Load CPCs cd /home/disk/hayes2/dvimont/csiro/matlab_data/20Deg_Isotherm %load 20deg_CEOF_LP10_yr101-550.mat; tit = 'Lowpass Filtered Data ( > 10 Years )'; load 20deg_CEOF_HP8_yr101-550.mat; tit = 'Highpass Filtered Data ( < 8 Years )'; %load 20deg_CEOF_RAW_yr101-550.mat; tit = 'Unfiltered Data'; % Get regressions lag = 0; lg = lag*pi/180; lg2 = 1; num = 1; lind = 1; nfrm = 6; % Store all regressions under one variable j = sqrt(-1); timeseries = sqrt(2)*pcs(:,num)./std(pcs(:,num)); clear temtim treg for i = 1:nfrm wgt = conj(exp(j * ((i-1) * pi/(lg2*nfrm) + lg) )); temtim(:,i) = squeeze(real(wgt .* timeseries)); treg(i, :) = temtim(:,i)' * dtdt ./ ntim; end treg = reshape(treg, 6, nlat, nlon); % Plot the data get_global; FRAME = ctlim; XAX = lont; YAX = latt; figure(1); figure_orient; cint = 2; clev = [-50:cint:-cint cint:cint:50]; for i = 1:nfrm tem = 1e9 * squeeze(treg(i,:,:)); subplot(3,2,i); gcont(tem, clev); dc if i > 4; xlabel(['Contour Interval: ' num2str(cint) ' x 10^-^9 K s^-^1']); end end subplot(3,2,3); ylabel([tit ': Centered Time Derivative of T; '... 'Depth = ' num2str(depth(lev)/100) 'm']);