//test of volume within volume to have
//single volume for meshing

//box 1
p1 = newp;
Point(p1) = {0.5,0.5,0.5};
l1[] = Extrude {0,0,1}{Point{p1};};
s1[] = Extrude {1,0,0}{Line{l1[1]};};
v1[] = Extrude {0,1,0}{Surface{s1[1]};};

//box2
p2 = newp;
Point(p2) = {0,0,0};
l2[] = Extrude{0,0,2}{Point{p2};};
s2[] = Extrude{2,0,0}{Line{l2[1]};};
v2[] = Extrude{0,2,0}{Surface{s2[1]};};

//line loops
ll1 = newll; 
Line Loop(ll1) = {l1[1],s1[0],s1[2],


//inner box surface loop
Printf ("isloop: &g", newsl);
ibsl[] = {s1[1],v1[0],v1[2],v1[3],v1[4],v1[5]};
isloop = newsl;
Surface Loop(isloop) = {ibsl[]};
//outer box surface loop
obsl[] = {s2[1],v2[0],v2[2],v2[3],v2[4],v2[5]};
osloop = newsl;
Surface Loop(osloop) = {obsl[]};

//inner box surface - for boundary conditions in ElmerFEM
isur = news;
Physical Surface(isur) ={ibsl[]};
//outer box surface
osur = news;
Physical Surface(osur) = {obsl[]};


//Printf("v1: %g,%g,%g,%g%g :: v2: %g,%g,%g,%g,%g", v1[0], v1[1],v1[2],v1[3],v1[4],v2[0],v2[1],v2[2],v2[3],v2[4]);

ivol = newv;
Volume(ivol) = {isloop};
gvol = newv;
Volume(gvol) = {osloop};
ipvol = newp;
Physical Volume(ipvol) = {ivol};
gpvol = newv;
Physical Volume(gpvol) = {gvol};
