function [x,y,z]=lava0(phi,theta)
% la lava au temps t=0

% le rayon initial est donné par:
r0=0.8-0.2*legendreP(2,cos((pi+theta)./2));

% coordonnées de la lava
[x,y,z]=sphere2(r0,phi,(pi+theta)./2);
