Point(1) = {-0.002, 0.00, 0.000, 0.0005};
Point(2) = {0.0, 0.002, 0.000, 0.0005};
Point(3) = {0.002, 0.00, 0.000, 0.0005};
Point(4) = {0.0, -0.002, 0.000, 0.0005};
Point(9) = {-0.001, 0.00, 0.00, 0.0005};
Point(10) = {0.0, 0.001, 0.00, 0.0005};
Point(11) = {0.001, 0.00, 0.00, 0.0005};
Point(12) = {0.0, -0.001, 0.00, 0.0005};

Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line Loop(13) = {1,2,3,4};
Line(19) = {9,10};
Line(20) = {10,11};
Line(21) = {11,12};
Line(22) = {12,9};
Line Loop(23) = {19,20,21,22};
Plane Surface(24) = {13,23};
Plane Surface(25) = {23};

e1[] = Extrude {0,0,0.004} {
  Surface{24, 25}; Recombine; Layers{10,1};
};

e2[] = Extrude {0,0,0.004} {
  Surface{e1[2]}; Recombine; Layers{10,1};
};

//Physical Line used to add Convection Boundary Condition
Physical Surface(1000) = {24,25};

//Physical Line used to Add Flux Boundary Condition
Physical Surface(1001) = {e2[0]};

//Physical Surface used to represent the Bar
Physical Volume(100) = {e1[1], e1[3]};
Physical Volume(101) = {e2[1]}; 

