// Vamos a crear una malla para GMSH
// La malla va a representar un campo con un único árbol

rt = 5; // Radio del árbol
xt = 7.5; // Posición x del árbol
yt = 22.5; // Posición y del árbol

tm = 1; // Tamaño de la malla en el contorno
tmr = .1; // Tamaño de la malla en la interfaz árbol / pasto

Point(1) = {0, 0, 0, tm};
Point(2) = {30, 0, 0, tm};
Point(3) = {30, 30, 0, tm};
Point(4) = {0, 30, 0, tm};

Point(5) = {xt, yt, 0, tm};
Point(6) = {xt+rt, yt, 0, tmr};
Point(7) = {xt-rt, yt, 0, tmr};


Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};

Circle(5) = {6, 5, 7};
Circle(6) = {7, 5, 6};

Curve Loop(1) = {1, 2, 3, 4};
Curve Loop(2) = {5, 6};

Plane Surface(1) = {1};
Curve{5,6} In Surface{1};
Point{5} In Surface{1};
Periodic Curve {1} = {3};
Periodic Curve {2} = {4};

Physical Line("Sur") = {1};
Physical Line("Este") = {2};
Physical Line("Norte") = {3};
Physical Line("Oeste") = {4};
Physical Line("Arbol") = {5,6};

Mesh.Algorithm = 5;
Mesh 2;

Mesh.SurfaceFaces = 1;
Mesh.Points = 1;

Save "ArbolUnico.msh";