function [x,y,z]=lavademi(phi,theta,t)
% une demi lava est construite comme suit:

% on prend la lava initiale
[x1,y1,z1]=lava0(phi,theta);

% on prend la lava au temps T
[x2,y2,z2]=lavac(phi,theta,t);

% le "morphing" entre la lava0 et la lava au temps T
% est donné par la fonction MORPH. Celle-ci est
% paramétrée par une fonction gaussienne du temps
% qui est étendue par périodicité pour permettre
% la continuité au départ et à la fin du cycle. 
% SIG donne la déviation standard.
sig=0.11;
pgt=exp(-(rem(t,1)./sig).^2)+exp(-((rem(t,1)-1)./sig).^2);

% le "morphing" a proprement parler
x=morph(x1,x2,pgt);
y=morph(y1,y2,pgt);
z=morph(z1,z2,pgt);
