Global Index (short | long) | Local contents | Local Index (short | long)
The following files are in the directory below: filin = 'temp_A_L1_1000_years.nc' % SST filin = 'temp_A_L5_1000_years.nc' % Layer 5 Temp. filin = 'temp_A_L1-5_1000_years.nc' % Averaged Pot. Temp., Top 160m
This script calls | |
---|---|
clear filin = 'temp_A_L1-10_1000_years.nc' % Averaged Pot. Temp., Top 620m cd /home/disk/hayes2/dvimont/csiro/data nc = netcdf(filin, 'nowrite'); lat = nc{'latitude'}(:); lon = nc{'longitude'}(:); % ctlim = [130 285 -50 65]; % ctlim = [180 270 -6 6]; % ctlim = [110 275 -75 65]; ctlim = [-.1 360 -90 90]; [xk, yk] = keep_var(ctlim, lon, lat); % temp = nc{'temp'}(:,1,yk,xk); temp = nc{'temp'}(:,yk,xk); mv = nc{'temp'}.missing_value(:); nc = close(nc); temp = squeeze(temp); [ntim, nlat, nlon] = size(temp); temp = reshape(temp, ntim, nlat*nlon); for i = 1:ntim temp(i,find(temp(i,:) == mv)) = NaN * ones(size(find(temp(i,:) == mv))); end; clim = mean(temp); kp = find(~isnan(clim)); temp = temp(:,kp); temp = detrend(temp); lat = lat(yk); lon = lon(xk); get_global default_global; FRAME = [0 360 -90 90]; cd /home/disk/tao/dvimont/matlab/CSIRO/Data load butter_4.5_ctstar.mat ctlow = detrend(ctlow); ctlow = (ctlow - mean(ctlow)) / std(ctlow); pos = find(ctlow > 0); neg = find(ctlow < 0); pospat = NaN * ones(1, nlat*nlon); negpat = NaN * ones(1, nlat*nlon); pospat(kp) = ctlow(pos)' * temp(pos, :) ./ sum(ctlow(pos)); negpat(kp) = ctlow(neg)' * temp(neg, :) ./ sum(ctlow(neg)); pospat = reshape(pospat, nlat, nlon); negpat = reshape(negpat, nlat, nlon); pospat(1,1) = -.4; pospat(1,2) = .4; negpat(1,1) = -.4; negpat(1,2) = .4; figure(1); figure_orient; sp(1) gshade(pospat, [-.4:.05:.4]); dc; colorbar sp(2) gshade(negpat, [-.4:.05:.4]); dc; colorbar