layer.edp

load "msh3" 
mesh Th=square(3,3); 

mesh3 Th3=buildlayers(Th,3); 
savemesh(Th3,"lll.mesh");