title Intersecting Pipes

parameter d 3             c Mesh density parameter
          i1 [  1+%d*2]
          i2 [%i1+%d*4]
          i3 [  1+%d*4]
          i4 [  1+%d*2];

sd 1 cy 0 0 0 1 0 0 1.8   c Inner cylinder of large pipe
sd 2 cy 0 0 0 1 0 0 2     c Outer cylinder of large pipe
sd 3 cy 0 0 0 0 0 1 .5    c Outer cylinder of small pipe
sd 4 cy 0 0 0 0 0 1 .4    c Inner cylinder of small pipe
sd 5 plan 0 0 0 0 1 0     c Planer side
sd 6 plan 0 0 0 1 0 0     c Planer side
sd 7 plan 1.5 0 0 1 0 0   c Planer side
sd 8 plan 0 0 0 0 -1 1    c Planer side
sd 9 plan 0 0 2.5 0 0 1   c Top of small pipe

c global duplication of the parts
gct 3 rxz; ryz;rxz ryz;
lev 1 grep 0 1 2 3;;
pslv 1 

c Part 1 - Larger pipe
block 1 %i1 %i2;1 %i1 %i2;1 %i4;0 .5 1.5 0 .5 1.5 1.8 2  

dei 1 2;1 2; ;            c Remove 1 region

sfi 1 -2;1 -2; ;sd 3      c Project faces to surfaces
sfi ; ;-1;sd 1 
sfi ; ;-2;sd 2 
sfi ;-1; ;sd 5 
sfi -1; ; ;sd 6
sfi -3; ; ;sd 7
sfi ;-3; ;sd 8 

endpart

c Part 2 - smaller pipe
block 1 %i3;1 %i4;1 %i1 %i2;0 .5 .4 .5 1.8 2 3  

sfi ;-1; ;sd 4            c Project faces to surfaces
sfi ;-2; ;sd 3 
sfi ; ;-1;sd 1 
sfi ; ;-2;sd 2 
sfi -2; ; ;sd 5
sfi -1; ; ;sd 6
sfi ; ;-3;sd 9 

endpart
pplv 
merge
rx 20 ry 30 zf 1.8 tvv
