load "PETSc" macro dimension()3// EOM include "macro_ddm.idp" load "iovtk" load "Curvature" load "isoline" load "msh3" load "medit" bool importMesh = true; string meshFileName = "mesh/"; string outputFileName = "output/"; int NN=1000, saveEach = 1; real muf=10, c1=1.e7, rhof=100, rhos=100, g=-0, rhod=rhos-rhof; real T=2, t=0, dt=T/NN, penal=0,H=20; int[int] orderOut = [1,1,1,1],orderOuts = [1,1,1,1,1,1]; func Pk=[P2,P2,P2,P1]; mesh3 Th, Ths, Ths0; if (importMesh) { Th=readmesh3(meshFileName+"Th_cylinder.mesh"); Ths=readmesh3(meshFileName+"Ths_cylinder.mesh"); } else {/* Th = square(160,40, [4*x,y]); //plot(Th, wait=1); Ths = square(1,40, [0.0212*x+2-0.0106,0.8*y]); //plot(Ths, Th, wait=1); savemesh(Th, "Th.msh"); savemesh(Ths, "Ths.msh"); */ cout<<"please import mesh..."<