Documentation of reemerge1

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

Help text

  Try the COADS data

Cross-Reference Information

Listing of script reemerge1

clean
filin = 'sstcoadsanom4by6.19501997.nc';
varn = 'data';
tim = 1:576;
lims = [90 300 -45 75];
lev = 1;

cd /home/disk/tao/data/coads/sstanom4by6
[sst, lat, lon] = getnc2(filin, varn, lims, lev, tim);

%  Do the same with Kaplan:
filin = 'sstanomldeo1856jul2000.nc';
varn = 'ssta';
tim = ((1950-1856)*12)+[1:588];
lims = [[90 300]-360 -45 75];
lev = 1;

cd /home/disk/tao/data/sst_ldeo
[sst, lat, lon] = getnc2(filin, varn, lims, lev, tim);
lims = [lims(1:2)+360 lims(3:4)];

sst = annave(sst);

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

sst = reshape(sst, ntim, nlat*nlon);

nstart = 12;
nlag = 12;

c1 = repmat(NaN, [nstart, nlag, nlat*nlon]);

ind = 1:12:(ntim-12);
for i = 1:nstart; disp(i)
  ind1 = ind+i*(12/nstart)-1;
  for j = 1:nlag;
    ind2 = ind1+j*(12/nlag);
    for k = 1:nlat*nlon;
      c1(i,j,k) = corr_nan(sst(ind1,k), sst(ind2,k));
    end
  end
end

cd ~/matlab/NMC/Data
save KAPLAN_reemerge1.mat nstart nlag c1 lat lon lims
%save COADS_reemerge1.mat nstart nlag c1 lat lon lims



%  Plot the data
clean

cd ~/matlab/NMC/Data
load COADS_reemerge1.mat
load KAPLAN_reemerge1.mat

default_global; FRAME = lims;
dlat = mean(diff(YAX))/2;
dlon = mean(diff(XAX))/2;
color_trans; colormap(cmap);
nlat = length(YAX); nlon = length(XAX);

mon1 = 3;
figure(1); fl(1); clf;

for mon2 = 3:3:12;
  subplot(3,2,mon2/3);
  tem = reshape(c1(mon1, mon2, :), nlat, nlon);
 tem = smooth([1 2 1], smooth([1 2 1], tem')');
  [cc1, hh1] = mcont2(tem, [-1:.1:1], 'mollweid');
  hold on;
    hh2 = surfacem((YAX-dlat)*ones(1, nlon), ...
		   ones(nlat, 1)*(XAX-dlon)', tem);
  hold off;
  zdat = get(hh2, 'ZData'); set(hh2, 'ZData', zdat-20);
  caxis([-1 1]);
  title(['Start:  ' char(get_month2(mon1*(12/nstart))) ';  End:  ' ...
	 char(get_month2(mon1*(12/nstart)+mon2*(12/nlag)))]);
  dcmfill(-1);
  gridm on;
  axis_limits(100)
  colorbar2
end