c contains necessary geometry
INCLUDE dblbtrf.geo

c governs density
  para d1 5
       d2 5
       d3 10
       d4 3
       d5 2
       d6 3
       d7 2
       d8 5;

c ********************************************
c *         PART I - The Elbow Pipe          *
c ********************************************
c 
c create a block
c --------------
  block 1 [1+%d1] [1+%d1+%d2] [1+2*%d1+%d2];
        1 [1+%d3] [1+%d3+%d4] [1+2*%d3+%d4];
        1 [1+%d5] [1+%d5+%d6] [1+%d5+%d6+%d1];
        0 0 0 0 ; 0 0 0 0 ; 0 0 0 0 ;

c delete regions
dei 1 2 0 3 4;    ; 3 4; 
dei          ; 2 3; 1 2;

c interactively reposition nodes
mbi ; -1;; xy 2.99784 -3.00015
mbi ; -2;; xy 2.75673 -0.338176
mbi ; -3;; xy 2.58395 -0.184205
mbi ; -4;; xy -0.639069e-02 0.129426e-03

c interactively place edges on curves
c     region      curve number
c     ------      ------------
cur 2 3 4 2 4 4      2
cur 2 2 4 2 3 4      2
cur 2 1 4 2 2 4      2

cur 1 3 3 1 4 3      2
cur 1 2 3 1 3 3      2
cur 1 1 3 1 2 3      2

cur 3 3 4 3 4 4      1
cur 3 2 4 3 3 4      1
cur 3 1 4 3 2 4      1

cur 4 1 3 4 2 3      1
cur 4 2 3 4 3 3      1
cur 4 3 3 4 4 3      1

cur 1 3 1 1 4 1      9
cur 1 1 1 1 2 1      9

cur 4 3 1 4 4 1     10
cur 4 1 1 4 2 1     10

cur 2 3 1 2 4 1      6
cur 2 1 1 2 2 1      6

cur 3 3 1 3 4 1      5
cur 3 1 1 3 2 1      5

cur 1 3 2 1 4 2      8
cur 1 1 2 1 2 2      8

cur 4 3 2 4 4 2      7
cur 4 1 2 4 2 2      7

cur 1 2 2 4 2 2     13
cur 1 3 2 4 3 2     14

c linear interpolation
lini -2 0 -3; -1 0 -4; 1 3;

c projection of regions onto surfaces
sfi     2  3   ;        ;       -3   ;sd 6
sfi     2 -3   ;        ;  1     3   ;sd 4
sfi    -2  3   ;        ;  1     3   ;sd 5
sfi -1       -4;        ;  1       -4;sd 1
sfi          -4;  2    3;    -2      ;sd 2
sfi          -4; -2 0 -3;  1  2      ;sd 2
sfi            ;        ; -1         ;sd 3
sfi  1     3   ; -2   -3;  1 -2      ;sd 7

mt 1 2 1 4 3 2 1
mate 2

c create the block boundaries
bb 1 3 1 4 3 2    1;
bb 1 2 1 4 2 2    2;
bb 1 2 2 4 3 2    3;

endpart

c ********************************************
c *       PART II - The Cylindrical Pipe     *
c ********************************************
c
block 1 [1+2*%d1+%d2] [1+2*%d1+%d2+%d8]; 
      1 [1+%d5] [1+%d5+%d7];
      1 [1+%d7] [1+%d7+%d4] [1+2*%d7+%d4];
      0 0 0; 0 0 0; 0 0 0 0;
dei ; 2 3; 1 2 0 3 4; 

c interactively reposition nodes
mbi -3;;; xy 5.15668 0.245911e-02
mbi -2;;; xy 2.92854 -0.153116e-02
mbi -1;;; xy 2.43177 -0.550366

mbi -1;;; xyz -0.242226e-01 0.241145e-01 0.805417e-03
mbi 1 2;; 1 2; xyz -0.449290e-01 0.447270e-01 0.192499e-04
mbi 1 2;; 3 4; xyz 0.472760e-01 -0.470627e-01 0.266617e-03
mbi 1 2; 2 3;; xyz 0.787020e-03 -0.731915e-03 0.515751e-01

mbi 1 2;; 3 4; xyz 0.385547e-02 -0.383824e-02 0.182632e-03
mbi 1 2;; 1 2; xyz -0.374436e-02 0.372767e-02 0.289891e-03
mbi 1 2;; 1 2; xyz -0.365448e-02 0.362697e-02 -0.253223e-03

mbi -3;; 3 4; xyz 0.476837e-06 -0.549551e-01 0.507498e-03
mbi -3;; 1 2; yz 0.480877e-01 0.548333e-03
mbi -3; 2 3;; xyz -0.762939e-05 -0.246647e-03 0.513116e-01

mbi -3;; 3 4; xyz 0.476837e-06 0.455779e-02 0.428945e-03

c attach parts of the part to previously defined block boundaries
bb 1 1 1 2 2 1 1 ; 
bb 1 1 4 2 2 4 2 ; 
bb 1 3 2 2 3 3 3 ; 

c final projections
sfi 2 3; 1 -3; -1 -4;sd 2
sfi -1;;;sd 1


endpart merge
rest rx 20 ry 30 tvv
