load FP_dens_clim
load Track_dens_clim_8100
load ACE_dens_clim_8100

titles = {'ACE','TD','FP',};

z(:,1) = acedens_cl(:);
z(:,3) = fpdens_cl(:);
z(:,2) = trdens_cl(:);

for ii=1:3
  for jj=1:3
    subplot(3,3,sub2ind([3,3],ii,jj))
    h=plot(z(:,ii),z(:,jj),'k.');
    set(h,'markersize',3)
    xlabel(titles{ii})
    ylabel(titles{jj})
    if ii==jj
      axis('equal')
    end
    axis('square')
  end
end
