function [x,y,z]=lavac(phi,theta,t)
global tc
% Forme lava dépendant du temps

% comme les pièces partagent une arête jusqu'au temps critique
% elles doivent être déplacées verticalement. La distance verticale
% est constante jusqu'au temps critique, puis décroit exponentiellement
% après le temps critique
gt=zeros(size(t));
i=find(t<tc);
gt=exp(-10*(t-tc));
gt(i)=t(i)./(t(i)+eps);

% La lava au temps t
rc0=rc(0,tc);
[x,y,z]=sphere2(rc(theta,t),phi,theta);
z=z-rc0*gt;
