/* boundary_phase.fe Evolver Code to model the shape of second phase which wet grain boundaries

Shapes of Grain Boundary Phases

The object is to calculate the shape and energy of a specified fixed volume inside the interstice of three spheres and contacts the three spheres with a specified wetting angle.
The arrangement of the three spheres is also specifiable. The opening angle (see figure below) can be specified as an angle between 0 (the sphere centers form an isosceles triangle) and 120 (the spheres are colinear). The gap between the three spheres can also be specified, a value of zero implies that the spheres are touching.
One particularly useful use of this code is the incremental calculation of the energy as a function of angle (a macro for this is provided at the bottom of the file: run_this) and then take numerical derivatives to get values for the torques that a meniscus applies to spheres. For more description of this file see this paper . The code comes with a rough approximation to the solution fo rthe default parameters. The mesh may not be robust for large incremental changes of the parameters--if it is not, you must resort to either changing the initial mesh for your problem, or use this one and make small incremental changes with intermittant relaxations towards equilibrium Three spheres with a meniscus. Wetting angle phi, opening angle theta, gap distance d and meniscus volume all variable parameters. Various macros at end for looping and recording results. Initial solution is rough approximation for phi=30,theta=0,d=0,volume=1.5 The sphere radii are all fixed at 1, and it should be straightforward to generalize. W. Craig Carter, NIST, wcraig@pruffle.nist.gov 11 march 1996 */ SYMMETRIC_CONTENT PARAMETER phi = @phi:30@ //Wetting angle (degrees) of liquid on solid PARAMETER theta = @theta:0@ //Opening Angle between bottom two spheres PARAMETER d = @d:0@//Gap between spheres (distances in units of 1 sphere radius) PARAMETER meniscus_volume = @meniscus_volume:1.5@//Volume of meniscus (in sphere radius^3) #define Rs 1 #define gamsl (-cos(phi*pi/180)) // virtual tension of facet on plane #define twoonrt3 1.1547005 #define circcent twoonrt3*(Rs + 0.5*d) #define radto2 theta*pi/360 #define ct2 cos( 60*pi/180 - theta*pi/360 ) #define st2 sin( 60*pi/180 - theta*pi/360 ) #define ct3 cos( 120*pi/180 + theta*pi/360 ) #define st3 sin( 120*pi/180 + theta*pi/360 ) #define r1sq (x)^2 + (y)^2 + z^2 #define r2sq (x - 2*Rs*ct2)^2 + (y - 2*Rs*st2)^2 + z^2 #define r3sq (x - 2*Rs*ct3)^2 + (y - 2*Rs*st3)^2 + z^2 #define w1 y/((x^2 + z^2)*sqrt(r1sq)) #define w2 (x*ct2 + y*st2 - 2*Rs)/(( (st2*x - ct2*y)^2 + z^2)*sqrt(r2sq)) #define w3 (x*ct3 + y*st3 - 2*Rs)/(( (st3*x - ct3*y)^2 + z^2)*sqrt(r3sq)) constraint 11 //bottom sphere formula: r1sq = Rs^2 constraint 1 //bottom sphere formula: r1sq = Rs^2 energy: e1: gamsl*(Rs^2)*w1*z e2: 0 e3: -gamsl*(Rs^2)*w1*x content: c1: ((Rs^3)/3)*w1*z c2: 0 c3: -((Rs^3)/3)*w1*x constraint 12 //northeast sphere theta = Pi/6 formula: r2sq = Rs^2 constraint 2 //northeast sphere theta = Pi/6 formula: r2sq = Rs^2 energy: e1: (gamsl*(Rs^2)*w2*(z*st2)) e2: -(gamsl*(Rs^2)*w2*(z*ct2)) e3: -(gamsl*(Rs^2)*w2*(x*st2 - y*ct2)) content: c1: -(-(((Rs^3)/3)*w2*(z*st2)) + z*(2*Rs)*st2/3) c2: -((((Rs^3)/3)*w2*(z*ct2)) - z*(2*Rs)*ct2/3) c3: -(((Rs^3)/3)*w2*(x*st2 - y*ct2)) constraint 13 formula: r3sq = Rs^2 constraint 3 formula: r3sq = Rs^2 energy: e1: (gamsl*(Rs^2)*w3*(z*st3)) e2: -(gamsl*(Rs^2)*w3*(z*ct3)) e3: -(gamsl*(Rs^2)*w3*(x*st3 - y*ct3)) content: c1: -(-(((Rs^3)/3)*w3*(z*st3)) + z*(2*Rs)*st3/3) c2: -((((Rs^3)/3)*w3*(z*ct3)) - z*(2*Rs)*ct3/3) c3: -(((Rs^3)/3)*w3*(x*st3 - y*ct3)) vertices /* coordinates */ 1 0.970462745604587 0.894248623238767 0.644330709717541 constraints 2 2 0.367518060798704 1.9845492403195 0.521808254978649 constraints 2 3 -0.367933385005334 1.98268445778052 0.523659384811855 constraints 3 4 -0.972181011159683 0.89254070544039 0.642823322569705 constraints 3 5 -0.663709681814714 0.407751757258772 0.627078912671918 constraints 1 6 0.665935461731492 0.40670416652271 0.625397219166772 constraints 1 7 0.972181011159683 0.89254070544039 -0.642823322569705 constraints 2 8 0.367933385005335 1.98268445778052 -0.523659384811855 constraints 2 9 -0.367518060798703 1.9845492403195 -0.521808254978649 constraints 3 10 -0.970462745604587 0.894248623238767 -0.644330709717541 constraints 3 11 -0.665935461731492 0.40670416652271 -0.625397219166772 constraints 1 12 0.663709681814714 0.407751757258771 -0.627078912671918 constraints 1 13 0.551333643089798 1.43718265406543 0.7775278864477 constraints 2 original 31 14 -0.552924595823299 1.43407021947537 0.77793398284754 constraints 3 original 33 15 0.00189114036368738 0.589839204141077 0.807518505575164 constraints 1 original 35 16 1.25276978667719 0.643897529886621 0.000814563675338598 constraints 2 original 37 17 0.329720933619116 2.21432852437754 -0.000935055537497777 constraints 2 original 38 18 -0.329720933619115 2.21432852437754 0.000935055537497808 constraints 3 original 39 19 -1.25276978667719 0.643897529886622 -0.000814563675339058 constraints 3 original 40 20 -0.970551771243615 0.240891177646197 0.000836581204673146 constraints 1 original 41 21 0.970551771243615 0.240891177646196 -0.000836581204673082 constraints 1 original 42 22 0.5529245958233 1.43407021947537 -0.77793398284754 constraints 2 original 43 23 -0.551333643089798 1.43718265406543 -0.7775278864477 constraints 3 original 45 24 -0.00189114036368764 0.589839204141077 -0.807518505575164 constraints 1 original 47 25 0.4322368803536 1.73322345608047 0.692737639772302 constraints 2 original 97 26 -0.742152459753616 1.14231697409822 0.768095988650899 constraints 3 original 98 27 0.359238563284428 0.53205517692794 0.766723511673874 constraints 1 original 99 28 1.16998799236389 0.712585783501712 0.377523743566876 constraints 2 original 100 29 0.33229014101403 2.14656235204332 -0.278697445494338 constraints 2 original 101 30 -0.332290141014029 2.14656235204332 0.278697445494338 constraints 3 original 102 31 -1.16998799236389 0.712585783501713 -0.377523743566876 constraints 3 original 103 32 -0.884954573638287 0.290817645905442 0.3637038622653 constraints 1 original 104 33 0.884954573638287 0.290817645905442 -0.3637038622653 constraints 1 original 105 34 0.742152459753616 1.14231697409822 -0.768095988650899 constraints 2 original 106 35 -0.432236880353599 1.73322345608047 -0.692737639772301 constraints 3 original 107 36 -0.359238563284429 0.532055176927941 -0.766723511673874 constraints 1 original 108 37 0.740071772414413 1.14495134522009 0.768692427704224 constraints 2 original 121 38 -0.433122348150238 1.7303824079405 0.694033135456501 constraints 3 original 125 39 -0.355914710077555 0.533082194980378 0.767559829945075 constraints 1 original 129 40 1.17077356580683 0.711947419159053 -0.375906083856486 constraints 2 original 133 41 0.332190027774218 2.14738145653237 0.276903945121237 constraints 2 original 135 42 -0.332190027774217 2.14738145653237 -0.276903945121238 constraints 3 original 137 43 -1.17077356580683 0.711947419159054 0.375906083856486 constraints 3 original 139 44 -0.885819161721006 0.290300115169157 -0.362008640588411 constraints 1 original 141 45 0.885819161721006 0.290300115169157 0.362008640588411 constraints 1 original 143 46 0.433122348150239 1.7303824079405 -0.694033135456501 constraints 2 original 145 47 -0.740071772414412 1.14495134522009 -0.768692427704224 constraints 3 original 149 48 0.355914710077555 0.533082194980378 -0.767559829945075 constraints 1 original 153 49 0.122171075180096 1.21448459334001 0.762494260589174 original 21 50 -0.122319095804889 1.21406058095504 0.762538040336152 original 22 51 0.000310161468331421 1.02138679581223 0.777230891307088 original 23 52 0.12231909580489 1.21406058095504 -0.762538040336152 original 31 53 -0.122171075180095 1.21448459334001 -0.762494260589174 original 32 54 -0.000310161468331302 1.02138679581223 -0.777230891307088 original 33 55 3.97617930747252e-16 2.14527915020098 -6.23550486379397e-17 original 19 56 -1.08329551380465 0.467683527519118 -1.17322854940636e-05 original 20 57 1.08329551380465 0.467683527519118 1.17322854940129e-05 original 21 58 -0.280320957897332 0.995642123220976 0.772319458790338 original 25 59 0.28074531219726 0.996190005146813 0.772298057290456 original 26 60 -0.000297038807248445 1.4451811310926 0.724512592331505 original 27 61 -0.28074531219726 0.996190005146813 -0.772298057290456 original 28 62 0.280320957897332 0.995642123220975 -0.772319458790339 original 29 63 0.000297038807248893 1.4451811310926 -0.724512592331506 original 30 64 -0.000295838759637342 1.86944411902946 0.536390639041099 original 32 65 -0.719014535407585 0.719543793196182 0.663376987911934 original 34 66 0.719397901999788 0.720134024243546 0.663386676067657 original 36 67 0.000295838759638036 1.86944411902946 -0.536390639041099 original 44 68 -0.719397901999788 0.720134024243546 -0.663386676067656 original 46 69 0.719014535407585 0.719543793196181 -0.663376987911935 original 48 70 0.524357676077355 1.07724550328381 0.76185711713097 original 49 71 0.22344683313521 1.59647927977215 0.692920719491509 original 50 72 -0.224232697054406 1.59505828197076 0.693473415403828 original 51 73 -0.525194576337507 1.07564449288542 0.761671988420135 original 52 74 -0.320274301988748 0.750026309752213 0.76819539994323 original 53 75 0.322155534337539 0.749904136760702 0.767913990540942 original 54 76 0.525194576337507 1.07564449288542 -0.761671988420135 original 55 77 0.224232697054407 1.59505828197076 -0.693473415403829 original 56 78 -0.22344683313521 1.59647927977215 -0.692920719491509 original 57 79 -0.524357676077354 1.07724550328382 -0.76185711713097 original 58 80 -0.322155534337539 0.749904136760702 -0.767913990540942 original 59 81 0.320274301988747 0.750026309752214 -0.76819539994323 original 60 82 -0.678948279639799 0.943890691503712 0.727099565310183 original 85 83 -0.505385571783192 0.663711561397643 0.724676768307211 original 86 84 0.507086695010654 0.663448118537346 0.724126051196973 original 87 85 0.678138758370102 0.945299398915565 0.727564618562068 original 88 86 0.201688104892039 1.76217145666071 0.622110784064762 original 89 87 -0.202342111693436 1.7609437871644 0.622880293460434 original 90 88 -0.507086695010655 0.663448118537346 -0.724126051196973 original 91 89 -0.678138758370102 0.945299398915566 -0.727564618562068 original 92 90 0.678948279639798 0.943890691503712 -0.727099565310183 original 93 91 0.505385571783193 0.663711561397643 -0.724676768307211 original 94 92 -0.201688104892038 1.76217145666071 -0.622110784064762 original 95 93 0.202342111693437 1.7609437871644 -0.622880293460433 original 96 94 5.43737543111337e-05 2.10054864914211 0.213154535041585 original 123 95 -0.00035143157835739 1.67077335665272 0.65187603351891 original 124 96 -0.994242821703562 0.529628380214522 0.36857858365851 original 127 97 -0.504890896224631 0.858504296316651 0.74453436540231 original 128 98 0.994298762635207 0.529753027556249 0.36860812071957 original 131 99 0.505321649374088 0.859155827055092 0.744519061077177 original 132 100 -5.43737543103433e-05 2.10054864914211 -0.213154535041585 original 147 101 0.000351431578358025 1.67077335665272 -0.65187603351891 original 148 102 -0.994298762635207 0.52975302755625 -0.36860812071957 original 151 103 -0.505321649374089 0.859155827055092 -0.744519061077176 original 152 104 0.994242821703563 0.529628380214521 -0.36857858365851 original 155 105 0.504890896224631 0.858504296316651 -0.744534365402309 original 156 106 0.337408111249011 1.16388924646497 0.768222203929928 original 69 107 -0.190719135099347 1.41019732751936 0.737785073685354 original 63 108 -0.152361939532519 0.871734651818738 0.781935828539253 original 61 109 0.190719135099347 1.41019732751936 -0.737785073685353 original 67 110 -0.337408111249011 1.16388924646497 -0.768222203929927 original 65 111 0.152361939532519 0.871734651818738 -0.781935828539252 original 71 112 0.334366415870822 1.3255953994641 0.757250681188982 original 72 113 -0.33519123425147 1.32394349711367 0.757447925670839 original 62 114 0.00102960182361643 0.803305456721938 0.788649619574357 original 64 115 0.33519123425147 1.32394349711367 -0.757447925670839 original 232 116 -0.334366415870821 1.32559539946411 -0.757250681188983 original 233 117 -0.0010296018236165 0.803305456721938 -0.788649619574357 original 234 118 0.190006066117117 1.41128936479615 0.737551061947448 original 235 119 -0.337904623537835 1.16259722184101 0.76828240969906 original 236 120 0.153838374794352 0.871832523946008 0.781864923582665 original 237 121 0.337904623537835 1.16259722184101 -0.76828240969906 original 238 122 -0.190006066117117 1.41128936479615 -0.737551061947447 original 239 123 -0.153838374794352 0.871832523946008 -0.781864923582664 original 240 124 -0.115588200608966 2.01899393688147 0.392557510101857 original 241 125 0.115253308593132 2.01928323684085 0.392082327083669 original 242 126 0.115588200608967 2.01899393688147 -0.392557510101857 original 243 127 -0.115253308593131 2.01928323684085 -0.39208232708367 original 244 128 -0.838920490536508 0.524874629801739 0.528494393817325 original 245 129 -0.942657958095673 0.680481670246711 0.532954799396314 original 246 130 -0.942499191682268 0.681086401925103 -0.533400609681691 original 247 131 -0.839532116841955 0.524939191225632 -0.528052102342574 original 248 132 0.942499191682268 0.681086401925103 0.533400609681691 original 249 133 0.839532116841955 0.524939191225632 0.528052102342575 original 250 134 0.838920490536508 0.524874629801739 -0.528494393817325 original 251 135 0.942657958095674 0.680481670246711 -0.532954799396314 original 252 136 -0.17281992598599 2.15882089918434 0.0979944517058819 original 325 137 0.172720666405476 2.15898577847545 0.0968092524653314 original 326 138 0.17281992598599 2.15882089918434 -0.0979944517058819 original 327 139 -0.172720666405475 2.15898577847545 -0.0968092524653314 original 328 140 -1.02105965158604 0.411938937377235 0.190431081318517 original 329 141 -1.11705726624741 0.54799741349895 0.192464026066464 original 330 142 -1.11698404966121 0.548156427992592 -0.192985409076999 original 331 143 -1.02122546798832 0.411943053563617 -0.189965889474037 original 332 144 1.11698404966121 0.548156427992592 0.192985409076999 original 333 145 1.02122546798832 0.411943053563616 0.189965889474037 original 334 146 1.02105965158604 0.411938937377234 -0.190431081318517 original 335 147 1.11705726624741 0.547997413498949 -0.192464026066464 original 336 148 0.395092551298509 1.86781373175128 0.617843504992005 constraints 2 149 -0.856298195766006 1.0106283747946 0.722098854195828 constraints 3 150 0.521233080877141 0.471986745789289 0.711016587147277 constraints 1 151 1.07973459847937 0.79257031766753 0.529328976197615 constraints 2 152 0.347749967471368 2.07834659695377 -0.409044723631817 constraints 2 153 -0.347749967471367 2.07834659695377 0.409044723631817 constraints 3 154 -1.07973459847937 0.79257031766753 -0.529328976197615 constraints 3 155 -0.78725193559462 0.344804217693993 0.511218584719841 constraints 1 156 0.78725193559462 0.344804217693992 -0.511218584719842 constraints 1 157 0.856298195766006 1.0106283747946 -0.722098854195828 constraints 2 158 -0.395092551298508 1.86781373175128 -0.617843504992005 constraints 3 159 -0.521233080877141 0.471986745789289 -0.711016587147276 constraints 1 160 0.637783103230827 1.28811875225052 0.786074212889604 constraints 2 161 -0.484725073325357 1.58452557136577 0.747192934722587 constraints 3 162 -0.180657103484788 0.574964416858922 0.797984292017435 constraints 1 163 1.23399617093618 0.661559028860935 -0.19602890059929 constraints 2 164 0.330523191825805 2.19825201145671 0.139908141205945 constraints 2 165 -0.330523191825804 2.19825201145671 -0.139908141205946 constraints 3 166 -1.23399617093618 0.661559028860936 0.196028900599289 constraints 3 167 -0.949373082502654 0.251622148488379 -0.188087864067572 constraints 1 168 0.949373082502654 0.251622148488379 0.188087864067572 constraints 1 169 0.484725073325357 1.58452557136577 -0.747192934722587 constraints 2 170 -0.637783103230826 1.28811875225052 -0.786074212889605 constraints 3 171 0.180657103484788 0.574964416858922 -0.797984292017435 constraints 1 172 0.48348554110535 1.58759847504925 0.746307318554452 constraints 2 173 -0.639680963640519 1.28516032378393 0.785978207649859 constraints 3 174 0.184317717571297 0.57436211206717 0.797580806696756 constraints 1 175 1.23362945346346 0.661857177875963 0.197669500834862 constraints 2 176 0.33054560602483 2.19782669711323 -0.141728693452689 constraints 2 177 -0.330545606024829 2.19782669711323 0.141728693452689 constraints 3 178 -1.23362945346346 0.661857177875964 -0.197669500834863 constraints 3 179 -0.948973582252083 0.251850399885149 0.189790717015715 constraints 1 180 0.948973582252084 0.251850399885149 -0.189790717015715 constraints 1 181 0.639680963640519 1.28516032378393 -0.785978207649859 constraints 2 182 -0.483485541105349 1.58759847504925 -0.746307318554452 constraints 3 183 -0.184317717571297 0.57436211206717 -0.797580806696756 constraints 1 184 0.854261807885837 1.01283615324939 0.723190629105899 constraints 2 185 -0.395701400603284 1.86539032406841 0.619477624076387 constraints 3 186 -0.518396325992008 0.473157601558487 0.712311121129951 constraints 1 187 1.0809193770334 0.791437314268184 -0.527664467840711 constraints 2 188 0.347507580847872 2.0796272345402 0.407187043081423 constraints 2 189 -0.347507580847871 2.0796272345402 -0.407187043081423 constraints 3 190 -1.0809193770334 0.791437314268184 0.52766446784071 constraints 3 191 -0.788731658196891 0.344106749884629 -0.509404471949178 constraints 1 192 0.788731658196891 0.344106749884629 0.509404471949178 constraints 1 193 0.395701400603284 1.86539032406841 -0.619477624076387 constraints 2 194 -0.854261807885836 1.01283615324939 -0.723190629105898 constraints 3 195 0.518396325992008 0.473157601558487 -0.712311121129951 constraints 1 196 -0.829903849651348 0.834242712402469 0.660156748960846 197 -0.658967442618513 0.569317067120127 0.654319294931162 198 0.660158087462342 0.569123147261463 0.653638890529105 199 0.829351313836781 0.835273372487837 0.660767466890109 200 0.196264067963671 1.90233296942796 0.53105955071766 201 -0.196802223490217 1.90146037281903 0.531909437120419 202 -0.660158087462343 0.569123147261464 -0.653638890529105 203 -0.82935131383678 0.835273372487837 -0.660767466890109 204 0.829903849651348 0.834242712402469 -0.660156748960846 205 0.658967442618513 0.569317067120126 -0.654319294931162 206 -0.19626406796367 1.90233296942796 -0.53105955071766 207 0.196802223490218 1.90146037281903 -0.531909437120418 208 0.476271884328847 1.2570937739244 0.771926085629617 209 0.352152662998479 1.47611574608606 0.742200193736253 210 -0.353181445388453 1.47408972978778 0.742698098075864 211 -0.477372782813292 1.25498532343639 0.77196851670094 212 -0.132407961122782 0.716122900926716 0.790780100996378 213 0.134976080133511 0.715988796555255 0.790604299369227 214 0.477372782813292 1.25498532343639 -0.77196851670094 215 0.353181445388454 1.47408972978778 -0.742698098075864 216 -0.352152662998478 1.47611574608606 -0.742200193736253 217 -0.476271884328847 1.2570937739244 -0.771926085629616 218 -0.134976080133511 0.715988796555256 -0.790604299369227 219 0.132407961122782 0.716122900926716 -0.790780100996378 edges 1 2 148 constraints 2 2 4 149 constraints 3 original 3 3 6 150 constraints 1 original 5 4 1 151 constraints 2 original 7 5 8 152 constraints 2 original 8 6 3 153 constraints 3 original 9 7 10 154 constraints 3 original 10 8 5 155 constraints 1 original 11 9 12 156 constraints 1 original 12 10 7 157 constraints 2 original 13 11 9 158 constraints 3 original 15 12 11 159 constraints 1 original 17 13 13 160 constraints 2 original 85 14 14 161 constraints 3 original 91 15 15 162 constraints 1 original 97 16 16 163 constraints 2 original 103 17 17 164 constraints 2 original 106 18 18 165 constraints 3 original 109 19 19 166 constraints 3 original 112 20 20 167 constraints 1 original 115 21 21 168 constraints 1 original 118 22 22 169 constraints 2 original 121 23 23 170 constraints 3 original 127 24 24 171 constraints 1 original 133 25 25 172 constraints 2 original 283 26 26 173 constraints 3 original 286 27 27 174 constraints 1 original 289 28 28 175 constraints 2 original 292 29 29 176 constraints 2 original 295 30 30 177 constraints 3 original 298 31 31 178 constraints 3 original 301 32 32 179 constraints 1 original 304 33 33 180 constraints 1 original 307 34 34 181 constraints 2 original 310 35 35 182 constraints 3 original 313 36 36 183 constraints 1 original 316 37 37 184 constraints 2 original 355 38 38 185 constraints 3 original 367 39 39 186 constraints 1 original 379 40 40 187 constraints 2 original 391 41 41 188 constraints 2 original 397 42 42 189 constraints 3 original 403 43 43 190 constraints 3 original 409 44 44 191 constraints 1 original 415 45 45 192 constraints 1 original 421 46 46 193 constraints 2 original 427 47 47 194 constraints 3 original 439 48 48 195 constraints 1 original 451 49 201 124 original 2 50 197 128 original 4 51 199 132 original 6 52 207 126 original 14 53 203 130 original 16 54 205 134 original 18 55 49 106 original 21 56 86 25 original 22 57 50 107 original 23 58 82 26 original 24 59 51 108 original 25 60 84 27 original 26 61 49 50 original 27 62 50 51 original 28 63 51 49 original 29 64 90 34 original 31 65 52 109 original 32 66 92 35 original 33 67 53 110 original 34 68 88 36 original 35 69 54 111 original 36 70 53 52 original 37 71 54 53 original 38 72 52 54 original 39 73 4 196 original 61 74 5 197 original 62 75 51 58 original 63 76 50 58 original 64 77 6 198 original 65 78 1 199 original 66 79 49 59 original 67 80 51 59 original 68 81 2 200 original 69 82 3 201 original 70 83 50 60 original 71 84 49 60 original 72 85 11 202 original 73 86 10 203 original 74 87 53 61 original 75 88 54 61 original 76 89 7 204 original 77 90 12 205 original 78 91 54 62 original 79 92 52 62 original 80 93 9 206 original 81 94 8 207 original 82 95 52 63 original 83 96 53 63 original 84 97 13 112 original 87 98 200 125 original 2 99 125 124 original 89 100 64 95 original 90 101 14 113 original 93 102 196 129 original 4 103 129 128 original 95 104 65 97 original 96 105 15 114 original 99 106 198 133 original 6 107 133 132 original 101 108 66 99 original 102 109 16 57 original 104 110 138 137 original 107 111 136 139 original 110 112 19 56 original 113 113 20 56 original 116 114 21 57 original 119 115 22 115 original 123 116 206 127 original 14 117 127 126 original 125 118 67 101 original 126 119 23 116 original 129 120 202 131 original 16 121 131 130 original 131 122 68 103 original 132 123 24 117 original 135 124 204 135 original 18 125 135 134 original 137 126 69 105 original 138 127 37 85 original 21 128 70 208 original 140 129 70 59 original 141 130 71 118 original 22 131 71 209 original 143 132 71 60 original 144 133 38 87 original 23 134 72 210 original 146 135 72 60 original 147 136 73 119 original 24 137 73 211 original 149 138 73 58 original 150 139 39 83 original 25 140 74 212 original 152 141 74 58 original 153 142 75 120 original 26 143 75 213 original 155 144 75 59 original 156 145 76 121 original 31 146 76 214 original 158 147 76 62 original 159 148 46 93 original 32 149 77 215 original 161 150 77 63 original 162 151 78 122 original 33 152 78 216 original 164 153 78 63 original 165 154 47 89 original 34 155 79 217 original 167 156 79 61 original 168 157 80 123 original 35 158 80 218 original 170 159 80 61 original 171 160 48 91 original 36 161 81 219 original 173 162 81 62 original 174 163 18 136 original 175 164 30 124 original 177 165 94 137 original 178 166 41 125 original 179 167 17 138 original 181 168 29 126 original 183 169 100 139 original 184 170 42 127 original 185 171 20 140 original 187 172 32 128 original 189 173 96 141 original 190 174 43 129 original 191 175 19 142 original 193 176 31 130 original 195 177 102 143 original 196 178 44 131 original 197 179 16 144 original 199 180 28 132 original 201 181 98 145 original 202 182 45 133 original 203 183 21 146 original 205 184 33 134 original 207 185 104 147 original 208 186 40 135 original 209 187 73 97 original 247 188 82 73 original 248 189 82 65 original 249 190 97 74 original 250 191 83 65 original 251 192 83 74 original 252 193 75 99 original 253 194 84 75 original 254 195 84 66 original 255 196 99 70 original 256 197 85 66 original 257 198 85 70 original 258 199 71 95 original 259 200 86 71 original 260 201 86 64 original 261 202 95 72 original 262 203 87 64 original 263 204 87 72 original 264 205 80 103 original 265 206 88 80 original 266 207 88 68 original 267 208 103 79 original 268 209 89 68 original 269 210 89 79 original 270 211 76 105 original 271 212 90 76 original 272 213 90 69 original 273 214 105 81 original 274 215 91 69 original 275 216 91 81 original 276 217 78 101 original 277 218 92 78 original 278 219 92 67 original 279 220 101 77 original 280 221 93 67 original 281 222 93 77 original 282 223 25 71 original 285 224 26 73 original 288 225 27 75 original 291 226 34 76 original 312 227 35 78 original 315 228 36 80 original 318 229 37 70 original 356 230 94 55 original 361 231 94 30 original 362 232 94 41 original 363 233 95 60 original 364 234 95 87 original 365 235 95 86 original 366 236 38 72 original 368 237 140 141 original 373 238 96 32 original 374 239 96 43 original 375 240 97 58 original 376 241 97 83 original 377 242 97 82 original 378 243 39 74 original 380 244 144 145 original 385 245 98 28 original 386 246 98 45 original 387 247 99 59 original 388 248 99 85 original 389 249 99 84 original 390 250 46 77 original 428 251 100 55 original 433 252 100 29 original 434 253 100 42 original 435 254 101 63 original 436 255 101 93 original 437 256 101 92 original 438 257 47 79 original 440 258 142 143 original 445 259 102 31 original 446 260 102 44 original 447 261 103 61 original 448 262 103 89 original 449 263 103 88 original 450 264 48 81 original 452 265 146 147 original 457 266 104 33 original 458 267 104 40 original 459 268 105 62 original 460 269 105 91 original 461 270 105 90 original 462 271 106 70 original 21 272 106 208 original 200 273 106 59 original 293 274 107 72 original 23 275 107 210 original 182 276 107 60 original 296 277 108 74 original 25 278 108 212 original 176 279 108 58 original 299 280 109 77 original 32 281 109 215 original 194 282 109 63 original 302 283 110 79 original 34 284 110 217 original 188 285 110 61 original 305 286 111 81 original 36 287 111 219 original 206 288 111 62 original 308 289 112 49 original 210 290 118 209 original 48 291 112 106 original 393 292 113 50 original 180 293 119 211 original 38 294 113 107 original 399 295 114 51 original 186 296 120 213 original 40 297 114 108 original 405 298 115 52 original 688 299 121 214 original 689 300 115 109 original 690 301 116 53 original 691 302 122 216 original 692 303 116 110 original 693 304 117 54 original 694 305 123 218 original 695 306 117 111 original 696 307 118 49 original 22 308 118 60 original 698 309 118 112 original 699 310 119 50 original 24 311 119 58 original 701 312 119 113 original 702 313 120 51 original 26 314 120 59 original 704 315 120 114 original 705 316 121 52 original 31 317 121 62 original 707 318 121 115 original 708 319 122 53 original 33 320 122 63 original 710 321 122 116 original 711 322 123 54 original 35 323 123 61 original 713 324 123 117 original 714 325 124 64 original 715 326 124 3 original 716 327 124 94 original 717 328 125 64 original 718 329 125 2 original 719 330 125 94 original 720 331 126 67 original 721 332 126 8 original 722 333 126 100 original 723 334 127 67 original 724 335 127 9 original 725 336 127 100 original 726 337 128 65 original 727 338 128 5 original 728 339 128 96 original 729 340 129 65 original 730 341 129 4 original 731 342 129 96 original 732 343 130 68 original 733 344 130 10 original 734 345 130 102 original 735 346 131 68 original 736 347 131 11 original 737 348 131 102 original 738 349 132 66 original 739 350 132 1 original 740 351 132 98 original 741 352 133 66 original 742 353 133 6 original 743 354 133 98 original 744 355 134 69 original 745 356 134 12 original 746 357 134 104 original 747 358 135 69 original 748 359 135 7 original 749 360 135 104 original 750 361 136 94 original 967 362 136 30 original 968 363 136 55 original 969 364 137 17 original 970 365 137 41 original 971 366 137 55 original 972 367 138 100 original 973 368 138 29 original 974 369 138 55 original 975 370 139 18 original 976 371 139 42 original 977 372 139 55 original 978 373 140 96 original 979 374 140 32 original 980 375 140 56 original 981 376 141 19 original 982 377 141 43 original 983 378 141 56 original 984 379 142 102 original 985 380 142 31 original 986 381 142 56 original 987 382 143 20 original 988 383 143 44 original 989 384 143 56 original 990 385 144 98 original 991 386 144 28 original 992 387 144 57 original 993 388 145 21 original 994 389 145 45 original 995 390 145 57 original 996 391 146 104 original 997 392 146 33 original 998 393 146 57 original 999 394 147 16 original 1000 395 147 40 original 1001 396 147 57 original 1002 397 148 25 constraints 2 original 1 398 148 86 399 149 26 constraints 3 original 3 400 149 82 401 150 27 constraints 1 original 5 402 150 84 403 151 28 constraints 2 original 7 404 151 132 405 152 29 constraints 2 original 8 406 152 126 407 153 30 constraints 3 original 9 408 153 124 409 154 31 constraints 3 original 10 410 154 130 411 155 32 constraints 1 original 11 412 155 128 413 156 33 constraints 1 original 12 414 156 134 415 157 34 constraints 2 original 13 416 157 90 417 158 35 constraints 3 original 15 418 158 92 419 159 36 constraints 1 original 17 420 159 88 421 160 37 constraints 2 original 85 422 160 70 423 161 38 constraints 3 original 91 424 161 72 425 162 39 constraints 1 original 97 426 162 74 427 163 40 constraints 2 original 103 428 163 147 429 164 41 constraints 2 original 106 430 164 137 431 165 42 constraints 3 original 109 432 165 139 433 166 43 constraints 3 original 112 434 166 141 435 167 44 constraints 1 original 115 436 167 143 437 168 45 constraints 1 original 118 438 168 145 439 169 46 constraints 2 original 121 440 169 77 441 170 47 constraints 3 original 127 442 170 79 443 171 48 constraints 1 original 133 444 171 81 445 172 13 constraints 2 original 283 446 172 71 447 173 14 constraints 3 original 286 448 173 73 449 174 15 constraints 1 original 289 450 174 75 451 175 16 constraints 2 original 292 452 175 144 453 176 17 constraints 2 original 295 454 176 138 455 177 18 constraints 3 original 298 456 177 136 457 178 19 constraints 3 original 301 458 178 142 459 179 20 constraints 1 original 304 460 179 140 461 180 21 constraints 1 original 307 462 180 146 463 181 22 constraints 2 original 310 464 181 76 465 182 23 constraints 3 original 313 466 182 78 467 183 24 constraints 1 original 316 468 183 80 469 184 1 constraints 2 original 355 470 184 85 471 185 3 constraints 3 original 367 472 185 87 473 186 5 constraints 1 original 379 474 186 83 475 187 7 constraints 2 original 391 476 187 135 477 188 2 constraints 2 original 397 478 188 125 479 189 9 constraints 3 original 403 480 189 127 481 190 4 constraints 3 original 409 482 190 129 483 191 11 constraints 1 original 415 484 191 131 485 192 6 constraints 1 original 421 486 192 133 487 193 8 constraints 2 original 427 488 193 93 489 194 10 constraints 3 original 439 490 194 89 491 195 12 constraints 1 original 451 492 195 91 493 196 82 original 61 494 196 65 495 196 149 496 197 83 original 62 497 197 186 498 197 65 499 198 84 original 65 500 198 66 501 198 150 502 199 85 original 66 503 199 184 504 199 66 505 200 86 original 69 506 200 64 507 200 148 508 201 87 original 70 509 201 185 510 201 64 511 202 88 original 73 512 202 68 513 202 159 514 203 89 original 74 515 203 194 516 203 68 517 204 90 original 77 518 204 69 519 204 157 520 205 91 original 78 521 205 195 522 205 69 523 206 92 original 81 524 206 67 525 206 158 526 207 93 original 82 527 207 193 528 207 67 529 208 13 original 140 530 208 112 531 208 160 532 209 13 original 143 533 209 172 534 209 112 535 210 14 original 146 536 210 113 537 210 161 538 211 14 original 149 539 211 173 540 211 113 541 212 15 original 152 542 212 114 543 212 162 544 213 15 original 155 545 213 174 546 213 114 547 214 22 original 158 548 214 181 549 214 115 550 215 22 original 161 551 215 115 552 215 169 553 216 23 original 164 554 216 182 555 216 116 556 217 23 original 167 557 217 116 558 217 170 559 218 24 original 170 560 218 183 561 218 117 562 219 24 original 173 563 219 117 564 219 171 faces /* edge loop */ 1 432 371 -431 tag 1 COLOR 2 2 436 383 -435 tag 1 COLOR 2 3 428 395 -427 tag 1 COLOR 2 4 503 470 -502 tag 3 COLOR 2 original 11 5 509 472 -508 tag 3 COLOR 2 original 12 6 497 474 -496 tag 3 COLOR 2 original 13 7 61 62 63 tag 3 COLOR 2 original 14 8 311 -76 -310 tag 4 COLOR 2 original 15 9 314 -80 -313 tag 4 COLOR 2 original 16 10 308 -84 -307 tag 4 COLOR 2 original 17 11 527 488 -526 tag 3 COLOR 2 original 21 12 515 490 -514 tag 3 COLOR 2 original 22 13 521 492 -520 tag 3 COLOR 2 original 23 14 72 71 70 tag 3 COLOR 2 original 24 15 323 -88 -322 tag 4 COLOR 2 original 25 16 317 -92 -316 tag 4 COLOR 2 original 26 17 320 -96 -319 tag 4 COLOR 2 original 27 18 110 364 167 tag 1 COLOR 2 original 1 19 430 365 -429 tag 1 COLOR 2 original 1 20 111 370 163 tag 1 COLOR 2 original 1 21 -378 376 112 tag 1 COLOR 2 original 2 22 434 377 -433 tag 1 COLOR 2 original 2 23 -384 382 113 tag 1 COLOR 2 original 2 24 -390 388 114 tag 1 COLOR 2 original 3 25 438 389 -437 tag 1 COLOR 2 original 3 26 -396 394 109 tag 1 COLOR 2 original 3 27 -187 -188 -242 tag 4 COLOR 2 original 15 28 -190 241 192 tag 4 COLOR 2 original 15 29 76 -75 -62 tag 4 COLOR 2 original 15 30 -193 -194 -249 tag 4 COLOR 2 original 16 31 -196 248 198 tag 4 COLOR 2 original 16 32 80 -79 -63 tag 4 COLOR 2 original 16 33 -199 -200 -235 tag 4 COLOR 2 original 17 34 -202 234 204 tag 4 COLOR 2 original 17 35 84 -83 -61 tag 4 COLOR 2 original 17 36 -205 -206 -263 tag 4 COLOR 2 original 25 37 -208 262 210 tag 4 COLOR 2 original 25 38 88 -87 -71 tag 4 COLOR 2 original 25 39 -211 -212 -270 tag 4 COLOR 2 original 26 40 -214 269 216 tag 4 COLOR 2 original 26 41 92 -91 -72 tag 4 COLOR 2 original 26 42 -217 -218 -256 tag 4 COLOR 2 original 27 43 -220 255 222 tag 4 COLOR 2 original 27 44 96 -95 -70 tag 4 COLOR 2 original 27 45 309 289 -307 tag 3 COLOR 2 original 11 46 362 -231 -361 tag 1 COLOR 2 original 1 47 -202 233 -135 tag 4 COLOR 2 original 17 48 312 292 -310 tag 3 COLOR 2 original 12 49 374 -238 -373 tag 1 COLOR 2 original 2 50 -190 240 -141 tag 4 COLOR 2 original 15 51 315 295 -313 tag 3 COLOR 2 original 13 52 386 -245 -385 tag 1 COLOR 2 original 3 53 -196 247 -129 tag 4 COLOR 2 original 16 54 -244 387 -390 tag 1 COLOR 2 original 3 55 -251 -367 369 tag 1 COLOR 2 original 1 56 -230 -361 363 tag 1 COLOR 2 original 1 57 -258 381 -384 tag 1 COLOR 2 original 2 58 -237 375 -378 tag 1 COLOR 2 original 2 59 -265 393 -396 tag 1 COLOR 2 original 3 60 318 298 -316 tag 3 COLOR 2 original 21 61 368 -252 -367 tag 1 COLOR 2 original 1 62 -220 254 -150 tag 4 COLOR 2 original 27 63 321 301 -319 tag 3 COLOR 2 original 22 64 380 -259 -379 tag 1 COLOR 2 original 2 65 -208 261 -156 tag 4 COLOR 2 original 25 66 324 304 -322 tag 3 COLOR 2 original 23 67 392 -266 -391 tag 1 COLOR 2 original 3 68 -214 268 -162 tag 4 COLOR 2 original 26 69 -530 529 97 tag 3 COLOR 2 original 11 70 -273 271 129 tag 4 COLOR 2 original 16 71 -533 532 -445 tag 3 COLOR 2 original 11 72 -199 132 -233 tag 4 COLOR 2 original 17 73 -536 535 101 tag 3 COLOR 2 original 12 74 -276 274 135 tag 4 COLOR 2 original 17 75 -539 538 -447 tag 3 COLOR 2 original 12 76 -187 138 -240 tag 4 COLOR 2 original 15 77 -542 541 105 tag 3 COLOR 2 original 13 78 -279 277 141 tag 4 COLOR 2 original 15 79 -545 544 -449 tag 3 COLOR 2 original 13 80 -193 144 -247 tag 4 COLOR 2 original 16 81 -548 547 -463 tag 3 COLOR 2 original 21 82 -211 147 -268 tag 4 COLOR 2 original 26 83 -551 550 115 tag 3 COLOR 2 original 21 84 -282 280 150 tag 4 COLOR 2 original 27 85 -554 553 -465 tag 3 COLOR 2 original 22 86 -217 153 -254 tag 4 COLOR 2 original 27 87 -557 556 119 tag 3 COLOR 2 original 22 88 -285 283 156 tag 4 COLOR 2 original 25 89 -560 559 -467 tag 3 COLOR 2 original 23 90 -205 159 -261 tag 4 COLOR 2 original 25 91 -563 562 123 tag 3 COLOR 2 original 23 92 -288 286 162 tag 4 COLOR 2 original 26 93 -49 -82 -326 tag 1 COLOR 2 original 1 94 98 328 -506 tag 1 COLOR 2 original 1 95 -52 -94 -332 tag 1 COLOR 2 original 1 96 116 334 -524 tag 1 COLOR 2 original 1 97 -50 -74 -338 tag 1 COLOR 2 original 2 98 102 340 -494 tag 1 COLOR 2 original 2 99 -53 -86 -344 tag 1 COLOR 2 original 2 100 120 346 -512 tag 1 COLOR 2 original 2 101 -51 -78 -350 tag 1 COLOR 2 original 3 102 106 352 -500 tag 1 COLOR 2 original 3 103 -54 -90 -356 tag 1 COLOR 2 original 3 104 124 358 -518 tag 1 COLOR 2 original 3 105 -58 188 -224 tag 4 COLOR 2 original 15 106 -494 493 189 tag 4 COLOR 2 original 15 107 -498 496 191 tag 4 COLOR 2 original 15 108 139 192 -243 tag 4 COLOR 2 original 15 109 -60 194 -225 tag 4 COLOR 2 original 16 110 -500 499 195 tag 4 COLOR 2 original 16 111 -504 502 197 tag 4 COLOR 2 original 16 112 127 198 -229 tag 4 COLOR 2 original 16 113 -56 200 -223 tag 4 COLOR 2 original 17 114 -506 505 201 tag 4 COLOR 2 original 17 115 -510 508 203 tag 4 COLOR 2 original 17 116 133 204 -236 tag 4 COLOR 2 original 17 117 -68 206 -228 tag 4 COLOR 2 original 25 118 -512 511 207 tag 4 COLOR 2 original 25 119 -516 514 209 tag 4 COLOR 2 original 25 120 154 210 -257 tag 4 COLOR 2 original 25 121 -64 212 -226 tag 4 COLOR 2 original 26 122 -518 517 213 tag 4 COLOR 2 original 26 123 -522 520 215 tag 4 COLOR 2 original 26 124 160 216 -264 tag 4 COLOR 2 original 26 125 -66 218 -227 tag 4 COLOR 2 original 27 126 -524 523 219 tag 4 COLOR 2 original 27 127 -528 526 221 tag 4 COLOR 2 original 27 128 148 222 -250 tag 4 COLOR 2 original 27 129 -398 397 -56 tag 3 COLOR 2 original 11 130 -400 399 -58 tag 3 COLOR 2 original 12 131 -402 401 -60 tag 3 COLOR 2 original 13 132 -416 415 -64 tag 3 COLOR 2 original 21 133 -418 417 -66 tag 3 COLOR 2 original 22 134 -420 419 -68 tag 3 COLOR 2 original 23 135 -422 421 229 tag 3 COLOR 2 original 11 136 99 325 -328 tag 1 COLOR 2 original 1 137 99 327 -330 tag 1 COLOR 2 original 1 138 234 203 100 tag 4 COLOR 2 original 17 139 235 201 100 tag 4 COLOR 2 original 17 140 -424 423 236 tag 3 COLOR 2 original 12 141 103 337 -340 tag 1 COLOR 2 original 2 142 103 339 -342 tag 1 COLOR 2 original 2 143 241 191 104 tag 4 COLOR 2 original 15 144 242 189 104 tag 4 COLOR 2 original 15 145 -426 425 243 tag 3 COLOR 2 original 13 146 107 349 -352 tag 1 COLOR 2 original 3 147 107 351 -354 tag 1 COLOR 2 original 3 148 248 197 108 tag 4 COLOR 2 original 16 149 249 195 108 tag 4 COLOR 2 original 16 150 -440 439 250 tag 3 COLOR 2 original 21 151 117 331 -334 tag 1 COLOR 2 original 1 152 117 333 -336 tag 1 COLOR 2 original 1 153 255 221 118 tag 4 COLOR 2 original 27 154 256 219 118 tag 4 COLOR 2 original 27 155 -442 441 257 tag 3 COLOR 2 original 22 156 121 343 -346 tag 1 COLOR 2 original 2 157 121 345 -348 tag 1 COLOR 2 original 2 158 262 209 122 tag 4 COLOR 2 original 25 159 263 207 122 tag 4 COLOR 2 original 25 160 -444 443 264 tag 3 COLOR 2 original 23 161 125 355 -358 tag 1 COLOR 2 original 3 162 125 357 -360 tag 1 COLOR 2 original 3 163 269 215 126 tag 4 COLOR 2 original 26 164 270 213 126 tag 4 COLOR 2 original 26 165 -291 289 55 tag 3 COLOR 2 original 11 166 273 -79 55 tag 4 COLOR 2 original 16 167 -294 292 57 tag 3 COLOR 2 original 12 168 276 -83 57 tag 4 COLOR 2 original 17 169 -297 295 59 tag 3 COLOR 2 original 13 170 279 -75 59 tag 4 COLOR 2 original 15 171 -300 298 65 tag 3 COLOR 2 original 21 172 282 -95 65 tag 4 COLOR 2 original 27 173 -303 301 67 tag 3 COLOR 2 original 22 174 285 -87 67 tag 4 COLOR 2 original 25 175 -306 304 69 tag 3 COLOR 2 original 23 176 288 -91 69 tag 4 COLOR 2 original 26 177 -534 532 97 tag 3 COLOR 2 original 11 178 272 530 291 tag 3 COLOR 2 original 11 179 -540 538 101 tag 3 COLOR 2 original 12 180 275 536 294 tag 3 COLOR 2 original 12 181 -546 544 105 tag 3 COLOR 2 original 13 182 278 542 297 tag 3 COLOR 2 original 13 183 -549 547 115 tag 3 COLOR 2 original 21 184 281 551 300 tag 3 COLOR 2 original 21 185 -555 553 119 tag 3 COLOR 2 original 22 186 284 557 303 tag 3 COLOR 2 original 22 187 -561 559 123 tag 3 COLOR 2 original 23 188 287 563 306 tag 3 COLOR 2 original 23 189 308 -132 130 tag 4 COLOR 2 original 17 190 290 -131 130 tag 3 COLOR 2 original 11 191 311 -138 136 tag 4 COLOR 2 original 15 192 293 -137 136 tag 3 COLOR 2 original 12 193 314 -144 142 tag 4 COLOR 2 original 16 194 296 -143 142 tag 3 COLOR 2 original 13 195 317 -147 145 tag 4 COLOR 2 original 26 196 299 -146 145 tag 3 COLOR 2 original 21 197 320 -153 151 tag 4 COLOR 2 original 27 198 302 -152 151 tag 3 COLOR 2 original 22 199 323 -159 157 tag 4 COLOR 2 original 25 200 305 -158 157 tag 3 COLOR 2 original 23 201 -408 407 164 tag 1 COLOR 2 original 1 202 327 231 164 tag 1 COLOR 2 original 1 203 478 329 -477 tag 1 COLOR 2 original 1 204 330 232 166 tag 1 COLOR 2 original 1 205 -406 405 168 tag 1 COLOR 2 original 1 206 333 252 168 tag 1 COLOR 2 original 1 207 480 335 -479 tag 1 COLOR 2 original 1 208 336 253 170 tag 1 COLOR 2 original 1 209 -412 411 172 tag 1 COLOR 2 original 2 210 339 238 172 tag 1 COLOR 2 original 2 211 482 341 -481 tag 1 COLOR 2 original 2 212 342 239 174 tag 1 COLOR 2 original 2 213 -410 409 176 tag 1 COLOR 2 original 2 214 345 259 176 tag 1 COLOR 2 original 2 215 484 347 -483 tag 1 COLOR 2 original 2 216 348 260 178 tag 1 COLOR 2 original 2 217 -404 403 180 tag 1 COLOR 2 original 3 218 351 245 180 tag 1 COLOR 2 original 3 219 486 353 -485 tag 1 COLOR 2 original 3 220 354 246 182 tag 1 COLOR 2 original 3 221 -414 413 184 tag 1 COLOR 2 original 3 222 357 266 184 tag 1 COLOR 2 original 3 223 476 359 -475 tag 1 COLOR 2 original 3 224 360 267 186 tag 1 COLOR 2 original 3 225 -456 455 163 tag 1 COLOR 2 original 1 226 -111 363 -372 tag 1 COLOR 2 original 1 227 365 -232 165 tag 1 COLOR 2 original 1 228 -230 165 366 tag 1 COLOR 2 original 1 229 -454 453 167 tag 1 COLOR 2 original 1 230 -110 369 -366 tag 1 COLOR 2 original 1 231 371 -253 169 tag 1 COLOR 2 original 1 232 -251 169 372 tag 1 COLOR 2 original 1 233 -460 459 171 tag 1 COLOR 2 original 2 234 375 -113 171 tag 1 COLOR 2 original 2 235 377 -239 173 tag 1 COLOR 2 original 2 236 -237 373 173 tag 1 COLOR 2 original 2 237 -458 457 175 tag 1 COLOR 2 original 2 238 381 -112 175 tag 1 COLOR 2 original 2 239 383 -260 177 tag 1 COLOR 2 original 2 240 -258 379 177 tag 1 COLOR 2 original 2 241 -452 451 179 tag 1 COLOR 2 original 3 242 387 -109 179 tag 1 COLOR 2 original 3 243 389 -246 181 tag 1 COLOR 2 original 3 244 -244 385 181 tag 1 COLOR 2 original 3 245 -462 461 183 tag 1 COLOR 2 original 3 246 393 -114 183 tag 1 COLOR 2 original 3 247 395 -267 185 tag 1 COLOR 2 original 3 248 -265 391 185 tag 1 COLOR 2 original 3 249 507 398 -505 tag 3 COLOR 2 original 11 250 495 400 -493 tag 3 COLOR 2 original 12 251 501 402 -499 tag 3 COLOR 2 original 13 252 404 350 4 tag 1 COLOR 2 original 3 253 406 332 5 tag 1 COLOR 2 original 1 254 408 326 6 tag 1 COLOR 2 original 1 255 410 344 7 tag 1 COLOR 2 original 2 256 412 338 8 tag 1 COLOR 2 original 2 257 414 356 9 tag 1 COLOR 2 original 3 258 519 416 -517 tag 3 COLOR 2 original 21 259 525 418 -523 tag 3 COLOR 2 original 22 260 513 420 -511 tag 3 COLOR 2 original 23 261 -531 529 13 tag 3 COLOR 2 original 11 262 -537 535 14 tag 3 COLOR 2 original 12 263 -543 541 15 tag 3 COLOR 2 original 13 264 428 394 16 tag 1 COLOR 2 original 3 265 430 364 17 tag 1 COLOR 2 original 1 266 432 370 18 tag 1 COLOR 2 original 1 267 434 376 19 tag 1 COLOR 2 original 2 268 436 382 20 tag 1 COLOR 2 original 2 269 438 388 21 tag 1 COLOR 2 original 3 270 -552 550 22 tag 3 COLOR 2 original 21 271 -558 556 23 tag 3 COLOR 2 original 22 272 -564 562 24 tag 3 COLOR 2 original 23 273 446 -223 25 tag 3 COLOR 2 original 11 274 448 -224 26 tag 3 COLOR 2 original 12 275 450 -225 27 tag 3 COLOR 2 original 13 276 452 386 28 tag 1 COLOR 2 original 3 277 454 368 29 tag 1 COLOR 2 original 1 278 456 362 30 tag 1 COLOR 2 original 1 279 458 380 31 tag 1 COLOR 2 original 2 280 460 374 32 tag 1 COLOR 2 original 2 281 462 392 33 tag 1 COLOR 2 original 3 282 464 -226 34 tag 3 COLOR 2 original 21 283 466 -227 35 tag 3 COLOR 2 original 22 284 468 -228 36 tag 3 COLOR 2 original 23 285 470 -127 37 tag 3 COLOR 2 original 11 286 472 -133 38 tag 3 COLOR 2 original 12 287 474 -139 39 tag 3 COLOR 2 original 13 288 476 -186 40 tag 1 COLOR 2 original 3 289 478 -166 41 tag 1 COLOR 2 original 1 290 480 -170 42 tag 1 COLOR 2 original 1 291 482 -174 43 tag 1 COLOR 2 original 2 292 484 -178 44 tag 1 COLOR 2 original 2 293 486 -182 45 tag 1 COLOR 2 original 3 294 488 -148 46 tag 3 COLOR 2 original 21 295 490 -154 47 tag 3 COLOR 2 original 22 296 492 -160 48 tag 3 COLOR 2 original 23 297 102 341 73 tag 4 COLOR 2 original 15 298 495 -2 73 tag 3 COLOR 2 original 12 299 497 473 74 tag 3 COLOR 2 original 13 300 -50 498 -337 tag 4 COLOR 2 original 15 301 106 353 77 tag 4 COLOR 2 original 16 302 501 -3 77 tag 3 COLOR 2 original 13 303 503 469 78 tag 3 COLOR 2 original 11 304 -51 504 -349 tag 4 COLOR 2 original 16 305 98 329 81 tag 4 COLOR 2 original 17 306 507 -1 81 tag 3 COLOR 2 original 11 307 509 471 82 tag 3 COLOR 2 original 12 308 -49 510 -325 tag 4 COLOR 2 original 17 309 120 347 85 tag 4 COLOR 2 original 25 310 513 -12 85 tag 3 COLOR 2 original 23 311 515 489 86 tag 3 COLOR 2 original 22 312 -53 516 -343 tag 4 COLOR 2 original 25 313 124 359 89 tag 4 COLOR 2 original 26 314 519 -10 89 tag 3 COLOR 2 original 21 315 521 491 90 tag 3 COLOR 2 original 23 316 -54 522 -355 tag 4 COLOR 2 original 26 317 116 335 93 tag 4 COLOR 2 original 27 318 525 -11 93 tag 3 COLOR 2 original 22 319 527 487 94 tag 3 COLOR 2 original 21 320 -52 528 -331 tag 4 COLOR 2 original 27 321 -272 271 128 tag 3 COLOR 2 original 11 322 531 422 128 tag 3 COLOR 2 original 11 323 533 446 131 tag 3 COLOR 2 original 11 324 290 534 -309 tag 3 COLOR 2 original 11 325 -275 274 134 tag 3 COLOR 2 original 12 326 537 424 134 tag 3 COLOR 2 original 12 327 539 448 137 tag 3 COLOR 2 original 12 328 293 540 -312 tag 3 COLOR 2 original 12 329 -278 277 140 tag 3 COLOR 2 original 13 330 543 426 140 tag 3 COLOR 2 original 13 331 545 450 143 tag 3 COLOR 2 original 13 332 296 546 -315 tag 3 COLOR 2 original 13 333 548 464 146 tag 3 COLOR 2 original 21 334 299 549 -318 tag 3 COLOR 2 original 21 335 -281 280 149 tag 3 COLOR 2 original 21 336 552 440 149 tag 3 COLOR 2 original 21 337 554 466 152 tag 3 COLOR 2 original 22 338 302 555 -321 tag 3 COLOR 2 original 22 339 -284 283 155 tag 3 COLOR 2 original 22 340 558 442 155 tag 3 COLOR 2 original 22 341 560 468 158 tag 3 COLOR 2 original 23 342 305 561 -324 tag 3 COLOR 2 original 23 343 -287 286 161 tag 3 COLOR 2 original 23 344 564 444 161 tag 3 COLOR 2 original 23 bodies /* facets */ 1 1 2 3 4 5 6 7 8 9 10 \ 11 12 13 14 15 16 17 18 19 20 \ 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 \ 41 42 43 44 -45 46 -47 -48 49 -50 \ -51 52 -53 -54 -55 -56 -57 -58 -59 -60 \ 61 -62 -63 64 -65 -66 67 -68 69 -70 \ 71 -72 73 -74 75 -76 77 -78 79 -80 \ 81 -82 83 -84 85 -86 87 -88 89 -90 \ 91 -92 -93 94 -95 96 -97 98 -99 100 \ -101 102 -103 104 105 -106 107 -108 109 -110 \ 111 -112 113 -114 115 -116 117 -118 119 -120 \ 121 -122 123 -124 125 -126 127 -128 -129 -130 \ -131 -132 -133 -134 -135 136 -137 -138 139 -140 \ 141 -142 -143 144 -145 146 -147 -148 149 -150 \ 151 -152 -153 154 -155 156 -157 -158 159 -160 \ 161 -162 -163 164 165 -166 167 -168 169 -170 \ 171 -172 173 -174 175 -176 -177 178 -179 180 \ -181 182 -183 184 -185 186 -187 188 -189 190 \ -191 192 -193 194 -195 196 -197 198 -199 200 \ -201 202 203 -204 -205 206 207 -208 -209 210 \ 211 -212 -213 214 215 -216 -217 218 219 -220 \ -221 222 223 -224 -225 226 -227 228 -229 230 \ -231 232 -233 234 -235 236 -237 238 -239 240 \ -241 242 -243 244 -245 246 -247 248 -249 -250 \ -251 -252 -253 -254 -255 -256 -257 -258 -259 -260 \ -261 -262 -263 -264 -265 -266 -267 -268 -269 -270 \ -271 -272 -273 -274 -275 -276 -277 -278 -279 -280 \ -281 -282 -283 -284 -285 -286 -287 -288 -289 -290 \ -291 -292 -293 -294 -295 -296 -297 298 -299 300 \ -301 302 -303 304 -305 306 -307 308 -309 310 \ -311 312 -313 314 -315 316 -317 318 -319 320 \ 321 -322 323 -324 325 -326 327 -328 329 -330 \ 331 -332 333 -334 335 -336 337 -338 339 -340 \ 341 -342 343 -344 volume meniscus_volume /*actual: 1.5000000893384*/ volconst -6.28318530717959 read rr := { refine edges where not fixed; u } checksmall := {delete edges where length < 0.0001; \ delete facets where area < 0.005; u; \ delete edges where length < 0.0001; u;g3} echeck_old := {refine edges where length > 0.33 ; u ; checksmall} edgemax := 0.4 echeck := {while max(edges,length) > edgemax do \ { refine edges where length > edgemax } ; \ checksmall } biggrind := { ig:= 1 ; change := 1.; olde := total_energy ;\ while (ig < 10000 && abs(change) > .00001) \ do {ig := ig+1; \ if ig%10 == 1 then checksmall ; \ g | "cat > null" ; \ if ig%3 == 1 then \ { change := total_energy - olde ; olde := total_energy;\ printf "sum of 3 dE: %12.10f\n", change } } } smallgrind := { ig:= 1 ; \ while ig < 10 do \ {ig := ig+1; g; change := total_energy - olde ; \ olde := total_energy } } goprint := \ { printf "%f %f %20.12f %f " \ ,phi, theta, total_energy, change | "cat >> results";\ foreach body do { printf "%f ", volume | "cat >> results" } ;\ printf "%f ", total_area | "cat >> results" ;\ printf "%g ", facet_count | "cat >> results" ; \ printf "%f \n", max(edges,length) | "cat >> results" } dump_to_file := {filename := sprintf "dump.v_%f_p_%f_t_%f",vol,phi,theta ; \ dump filename} inc := {echeck; D; D ; O ; o ;u; theta := theta + 1 ; \ printf "%f", theta ;u;g5; biggrind ; \ if theta%10 = 0 then dump_to_file; goprint } run_this := { while theta < 120 do {inc} } //rr //echeck