Include "geometry.geo";

Physical Surface ("R_Boundary" ) ={21,20,19,22,18,17,16,15};
Physical Volume("P") = {1};
Physical Volume("R") = {2};
Physical Volume("H") = {3};