Merge "asdf2.vtk";
CreateTopology;


// We can now define a compound line (resp. surface) for each discrete
// line (resp. surface) in the model
ll[] = Line "*";
For j In {0 : #ll[]-1}
Compound Line(newl) = ll[j];
EndFor
ss[] = Surface "*";
s = news;
For i In {0 : #ss[]-1}
Compound Surface(s+i) = ss[i];
EndFor
// And we can create the volume based on the new compound entities
Surface Loop(1) = {s : s + #ss[]-1};
//Volume(1) = {1};

Physical Surface(1) = {s : s + #ss[]-1};
Physical Volume(1) = 1;

Mesh.RemeshAlgorithm = 1; // (0) no split (1) automatic (2) automatic only with metis
Mesh.RemeshParametrization = 0; // (0) harmonic (1) conformal
//Mesh.Algorithm=6;
Mesh.SubdivisionAlgorithm=1;

Field[1] = MathEval;
Field[1].F="2.4";
Background Field =1;


//Geometry.HideCompounds = 0; // don’t hide the compound entities


//Mesh.CharacteristicLengthFromCurvature=3;

Hide {Volume{1}; }
//mesh;
