pycad/examples/lava_dom.geo

lc_border = 0.4;
lc_surf = 0.4;

N=8;
DTheta = 2*Pi/N;
base_z = 0.0;
top_z = 4.0;
surf_top_z = 7.0;
Theta = 0;
rad_conduit = 1.5;
rad_surf = 3.0;

CenterPoint1 = newp;
Point(CenterPoint1) = { 0.0,0.0,base_z,lc_border};

// Conduit bottom Circle
For I In {1:N}
   PointList[I] = newp;
   Theta+=DTheta;
   Point( PointList[ I] ) = { rad_conduit*Cos(Theta)  
,rad_conduit*Sin(Theta), base_z, lc_border};
EndFor

// Conduit Top Circle
CenterPoint2 = newp;
Point(CenterPoint2) = {0.0,0.0,top_z,lc_border};
Theta = 0.0;
For I In {1:N}
   PointList2[I] = newp;
   Theta+=DTheta;
   Point( PointList2[ I] ) = { rad_conduit*Cos(Theta)  
,rad_conduit*Sin(Theta), top_z, lc_border};
EndFor

// Bottom Surface Lines
For I In {1:N}
   LineList[I] = newl;
   Line( LineList[I] ) = { PointList[I], CenterPoint1 };
EndFor
For I In {1:N}
   CircleList[I] = newc;
   Circle( CircleList[I] ) = { PointList[I], CenterPoint1,  
PointList[(I%N) + 1] };
EndFor


// Top Surface Lines
//For I In {1:N}
//  LineList2[I] = newl;
//  Line( LineList2[I] ) = { PointList2[I], CenterPoint2 };
//EndFor
For I In {1:N}
   CircleList2[I] = newc;
   Circle( CircleList2[I] ) = { PointList2[I], CenterPoint2,  
PointList2[(I%N) + 1] };
EndFor


// Lines Connecting Top and Bottom Surfaces
For I In {1:N}
   Line(newl) = { PointList[I], PointList2[I] };
EndFor


// ******   DOME ****************
// ******************************

// Surf Bot
Theta = 0.0;
For I In {1:N}
   SurfPoint1[I] =  newp;
   Theta += DTheta;
   Point( SurfPoint1[I] ) = { rad_surf*Cos(Theta)  
,rad_surf*Sin(Theta), top_z, lc_surf };
EndFor


// Surf Top
CenterPoint3 = newp;
Point(CenterPoint3) = {0.0,0.0,surf_top_z,lc_border};
Theta = 0.0;
For I In {1:N}
   SurfPoint2[I] = newp;
   Theta+=DTheta;
   Point( SurfPoint2[ I] ) = { rad_surf*Cos(Theta)  
,rad_surf*Sin(Theta), surf_top_z, lc_surf};
EndFor

// Lines Connecting Top and bottom
For I In {1:N}
   Line(newl) = { SurfPoint1[I], SurfPoint2[I] };
EndFor

// Lines connecting outer to inner circle on surface

For I In {1:N}
   Line(newl) = { PointList2[I], SurfPoint1[I]};
EndFor


// Bottom Surface Lines
For I In {1:N}
   CircleList4[I] = newc;
   Circle(CircleList4[I] ) = { SurfPoint1[I], CenterPoint2,  
SurfPoint1[(I%N) + 1] };
EndFor


// Top Surface Lid
For I In {1:N}
   LineList3[I] = newl;
   Line( LineList3[I] ) = { SurfPoint2[I], CenterPoint3 };
EndFor
For I In {1:N}
   CircleList3[I] = newc;
   Circle( CircleList3[I] ) = { SurfPoint2[I], CenterPoint3,  
SurfPoint2[(I%N) + 1] };
EndFor

// ** BEGIN **
// conduit Bottom Surface
Line Loop(73) = {4,-5,-12};
Plane Surface(74) = {73};
Line Loop(75) = {3,-4,-11};
Plane Surface(76) = {75};
Line Loop(77) = {10,3,-2};
Plane Surface(78) = {77};
Line Loop(79) = {9,2,-1};
Plane Surface(80) = {79};
Line Loop(81) = {16,1,-8};
Plane Surface(82) = {81};
Line Loop(83) = {15,8,-7};
Plane Surface(84) = {83};
Line Loop(85) = {14,7,-6};
Plane Surface(86) = {85};
Line Loop(87) = {5,-6,-13};
Plane Surface(88) = {87};


Physical Surface(153) = {82,80,78,76,74,88,86,84};
//  ** END **

// ** BEGIN **
// Volcano Top Surface
Line Loop(89) = {54,-47,-22,46};
Plane Surface(90) = {89};
Line Loop(91) = {48,-55,-47,23};
Plane Surface(92) = {91};
Line Loop(93) = {24,41,-56,-48};
Plane Surface(94) = {93};
Line Loop(95) = {49,-42,-17,41};
Plane Surface(96) = {95};
Line Loop(97) = {50,-43,-18,42};
Plane Surface(98) = {97};
Line Loop(99) = {51,-44,-19,43};
Plane Surface(100) = {99};
Line Loop(101) = {45,-52,-44,20};
Plane Surface(102) = {101};
Line Loop(103) = {53,-46,-21,45};
Plane Surface(104) = {103};


Physical Surface(154) = {90,92,94,96,104,102,100,116,98};
// ** END **

// ** BEGIN **
// Volcano Free Surface Top
Line Loop(105) = {70,63,-62};
Plane Surface(106) = {105};
Line Loop(107) = {62,-61,69};
Plane Surface(108) = {107};
Line Loop(109) = {61,-60,68};
Plane Surface(110) = {109};
Line Loop(111) = {60,-59,67};
Plane Surface(112) = {111};
Line Loop(113) = {58,-59,-66};
Plane Surface(114) = {113};
Line Loop(115) = {57,-58,-65};
Plane Surface(116) = {115};
Line Loop(117) = {72,57,-64};
Plane Surface(118) = {117};
Line Loop(119) = {63,-64,-71};
Plane Surface(120) = {119};


Physical Surface(155) = {106,110,112,114,116,118,120,108};
// ** END **


// ** BEGIN **
// Conduit Outer Surface
Line Loop(121) = {13,30,-21,-29};
Ruled Surface(122) = {121};
Line Loop(123) = {29,-20,-28,12};
Ruled Surface(124) = {123};
Line Loop(125) = {28,-19,-27,11};
Ruled Surface(126) = {125};
Line Loop(127) = {27,-18,-26,10};
Ruled Surface(128) = {127};
Line Loop(129) = {26,-17,-25,9};
Ruled Surface(130) = {129};
Line Loop(131) = {25,-24,-32,16};
Ruled Surface(132) = {131};
Line Loop(133) = {32,-23,-31,15};
Ruled Surface(134) = {133};
Line Loop(135) = {30,22,-31,-14};
Ruled Surface(136) = {135};


Physical Surface(156) = {134,136,122,126,128,130,132,124};
// ** END **


// ** BEGIN **
// Outer Free Surface Cylinder
Line Loop(137) = {53,38,-69,-37};
Ruled Surface(138) = {137};
Line Loop(139) = {52,37,-68,-36};
Ruled Surface(140) = {139};
Line Loop(141) = {51,36,-67,-35};
Ruled Surface(142) = {141};
Line Loop(143) = {50,35,-66,-34};
Ruled Surface(144) = {143};
Line Loop(145) = {49,34,-65,-33};
Ruled Surface(146) = {145};
Line Loop(147) = {56,33,-72,-40};
Ruled Surface(148) = {147};
Line Loop(149) = {55,40,-71,-39};
Ruled Surface(150) = {149};
Line Loop(151) = {39,-70,-38,54};
Ruled Surface(152) = {151};

Physical Surface(157) = {138,148,144,142,140,150,146,152};
// ** End **


Surface Loop(158) =  
{134,132,130,128,126,124,122,88,74,76,-78,-80,-82,-84,-86,-136,-90,152,150,92,94,-96,146,144,-98,-100,142,140,102,-104,138,108,106,-120,118,148,-116,-114,112,110};

Volume(159) = {158};
Physical Volume(160)  = { 159};
1	lc_border = 0.4;
2	lc_surf = 0.4;
3
4	N=8;
5	DTheta = 2*Pi/N;
6	base_z = 0.0;
7	top_z = 4.0;
8	surf_top_z = 7.0;
9	Theta = 0;
10	rad_conduit = 1.5;
11	rad_surf = 3.0;
12
13	CenterPoint1 = newp;
14	Point(CenterPoint1) = { 0.0,0.0,base_z,lc_border};
15
16	// Conduit bottom Circle
17	For I In {1:N}
18	PointList[I] = newp;
19	Theta+=DTheta;
20	Point( PointList[ I] ) = { rad_conduit*Cos(Theta)
21	,rad_conduit*Sin(Theta), base_z, lc_border};
22	EndFor
23
24	// Conduit Top Circle
25	CenterPoint2 = newp;
26	Point(CenterPoint2) = {0.0,0.0,top_z,lc_border};
27	Theta = 0.0;
28	For I In {1:N}
29	PointList2[I] = newp;
30	Theta+=DTheta;
31	Point( PointList2[ I] ) = { rad_conduit*Cos(Theta)
32	,rad_conduit*Sin(Theta), top_z, lc_border};
33	EndFor
34
35	// Bottom Surface Lines
36	For I In {1:N}
37	LineList[I] = newl;
38	Line( LineList[I] ) = { PointList[I], CenterPoint1 };
39	EndFor
40	For I In {1:N}
41	CircleList[I] = newc;
42	Circle( CircleList[I] ) = { PointList[I], CenterPoint1,
43	PointList[(I%N) + 1] };
44	EndFor
45
46
47	// Top Surface Lines
48	//For I In {1:N}
49	// LineList2[I] = newl;
50	// Line( LineList2[I] ) = { PointList2[I], CenterPoint2 };
51	//EndFor
52	For I In {1:N}
53	CircleList2[I] = newc;
54	Circle( CircleList2[I] ) = { PointList2[I], CenterPoint2,
55	PointList2[(I%N) + 1] };
56	EndFor
57
58
59
60	// Lines Connecting Top and Bottom Surfaces
61	For I In {1:N}
62	Line(newl) = { PointList[I], PointList2[I] };
63	EndFor
64
65
66
67	// **** DOME **************
68	// ******************************
69
70	// Surf Bot
71	Theta = 0.0;
72	For I In {1:N}
73	SurfPoint1[I] = newp;
74	Theta += DTheta;
75	Point( SurfPoint1[I] ) = { rad_surf*Cos(Theta)
76	,rad_surf*Sin(Theta), top_z, lc_surf };
77	EndFor
78
79
80
81	// Surf Top
82	CenterPoint3 = newp;
83	Point(CenterPoint3) = {0.0,0.0,surf_top_z,lc_border};
84	Theta = 0.0;
85	For I In {1:N}
86	SurfPoint2[I] = newp;
87	Theta+=DTheta;
88	Point( SurfPoint2[ I] ) = { rad_surf*Cos(Theta)
89	,rad_surf*Sin(Theta), surf_top_z, lc_surf};
90	EndFor
91
92	// Lines Connecting Top and bottom
93	For I In {1:N}
94	Line(newl) = { SurfPoint1[I], SurfPoint2[I] };
95	EndFor
96
97	// Lines connecting outer to inner circle on surface
98
99	For I In {1:N}
100	Line(newl) = { PointList2[I], SurfPoint1[I]};
101	EndFor
102
103
104
105	// Bottom Surface Lines
106	For I In {1:N}
107	CircleList4[I] = newc;
108	Circle(CircleList4[I] ) = { SurfPoint1[I], CenterPoint2,
109	SurfPoint1[(I%N) + 1] };
110	EndFor
111
112
113
114	// Top Surface Lid
115	For I In {1:N}
116	LineList3[I] = newl;
117	Line( LineList3[I] ) = { SurfPoint2[I], CenterPoint3 };
118	EndFor
119	For I In {1:N}
120	CircleList3[I] = newc;
121	Circle( CircleList3[I] ) = { SurfPoint2[I], CenterPoint3,
122	SurfPoint2[(I%N) + 1] };
123	EndFor
124
125	// BEGIN
126	// conduit Bottom Surface
127	Line Loop(73) = {4,-5,-12};
128	Plane Surface(74) = {73};
129	Line Loop(75) = {3,-4,-11};
130	Plane Surface(76) = {75};
131	Line Loop(77) = {10,3,-2};
132	Plane Surface(78) = {77};
133	Line Loop(79) = {9,2,-1};
134	Plane Surface(80) = {79};
135	Line Loop(81) = {16,1,-8};
136	Plane Surface(82) = {81};
137	Line Loop(83) = {15,8,-7};
138	Plane Surface(84) = {83};
139	Line Loop(85) = {14,7,-6};
140	Plane Surface(86) = {85};
141	Line Loop(87) = {5,-6,-13};
142	Plane Surface(88) = {87};
143
144
145	Physical Surface(153) = {82,80,78,76,74,88,86,84};
146	// END
147
148	// BEGIN
149	// Volcano Top Surface
150	Line Loop(89) = {54,-47,-22,46};
151	Plane Surface(90) = {89};
152	Line Loop(91) = {48,-55,-47,23};
153	Plane Surface(92) = {91};
154	Line Loop(93) = {24,41,-56,-48};
155	Plane Surface(94) = {93};
156	Line Loop(95) = {49,-42,-17,41};
157	Plane Surface(96) = {95};
158	Line Loop(97) = {50,-43,-18,42};
159	Plane Surface(98) = {97};
160	Line Loop(99) = {51,-44,-19,43};
161	Plane Surface(100) = {99};
162	Line Loop(101) = {45,-52,-44,20};
163	Plane Surface(102) = {101};
164	Line Loop(103) = {53,-46,-21,45};
165	Plane Surface(104) = {103};
166
167
168	Physical Surface(154) = {90,92,94,96,104,102,100,116,98};
169	// END
170
171	// BEGIN
172	// Volcano Free Surface Top
173	Line Loop(105) = {70,63,-62};
174	Plane Surface(106) = {105};
175	Line Loop(107) = {62,-61,69};
176	Plane Surface(108) = {107};
177	Line Loop(109) = {61,-60,68};
178	Plane Surface(110) = {109};
179	Line Loop(111) = {60,-59,67};
180	Plane Surface(112) = {111};
181	Line Loop(113) = {58,-59,-66};
182	Plane Surface(114) = {113};
183	Line Loop(115) = {57,-58,-65};
184	Plane Surface(116) = {115};
185	Line Loop(117) = {72,57,-64};
186	Plane Surface(118) = {117};
187	Line Loop(119) = {63,-64,-71};
188	Plane Surface(120) = {119};
189
190
191	Physical Surface(155) = {106,110,112,114,116,118,120,108};
192	// END
193
194
195	// BEGIN
196	// Conduit Outer Surface
197	Line Loop(121) = {13,30,-21,-29};
198	Ruled Surface(122) = {121};
199	Line Loop(123) = {29,-20,-28,12};
200	Ruled Surface(124) = {123};
201	Line Loop(125) = {28,-19,-27,11};
202	Ruled Surface(126) = {125};
203	Line Loop(127) = {27,-18,-26,10};
204	Ruled Surface(128) = {127};
205	Line Loop(129) = {26,-17,-25,9};
206	Ruled Surface(130) = {129};
207	Line Loop(131) = {25,-24,-32,16};
208	Ruled Surface(132) = {131};
209	Line Loop(133) = {32,-23,-31,15};
210	Ruled Surface(134) = {133};
211	Line Loop(135) = {30,22,-31,-14};
212	Ruled Surface(136) = {135};
213
214
215	Physical Surface(156) = {134,136,122,126,128,130,132,124};
216	// END
217
218
219	// BEGIN
220	// Outer Free Surface Cylinder
221	Line Loop(137) = {53,38,-69,-37};
222	Ruled Surface(138) = {137};
223	Line Loop(139) = {52,37,-68,-36};
224	Ruled Surface(140) = {139};
225	Line Loop(141) = {51,36,-67,-35};
226	Ruled Surface(142) = {141};
227	Line Loop(143) = {50,35,-66,-34};
228	Ruled Surface(144) = {143};
229	Line Loop(145) = {49,34,-65,-33};
230	Ruled Surface(146) = {145};
231	Line Loop(147) = {56,33,-72,-40};
232	Ruled Surface(148) = {147};
233	Line Loop(149) = {55,40,-71,-39};
234	Ruled Surface(150) = {149};
235	Line Loop(151) = {39,-70,-38,54};
236	Ruled Surface(152) = {151};
237
238	Physical Surface(157) = {138,148,144,142,140,150,146,152};
239	// End
240
241
242	Surface Loop(158) =
243	{134,132,130,128,126,124,122,88,74,76,-78,-80,-82,-84,-86,-136,-90,152,150,92,94,-96,146,144,-98,-100,142,140,102,-104,138,108,106,-120,118,148,-116,-114,112,110};
244
245	Volume(159) = {158};
246	Physical Volume(160) = { 159};