// #################################################################
// #                       Definition of Parameters                #
// #################################################################
x_min = 0;
x_max = 1;
y_min = 0;
y_max = .5;
elementSize = .5;

// #################################################################
// #               Definition of  main points                      #
// #################################################################
Point (1) = {x_min, y_min,        0};
Point (2) = {x_max, y_min,        0};
Point (3) = {x_min, y_max,        0};
Point (4) = {x_max, y_max,        0};

// #################################################################
// #                    Define main lines                          #
// #################################################################
Line(1) = {1, 2};
Line(2) = {2, 4};
Line(3) = {4, 3};
Line(4) = {3, 1};

// #################################################################
// #                Define layer (material)                        #
// #################################################################
Line Loop (1) = {1, 2, 3, 4};
Plane Surface (1) = {1};

// #################################################################
// #                    Add layer tag                              #
// #################################################################
Physical Surface ('Material_1', 100) = {1};

// #################################################################
// #             Define mesh size for each layer                   #
// #################################################################
// First material
Characteristic Length {1, 2, 3, 4} = elementSize;

// #################################################################
// #                          Set periodic mesh                    #
// #################################################################
Periodic Line {4} = {2};
Periodic Line {1} = {3};
