%   Compute profile of pendant drop in dimensional units.
    clear all
    clf reset
    sig=25;  % dyne/cm
    rho=1;   % g/cm^3
    g=980;   % cm/s^2
    L=sqrt(sig/(rho*g))
    load profile.dat;
    load shape.dat
    maxz=2;
    for i=1:length(profile(:,1))
        if profile(i,2)<2
            z(i)=L.*profile(i,2);
            r(i)=L.*profile(i,3);
        end
    end
    mr=-r;
    plot(r,z,'b-',mr,z,'b-',shape(:,1),shape(:,2),'r:',...
        -shape(:,1),shape(:,2),'r:')
    xlabel('cm')
    ylabel('cm')
    title('Shape of Pendant Drop')
    axis equal
%    shape=[r' z'];
%    save shape.dat shape -ascii