function mass_plot(scaling, delta, trans, dt, wdir)

if (nargin<1),                    scaling = [];    end
if (nargin<2) | isempty(delta),   delta = .035;    end
if (nargin<3) | isempty(trans),   trans = [.35 0]; end
if (nargin<4) | isempty(dt),      dt = .01;        end

if (nargin<5) | isempty(wdir),
  if isempty(dir('./packet0_r')),
    wdir = [];
  else
    wdir = '.';
  end
end

figure(1);

set_contour(.3);

for i = 0:8,
  n = 50*i;
  subplot(3,3,i+1);
  disp(sprintf('plot %d', n));
  pack_circle(packet(n, wdir), scaling, delta, trans);
  set(gca,'XTick',[],'YTick',[],'Box','on')
  title(sprintf('t = %f', n*dt));
  drawnow;
end

disp('print');
print -deps2 /var/tmp/prop_mosaic.eps
