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

%

if (nargin<1),                    scaling = .19;    end
if (nargin<2) | isempty(delta),   delta = .046;     end
if (nargin<3) | isempty(trans),   trans = [.322 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);

% clear the energy contour
%global energy_contour
%energy_contour=[];
set_contour(.3);

for n = 0:25:1500,
  disp(sprintf('plot %d', n));
  pack_circle(packet(n, wdir), scaling, delta, trans);
  set(gca,'YLim',[-1 1],'Visible','off');
  disp(sprintf('print %d', n));
  print('-deps2', sprintf('prop_%d.eps', n));
end

%delete(1);
