c = 50e-3;

// Important points
p1 = newp; Point (p1) = {  0,  0, 0, c };				// Main channel valve
p2 = newp; Point (p2) = {  0,  0, 13.5e-3, c };			// Valve
p3 = newp; Point (p3) = {  0,  22e-3, 0, c };			// Prerestriction
p4 = newp; Point (p4) = {  0, -12.5e-3, 0, c }; 		// Main channel sensor
p5 = newp; Point (p5) = {  0, -12.5e-3, 13.5e-3, c }; 	// Sensor
p6 = newp; Point (p6) = {  0, -56.3e-3, 0, c }; 		// First t-junction
p7 = newp; Point (p7) = {  0, -76.3e-3, 0, c }; 		// Second t-junction
p8 = newp; Point (p8) = {  0, -96.3e-3, 0, c }; 		// Corner
p9 = newp; Point (p9) = { 65.5e-3, -56.3e-3, 0, c }; 	// Branch 1
p10 = newp; Point (p10) = { 65.5e-3, -76.3e-3, 0, c }; 	// Branch 2
p11 = newp; Point (p11) = { 65.5e-3, -96.3e-3, 0, c }; 	// Branch 3

// Lines
l1 = newreg; Line (l1) = { p2, p1 }; 					// Valve channel
l2 = newreg; Line (l2) = { p1, p3 }; 					// Main channel
l3 = newreg; Line (l3) = { p1, p4 };
l4 = newreg; Line (l4) = { p4, p6 };
l5 = newreg; Line (l5) = { p6, p7 };
l6 = newreg; Line (l6) = { p7, p8 };
l7 = newreg; Line (l7) = { p4, p5 }; 					// Sensor channel
l8 = newreg; Line (l8) = { p6, p9 }; 					// Branch 1
l9 = newreg; Line (l9) = { p7, p10 }; 					// Branch 2
l10 = newreg; Line (l10) = { p8, p11 };					// Branch 3

// Domains
Physical Point ( "valve" ) = { p2 };
Physical Point ( "sensor" ) = { p5 };
Physical Line ( "supply" ) = { l2, l3, l4, l5, l6, l8, l9, l10 };
Physical Line ( "taps" ) = { l2, l7 };


