diff --git a/test/porousmediumflow/2p/implicit/CMakeLists.txt b/test/porousmediumflow/2p/implicit/CMakeLists.txt index 069135d12e5e6bf1150eda5ab6d8bda8df308c8d..8b02698b0cc973b0bd77887187a2dd8d1b86c90a 100644 --- a/test/porousmediumflow/2p/implicit/CMakeLists.txt +++ b/test/porousmediumflow/2p/implicit/CMakeLists.txt @@ -4,3 +4,4 @@ add_subdirectory(cornerpoint) add_subdirectory(fracture) add_subdirectory(incompressible) add_subdirectory(nonisothermal) +add_subdirectory(rotationsymmetry) diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/CMakeLists.txt b/test/porousmediumflow/2p/implicit/rotationsymmetry/CMakeLists.txt new file mode 100644 index 0000000000000000000000000000000000000000..8dd8521ce3a640e5400b483add39e31f51186966 --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/CMakeLists.txt @@ -0,0 +1,11 @@ +dune_symlink_to_source_files(FILES "params.input" "dome.msh") + +dumux_add_test(NAME test_2p_rotationsymmetry_dome + SOURCES main.cc + LABELS porousmediumflow 2p + CMAKE_GUARD dune-alugrid_FOUND + COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py + CMD_ARGS --script fuzzy + --files ${CMAKE_SOURCE_DIR}/test/references/test_2p_rotationsymmetry_dome-reference.vtu + ${CMAKE_CURRENT_BINARY_DIR}/test_2p_rotationsymmetry_dome-00002.vtu + --command "${CMAKE_CURRENT_BINARY_DIR}/test_2p_rotationsymmetry_dome -TimeLoop.TEnd 1e6") diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.geo b/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.geo new file mode 100644 index 0000000000000000000000000000000000000000..00a21b18bb760a55f3bcbaccc4ee1ba156d4f28d --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.geo @@ -0,0 +1,8 @@ +Point(1) = {0, 0, 0, 1.0}; +Point(2) = {0, 100, 0, 1.0}; +Point(3) = {100, 0, 0, 1.0}; +Circle(1) = {2, 1, 3}; +Line(2) = {2, 1}; +Line(3) = {1, 3}; +Curve Loop(1) = {2, 3, -1}; +Plane Surface(1) = {1}; diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.msh b/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.msh new file mode 100644 index 0000000000000000000000000000000000000000..c689d06a4fc10dcaaf182ee7a2d6dd304faed3fc --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/dome.msh @@ -0,0 +1,1264 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +418 +1 0 0 0 +2 0 100 0 +3 100 0 0 +4 4.906767442902359 99.87954562001808 0 +5 9.801714051886403 99.51847266535522 0 +6 14.67304747598722 98.91765099196111 0 +7 19.50903224507124 98.07852803167863 0 +8 24.29801804552329 97.00312530562829 0 +9 29.02846779285695 95.69403355277207 0 +10 33.68898541848564 94.1544064899411 0 +11 38.26834333025619 92.3879532122973 0 +12 42.75550944722013 90.39892926306521 0 +13 47.13967379803581 88.19212637313365 0 +14 51.41027454417063 85.77286092519591 0 +15 55.55702343687956 83.14696114010431 0 +16 59.56993059052324 80.32075304328417 0 +17 63.43932856700083 77.30104521264968 0 +18 67.15589564692182 74.09511238846829 0 +19 70.71067828870152 70.71067794860799 0 +20 74.09511269245844 67.15589531152037 0 +21 77.30104547792956 63.43932824375638 0 +22 80.32075327125473 59.56993028314053 0 +23 83.14696133638893 55.55702314311887 0 +24 85.77286109258662 51.4102742648961 0 +25 88.19212651248949 47.13967353731935 0 +26 90.39892937667582 42.75550920701072 0 +27 92.38795330229983 38.2683431129709 0 +28 94.15440655902738 33.68898522540241 0 +29 95.69403360237926 29.02846762932397 0 +30 97.00312533880519 24.29801791307371 0 +31 98.07852805328298 19.50903213645886 0 +32 98.91765100396263 14.6730473950796 0 +33 99.51847267037007 9.801714000969689 0 +34 99.87954562120397 4.906767418762911 0 +35 0 95.00000000003553 0 +36 0 90.00000000007105 0 +37 0 85.00000000010658 0 +38 0 80.00000000014211 0 +39 0 75.00000000017764 0 +40 0 70.00000000021316 0 +41 0 65.00000000024869 0 +42 0 60.00000000025813 0 +43 0 55.00000000025813 0 +44 0 50.00000000025813 0 +45 0 45.00000000025813 0 +46 0 40.00000000025814 0 +47 0 35.0000000002487 0 +48 0 30.00000000021316 0 +49 0 25.00000000017765 0 +50 0 20.00000000014211 0 +51 0 15.00000000010658 0 +52 0 10.00000000007105 0 +53 0 5.000000000035527 0 +54 4.999999999979755 0 0 +55 9.999999999961581 0 0 +56 14.99999999994008 0 0 +57 19.99999999991413 0 0 +58 24.99999999988819 0 0 +59 29.99999999986224 0 0 +60 34.99999999984823 0 0 +61 39.99999999984005 0 0 +62 44.99999999983187 0 0 +63 49.99999999982369 0 0 +64 54.99999999981551 0 0 +65 59.99999999980733 0 0 +66 64.99999999980859 0 0 +67 69.99999999983592 0 0 +68 74.99999999986328 0 0 +69 79.99999999989062 0 0 +70 84.99999999991796 0 0 +71 89.9999999999453 0 0 +72 94.99999999997266 0 0 +73 65.74764107637027 70.07739152843438 0 +74 42.49999999983596 4.330127018915107 0 +75 4.329610735593052 42.50124369698572 0 +76 81.06893144331036 51.28834902471475 0 +77 47.26779913280121 82.68970037615367 0 +78 62.49722997421908 4.328884071679592 0 +79 4.33012701890272 62.50000000025341 0 +80 27.29767338392531 4.213264007515257 0 +81 89.4855067614475 34.74757790418513 0 +82 4.669985493878081 26.77382104369122 0 +83 30.03215592899558 90.97009880298751 0 +84 77.63286468266325 5.202060439815369 0 +85 4.32927840926416 77.50924710018683 0 +86 93.88135826884982 21.26065778235738 0 +87 16.36433806022532 94.31000385062856 0 +88 52.50129042238162 4.293189013966979 0 +89 4.328872519847133 52.59228124833739 0 +90 17.50079587705269 4.330082909154775 0 +91 4.445100235827688 17.45352014180569 0 +92 59.42419838353013 75.28859411333592 0 +93 73.17300819681915 62.33247967861506 0 +94 38.64026815173551 87.40345033088728 0 +95 85.47490603781043 42.96807886433079 0 +96 87.36865184104705 3.916335428142458 0 +97 4.279210640981523 87.4195603816632 0 +98 95.0156209575843 11.83714149663033 0 +99 4.529962364585987 7.959337617259965 0 +100 7.041542396036847 95.45986102913797 0 +101 94.28456122425735 5.270244443758822 0 +102 4.316107764688089 82.48940523130123 0 +103 8.640588398377419 80.01186444862348 0 +104 8.653872498676236 75.02756608980934 0 +105 12.95535640490645 77.5285657904296 0 +106 12.9692534815371 72.55279725433547 0 +107 17.26275297620415 75.04673797180237 0 +108 17.28099487024959 70.08603133499213 0 +109 21.56702704636393 72.572647466658 0 +110 21.58552871737995 67.61737812528918 0 +111 25.87799976615322 70.08839122550148 0 +112 25.88343226160008 65.14192912999911 0 +113 30.1769107163653 67.61555513166105 0 +114 30.18180113179168 62.6733460051151 0 +115 34.49682343608109 65.14863897781237 0 +116 34.47608246845549 60.21018296930308 0 +117 30.1917206301029 57.72032042627465 0 +118 34.49098290180054 55.25332503180154 0 +119 30.20341663530289 52.76288770533814 0 +120 34.50552522680265 50.2938646462994 0 +121 38.82052927169608 52.76892270325682 0 +122 38.82381267801115 47.81793695889353 0 +123 34.52135203845079 45.33066720941083 0 +124 38.82819650996136 42.86351335388304 0 +125 34.53734965606475 40.36641725124647 0 +126 43.13040868830476 45.34721920820142 0 +127 43.13145091893426 40.39863838329646 0 +128 47.42245663479966 42.88998642096508 0 +129 47.43642938294413 37.93336994694638 0 +130 51.72520791048775 40.4236886896651 0 +131 51.73995344756116 35.46894222492031 0 +132 47.45529633456669 32.97299833344081 0 +133 56.02735377240293 37.95802111253326 0 +134 56.04234314039306 33.00478434423913 0 +135 51.73235872889467 45.36473797279909 0 +136 60.32683456446862 35.49692040366982 0 +137 60.34317482184373 30.5409414383875 0 +138 64.62502139360612 33.03550272072999 0 +139 64.64455396513011 28.07953498099611 0 +140 60.35852210326364 25.58202264331618 0 +141 64.60773395819372 37.98486168674243 0 +142 30.22483376827028 42.83740625162804 0 +143 30.23841209075118 37.8688461083787 0 +144 38.81771362748392 62.66519893387589 0 +145 21.59254986772645 62.67359588016453 0 +146 25.92506910462838 40.34184675154989 0 +147 25.93830986323014 35.366993346329 0 +148 30.25894453426202 32.89157995866493 0 +149 25.953621552062 30.38919453952801 0 +150 30.2755013290433 27.9129284906764 0 +151 25.96862860404781 25.40822909270841 0 +152 30.29153490596536 22.93178073993816 0 +153 64.6643656961804 23.11757894598578 0 +154 68.95142541992786 25.61516312266223 0 +155 68.90109007412715 35.52736341413807 0 +156 21.64548094190252 32.86768180576562 0 +157 34.5616187354652 70.05318790647904 0 +158 68.97121880444675 20.65367636032717 0 +159 73.25859053168384 23.14913311758746 0 +160 60.30256476871432 20.74129787294756 0 +161 25.98365924331013 20.4250474559013 0 +162 21.65802750679567 22.90349155766883 0 +163 25.89412687297853 75.0211934529064 0 +164 68.88528692493487 40.47098291199519 0 +165 64.59721709972784 42.92907571667419 0 +166 34.5983199980419 25.4372416932837 0 +167 34.61528844723193 20.455433187418 0 +168 38.91435114902745 22.97276341683611 0 +169 56.0527634478957 23.10669605886711 0 +170 17.24987407365686 80.0181583352759 0 +171 68.86621609925928 45.41130222350336 0 +172 64.58534310708261 47.86622573372247 0 +173 73.27791298000245 18.19141890587875 0 +174 77.56096933993875 20.68712963902767 0 +175 77.53356956738668 25.6472371390876 0 +176 55.99392493671812 18.26287100377937 0 +177 21.69044790280134 17.94703327184519 0 +178 38.92632293577518 18.005131769814 0 +179 73.17051715368531 38.01527111090213 0 +180 17.52192371888526 20.40273195369127 0 +181 34.66152524192101 15.53012203790597 0 +182 51.75882336578965 20.61608175692252 0 +183 69.54321269493596 16.42541123447607 0 +184 68.85138419634458 50.34365178953744 0 +185 64.57301221640233 52.79834284268286 0 +186 51.87005557764096 15.8255694460083 0 +187 77.58151030569566 15.72879849639907 0 +188 81.8637001490448 18.22445942944843 0 +189 17.51587216857569 15.6527413058495 0 +190 56.82194829175987 13.18422510146814 0 +191 81.84100777787486 13.34145554614948 0 +192 12.98474976450527 67.58223606211624 0 +193 60.30617478363577 50.32123722015436 0 +194 60.29267946633962 55.25672972992788 0 +195 17.35948301706967 25.3866172348808 0 +196 25.91432996835665 45.31089289670485 0 +197 21.60968472455217 42.81730603649814 0 +198 21.59790186910651 47.7905059035465 0 +199 17.30289416364726 45.29323720547636 0 +200 17.27925245768088 50.27306735269514 0 +201 26.22297600372288 15.75781855479648 0 +202 17.328294004548 40.26683201758286 0 +203 56.02698079697925 52.77928147908927 0 +204 56.01248631144276 57.71491317753299 0 +205 51.74555319155682 55.23640690223752 0 +206 51.73404832430953 60.17110408819899 0 +207 55.99914231226552 62.64920276773076 0 +208 51.72337204407457 65.10459812430139 0 +209 47.46312936551016 57.69136861418507 0 +210 25.90153166144632 55.23859957533934 0 +211 38.89787040099515 27.94393731785816 0 +212 43.20724718170744 25.49352967837009 0 +213 12.99451589576677 47.72753483229309 0 +214 12.97349960247049 52.74619927427995 0 +215 12.64163987569352 22.32574500438168 0 +216 13.0307180636107 27.74167566401164 0 +217 64.5307821026702 57.68177834214084 0 +218 86.10173404010602 15.69844367247996 0 +219 38.84745348665648 13.22505312165234 0 +220 34.64843587343802 10.48419511339939 0 +221 73.00811196017266 47.88953131758451 0 +222 72.96077003504584 52.80982880181089 0 +223 73.19689215142554 33.07019996349504 0 +224 77.45743663731901 35.56229784123614 0 +225 77.33362413788674 40.49849268834229 0 +226 47.47394741982909 52.75119612862419 0 +227 51.75637778843235 25.56295626781408 0 +228 17.27567590672354 55.24471522224147 0 +229 12.96844718185902 57.70656835847542 0 +230 55.95142792497715 67.51669744675027 0 +231 51.67219591997037 69.96147558334161 0 +232 47.44309908333406 67.54800986230839 0 +233 47.43173719954036 72.47566491365305 0 +234 43.17739317086539 70.01887495440815 0 +235 43.16462577926845 74.94494372945672 0 +236 47.49567843667323 18.11426164437076 0 +237 47.50080787932236 13.00274038383438 0 +238 86.14685605935694 20.71484994225121 0 +239 51.61839054970432 75.04315169136586 0 +240 60.13654755904322 64.8988892443184 0 +241 82.29044762479515 38.49898198775259 0 +242 81.73919483500615 33.11680037232429 0 +243 85.9172632043704 11.11791658516535 0 +244 77.62577677528989 10.76201567767068 0 +245 60.55745327251269 69.90923766853741 0 +246 8.660545076515575 55.16344797469721 0 +247 38.89810271021763 72.49855065807584 0 +248 38.89797072757375 77.41652768413617 0 +249 43.070978753521 79.96595733822934 0 +250 38.85662253663937 82.35031400283718 0 +251 34.62389689053852 79.8753210944318 0 +252 34.61737559638173 84.78863738791047 0 +253 30.32405534835639 82.27795234401098 0 +254 30.25300433719519 77.44973537230723 0 +255 39.53141469134343 8.246997760087364 0 +256 34.81397781951908 5.803735709222786 0 +257 43.16878831128148 65.09893945233412 0 +258 43.14507045148686 55.23121575518452 0 +259 56.03852589264197 47.84052467403419 0 +260 4.332911793485893 72.5175950803992 0 +261 4.332431201837555 67.51160896223899 0 +262 8.665458256500978 65.07204565003953 0 +263 73.23993361046197 28.1064722036555 0 +264 81.83926112130763 23.18078196202001 0 +265 86.05661737788525 25.76303739466323 0 +266 81.79806952936998 28.14876039775954 0 +267 86.57325226943968 30.9961737839952 0 +268 90.4811502938928 27.91189120038829 0 +269 8.626813138971411 84.98446155817533 0 +270 8.871573180563942 90.59022479940126 0 +271 12.93861975228699 87.48445340508485 0 +272 12.94341839157717 82.50769912595965 0 +273 17.26919313454631 84.94029503036117 0 +274 21.54149460191296 82.49431391130604 0 +275 22.19992475414031 88.134469712591 0 +276 25.96660177335635 84.87366708645982 0 +277 25.92020729820019 79.94239491673613 0 +278 82.70655500598889 4.22277106080598 0 +279 47.37403472576841 77.50052899822427 0 +280 51.49221661078138 80.4223379061696 0 +281 68.77274093600489 55.20422115728036 0 +282 73.11411638969841 42.95251662288733 0 +283 76.86698076098473 45.41474480071059 0 +284 60.25325570141528 60.14391862046401 0 +285 64.47730302903805 62.51752355179332 0 +286 57.5001958552845 4.287162030131886 0 +287 59.88454548249744 8.64592589550001 0 +288 64.97144316076027 8.47247428496201 0 +289 62.49406105670936 12.97261483043062 0 +290 67.48770406055877 4.300252661153177 0 +291 69.98688334662016 8.409983187581334 0 +292 72.47740129132599 4.341191207251437 0 +293 60.31427289694981 40.44547338338117 0 +294 34.58042337641118 30.41634276857175 0 +295 38.87080275264356 32.92453220369988 0 +296 47.50060944040172 4.322056909026594 0 +297 49.98962347164866 8.410032656232989 0 +298 44.94359464651194 8.342298228498993 0 +299 34.55850761012051 35.39454149167496 0 +300 4.33838421479682 47.55750648939694 0 +301 8.663891220570918 45.03623195500357 0 +302 8.670865920514558 40.05859630274539 0 +303 4.332614033575331 37.51278334224497 0 +304 8.662617807733483 35.04641078422374 0 +305 13.0330512694825 37.70025488989715 0 +306 4.371806070467909 32.3888312495379 0 +307 38.82169584487648 57.71285312152594 0 +308 68.92648610244646 30.57237273427949 0 +309 22.46634943595914 4.31026004921147 0 +310 20.00062990747194 8.659265610237508 0 +311 15.00205668538001 8.423800593832368 0 +312 12.54215616321364 4.267438943148594 0 +313 9.971816199266264 8.469956762080136 0 +314 13.08346492033741 13.18835413903467 0 +315 7.564905907339966 12.81490883230902 0 +316 9.853704409774153 17.56834745913535 0 +317 7.322067298700945 21.84989356940932 0 +318 24.76734896193813 8.462520351099283 0 +319 30.02554064857317 8.447275617737059 0 +320 21.59895248406224 52.75551323314012 0 +321 17.26663122185808 60.21550724096004 0 +322 30.21494366805078 47.80103655796881 0 +323 17.36112929046558 30.31982742183579 0 +324 30.33542219261327 17.99172647384007 0 +325 21.57258047821944 77.5159076757808 0 +326 21.65772848539054 27.87917360873124 0 +327 21.64017550938289 37.83445454583325 0 +328 25.8958932315126 60.19329861534246 0 +329 8.664763616065917 70.04884522776555 0 +330 77.49418272182822 30.60862798625968 0 +331 56.0598380838785 28.04865196626943 0 +332 47.45203891280764 62.62704891284362 0 +333 60.31924108777795 45.38898532351842 0 +334 43.16702049941901 60.15659359663918 0 +335 64.74895775813718 17.6558989680275 0 +336 38.84396828643512 37.90476482784889 0 +337 25.90517938105423 50.27657264533961 0 +338 56.02213889103706 42.90218677548332 0 +339 43.15224162307582 35.4378343824732 0 +340 72.81448561878584 12.80429937064667 0 +341 51.75693664363422 30.51653325684634 0 +342 43.17854489332703 30.4689497200394 0 +343 51.73990888527399 50.29916591542062 0 +344 47.43679833180885 47.82301091498905 0 +345 43.1502419749551 50.28752281500319 0 +346 21.58853906222154 57.72020496119799 0 +347 47.48085309615976 28.02101556132072 0 +348 17.29146773680154 65.14414937889995 0 +349 38.85340666523244 67.58056514716421 0 +350 72.89010109052971 57.53205163450111 0 +351 77.63801694498021 55.25319404376224 0 +352 47.49885189214582 23.06170976272566 0 +353 34.57266030184984 74.96750960998897 0 +354 4.334432580230314 57.56208511394772 0 +355 8.659951421661283 50.13720029566802 0 +356 12.97730083794022 62.6395480250622 0 +357 43.21068471390196 20.54250097572984 0 +358 30.22272012166788 72.53259544980735 0 +359 13.02635825629062 42.73565503840436 0 +360 43.10334330814145 15.7719204230498 0 +361 90.37944169968721 13.36366051182722 0 +362 90.686357822248 8.090482484035171 0 +363 17.34925090044651 35.27991851500461 0 +364 8.565619087620563 29.99997232660565 0 +365 8.664908694122419 60.10999650311816 0 +366 13.02524820808494 32.68394816469641 0 +367 55.00265356557755 8.447509116850913 0 +368 22.10648788207682 12.83849508407011 0 +369 42.95514196999216 85.16674628071772 0 +370 68.66502869976063 65.57326477220529 0 +371 87.71876369341878 38.88305359235953 0 +372 83.60097907913558 47.2642313465323 0 +373 34.32589646154753 89.38807738838014 0 +374 25.59118761813957 92.47466935611492 0 +375 20.89319622428411 93.28335780543536 0 +376 17.17899435776213 89.94559290493416 0 +377 11.74534950769584 95.22877780251964 0 +378 94.50452951873228 16.40489240658113 0 +379 7.654236161452737 4.096638941984686 0 +380 68.67984172554959 59.96475267009521 0 +381 30.51133473185923 13.29229023833617 0 +382 55.76279789469748 72.13679577879026 0 +383 3.973754479272422 92.35670409051875 0 +384 90.18026354541556 18.54026640581831 0 +385 81.64418079406788 8.759235097401891 0 +386 80.92945177428504 43.74135498904107 0 +387 90.0526736420735 23.62578288999698 0 +388 30.33762622837198 86.84152579437286 0 +389 55.57075957768602 77.72643227884168 0 +390 9.24600596390076 25.73822152161991 0 +391 76.79002758854654 50.2197451676239 0 +392 85.22193164860604 35.48873841432025 0 +393 60.04104456053673 16.51316902596032 0 +394 63.21285348038026 73.27602967870479 0 +395 74.10748234293703 8.303909976593097 0 +396 3.307993701550728 12.64555331831046 0 +397 38.36907850207734 3.676172097645051 0 +398 13.41977497170684 91.51181055251369 0 +399 14.11105409719031 17.84594513373735 0 +400 3.265649226084689 22.50000000015988 0 +401 67.46844156850676 12.12597601391473 0 +402 26.82549926060076 88.65888615050523 0 +403 76.36947525299063 59.63865655034557 0 +404 52.43001417990796 11.73998363237946 0 +405 93.42246822920011 25.22496348302806 0 +406 3.660254037864924 3.660254037879868 0 +407 42.63911239651951 11.74183464289752 0 +408 31.29012826660029 3.631607248023507 0 +409 17.50165397260112 11.57134348584253 0 +410 63.85803716006139 66.45842255892931 0 +411 85.49615100016553 7.217997416425657 0 +412 3.441333151781442 96.44631184511105 0 +413 27.04648625742664 11.81874300882031 0 +414 91.17208348027135 3.563454058842772 0 +415 91.33595260406413 31.55519826095835 0 +416 79.92890623170726 47.69241903185679 0 +417 96.9079429549133 2.982263075779672 0 +418 96.21848417470137 8.107766871032911 0 +$EndNodes +$Elements +837 +1 15 2 0 1 1 +2 15 2 0 2 2 +3 15 2 0 3 3 +4 1 2 0 1 2 4 +5 1 2 0 1 4 5 +6 1 2 0 1 5 6 +7 1 2 0 1 6 7 +8 1 2 0 1 7 8 +9 1 2 0 1 8 9 +10 1 2 0 1 9 10 +11 1 2 0 1 10 11 +12 1 2 0 1 11 12 +13 1 2 0 1 12 13 +14 1 2 0 1 13 14 +15 1 2 0 1 14 15 +16 1 2 0 1 15 16 +17 1 2 0 1 16 17 +18 1 2 0 1 17 18 +19 1 2 0 1 18 19 +20 1 2 0 1 19 20 +21 1 2 0 1 20 21 +22 1 2 0 1 21 22 +23 1 2 0 1 22 23 +24 1 2 0 1 23 24 +25 1 2 0 1 24 25 +26 1 2 0 1 25 26 +27 1 2 0 1 26 27 +28 1 2 0 1 27 28 +29 1 2 0 1 28 29 +30 1 2 0 1 29 30 +31 1 2 0 1 30 31 +32 1 2 0 1 31 32 +33 1 2 0 1 32 33 +34 1 2 0 1 33 34 +35 1 2 0 1 34 3 +36 1 2 0 2 2 35 +37 1 2 0 2 35 36 +38 1 2 0 2 36 37 +39 1 2 0 2 37 38 +40 1 2 0 2 38 39 +41 1 2 0 2 39 40 +42 1 2 0 2 40 41 +43 1 2 0 2 41 42 +44 1 2 0 2 42 43 +45 1 2 0 2 43 44 +46 1 2 0 2 44 45 +47 1 2 0 2 45 46 +48 1 2 0 2 46 47 +49 1 2 0 2 47 48 +50 1 2 0 2 48 49 +51 1 2 0 2 49 50 +52 1 2 0 2 50 51 +53 1 2 0 2 51 52 +54 1 2 0 2 52 53 +55 1 2 0 2 53 1 +56 1 2 0 3 1 54 +57 1 2 0 3 54 55 +58 1 2 0 3 55 56 +59 1 2 0 3 56 57 +60 1 2 0 3 57 58 +61 1 2 0 3 58 59 +62 1 2 0 3 59 60 +63 1 2 0 3 60 61 +64 1 2 0 3 61 62 +65 1 2 0 3 62 63 +66 1 2 0 3 63 64 +67 1 2 0 3 64 65 +68 1 2 0 3 65 66 +69 1 2 0 3 66 67 +70 1 2 0 3 67 68 +71 1 2 0 3 68 69 +72 1 2 0 3 69 70 +73 1 2 0 3 70 71 +74 1 2 0 3 71 72 +75 1 2 0 3 72 3 +76 2 2 0 1 335 289 401 +77 2 2 0 1 183 335 401 +78 2 2 0 1 275 374 375 +79 2 2 0 1 370 285 380 +80 2 2 0 1 190 287 289 +81 2 2 0 1 95 241 371 +82 2 2 0 1 244 84 385 +83 2 2 0 1 316 215 317 +84 2 2 0 1 99 313 315 +85 2 2 0 1 374 275 402 +86 2 2 0 1 239 382 389 +87 2 2 0 1 100 270 377 +88 2 2 0 1 82 317 390 +89 2 2 0 1 14 77 280 +90 2 2 0 1 275 375 376 +91 2 2 0 1 317 82 400 +92 2 2 0 1 77 279 280 +93 2 2 0 1 84 244 395 +94 2 2 0 1 93 370 380 +95 2 2 0 1 190 289 393 +96 2 2 0 1 241 95 386 +97 2 2 0 1 99 315 396 +98 2 2 0 1 76 23 351 +99 2 2 0 1 382 92 389 +100 2 2 0 1 215 316 399 +101 2 2 0 1 77 249 279 +102 2 2 0 1 285 370 410 +103 2 2 0 1 238 384 387 +104 2 2 0 1 255 74 298 +105 2 2 0 1 15 14 280 +106 2 2 0 1 23 22 351 +107 2 2 0 1 91 315 316 +108 2 2 0 1 291 340 401 +109 2 2 0 1 377 270 398 +110 2 2 0 1 84 278 385 +111 2 2 0 1 91 316 317 +112 2 2 0 1 287 190 367 +113 2 2 0 1 13 12 369 +114 2 2 0 1 371 241 392 +115 2 2 0 1 255 298 407 +116 2 2 0 1 384 86 387 +117 2 2 0 1 98 362 418 +118 2 2 0 1 313 99 379 +119 2 2 0 1 256 220 319 +120 2 2 0 1 340 183 401 +121 2 2 0 1 98 361 362 +122 2 2 0 1 270 100 383 +123 2 2 0 1 73 245 410 +124 2 2 0 1 245 73 394 +125 2 2 0 1 302 305 359 +126 2 2 0 1 245 92 382 +127 2 2 0 1 253 276 277 +128 2 2 0 1 163 254 277 +129 2 2 0 1 254 253 277 +130 2 2 0 1 300 301 355 +131 2 2 0 1 249 77 369 +132 2 2 0 1 305 202 359 +133 2 2 0 1 317 215 390 +134 2 2 0 1 77 13 369 +135 2 2 0 1 12 94 369 +136 2 2 0 1 310 318 368 +137 2 2 0 1 157 115 349 +138 2 2 0 1 301 213 355 +139 2 2 0 1 340 291 395 +140 2 2 0 1 250 94 252 +141 2 2 0 1 304 364 366 +142 2 2 0 1 362 101 418 +143 2 2 0 1 289 335 393 +144 2 2 0 1 247 157 349 +145 2 2 0 1 297 237 298 +146 2 2 0 1 257 144 334 +147 2 2 0 1 144 307 334 +148 2 2 0 1 121 258 307 +149 2 2 0 1 307 258 334 +150 2 2 0 1 242 266 267 +151 2 2 0 1 92 245 394 +152 2 2 0 1 301 302 359 +153 2 2 0 1 254 163 358 +154 2 2 0 1 256 319 408 +155 2 2 0 1 370 73 410 +156 2 2 0 1 89 300 355 +157 2 2 0 1 89 246 354 +158 2 2 0 1 262 192 329 +159 2 2 0 1 266 242 330 +160 2 2 0 1 261 262 329 +161 2 2 0 1 175 266 330 +162 2 2 0 1 121 122 345 +163 2 2 0 1 258 121 345 +164 2 2 0 1 76 351 391 +165 2 2 0 1 321 228 346 +166 2 2 0 1 163 111 358 +167 2 2 0 1 368 318 413 +168 2 2 0 1 115 144 349 +169 2 2 0 1 353 254 358 +170 2 2 0 1 375 87 376 +171 2 2 0 1 105 170 272 +172 2 2 0 1 276 253 388 +173 2 2 0 1 202 305 363 +174 2 2 0 1 43 89 354 +175 2 2 0 1 272 170 273 +176 2 2 0 1 271 272 273 +177 2 2 0 1 213 301 359 +178 2 2 0 1 68 84 292 +179 2 2 0 1 163 277 325 +180 2 2 0 1 154 263 308 +181 2 2 0 1 263 223 308 +182 2 2 0 1 103 105 272 +183 2 2 0 1 259 135 338 +184 2 2 0 1 364 216 366 +185 2 2 0 1 145 321 346 +186 2 2 0 1 315 91 396 +187 2 2 0 1 109 163 325 +188 2 2 0 1 188 238 264 +189 2 2 0 1 266 265 267 +190 2 2 0 1 41 79 261 +191 2 2 0 1 260 40 261 +192 2 2 0 1 264 238 265 +193 2 2 0 1 320 198 337 +194 2 2 0 1 246 89 355 +195 2 2 0 1 264 265 266 +196 2 2 0 1 261 79 262 +197 2 2 0 1 102 85 103 +198 2 2 0 1 39 40 260 +199 2 2 0 1 157 353 358 +200 2 2 0 1 250 249 369 +201 2 2 0 1 102 103 269 +202 2 2 0 1 263 175 330 +203 2 2 0 1 179 225 282 +204 2 2 0 1 97 269 270 +205 2 2 0 1 260 261 329 +206 2 2 0 1 135 130 338 +207 2 2 0 1 282 225 283 +208 2 2 0 1 270 269 271 +209 2 2 0 1 157 247 353 +210 2 2 0 1 38 85 102 +211 2 2 0 1 273 170 274 +212 2 2 0 1 194 217 284 +213 2 2 0 1 221 282 283 +214 2 2 0 1 284 217 285 +215 2 2 0 1 97 102 269 +216 2 2 0 1 240 284 285 +217 2 2 0 1 84 69 278 +218 2 2 0 1 210 320 337 +219 2 2 0 1 67 68 292 +220 2 2 0 1 269 103 272 +221 2 2 0 1 144 257 349 +222 2 2 0 1 40 41 261 +223 2 2 0 1 343 226 344 +224 2 2 0 1 267 265 268 +225 2 2 0 1 279 239 280 +226 2 2 0 1 174 188 264 +227 2 2 0 1 164 179 282 +228 2 2 0 1 139 154 308 +229 2 2 0 1 273 274 275 +230 2 2 0 1 188 218 238 +231 2 2 0 1 233 239 279 +232 2 2 0 1 143 148 299 +233 2 2 0 1 105 107 170 +234 2 2 0 1 217 185 281 +235 2 2 0 1 97 37 102 +236 2 2 0 1 103 104 105 +237 2 2 0 1 36 37 97 +238 2 2 0 1 141 165 293 +239 2 2 0 1 48 82 306 +240 2 2 0 1 204 194 284 +241 2 2 0 1 159 174 175 +242 2 2 0 1 103 85 104 +243 2 2 0 1 175 264 266 +244 2 2 0 1 66 67 290 +245 2 2 0 1 136 141 293 +246 2 2 0 1 70 96 278 +247 2 2 0 1 175 174 264 +248 2 2 0 1 148 294 299 +249 2 2 0 1 125 143 299 +250 2 2 0 1 74 296 298 +251 2 2 0 1 154 159 263 +252 2 2 0 1 70 71 96 +253 2 2 0 1 275 274 276 +254 2 2 0 1 88 64 286 +255 2 2 0 1 109 111 163 +256 2 2 0 1 105 106 107 +257 2 2 0 1 296 297 298 +258 2 2 0 1 223 263 330 +259 2 2 0 1 256 60 397 +260 2 2 0 1 85 39 260 +261 2 2 0 1 141 155 164 +262 2 2 0 1 290 67 292 +263 2 2 0 1 74 62 296 +264 2 2 0 1 179 224 225 +265 2 2 0 1 144 116 307 +266 2 2 0 1 113 115 157 +267 2 2 0 1 136 138 141 +268 2 2 0 1 287 288 289 +269 2 2 0 1 225 224 241 +270 2 2 0 1 164 155 179 +271 2 2 0 1 79 354 365 +272 2 2 0 1 184 221 222 +273 2 2 0 1 274 170 325 +274 2 2 0 1 323 195 326 +275 2 2 0 1 141 138 155 +276 2 2 0 1 230 240 245 +277 2 2 0 1 148 150 294 +278 2 2 0 1 38 39 85 +279 2 2 0 1 286 78 287 +280 2 2 0 1 107 108 109 +281 2 2 0 1 165 164 171 +282 2 2 0 1 122 126 345 +283 2 2 0 1 128 130 135 +284 2 2 0 1 105 104 106 +285 2 2 0 1 184 171 221 +286 2 2 0 1 235 233 279 +287 2 2 0 1 208 207 230 +288 2 2 0 1 287 78 288 +289 2 2 0 1 248 235 249 +290 2 2 0 1 135 343 344 +291 2 2 0 1 172 171 184 +292 2 2 0 1 47 48 306 +293 2 2 0 1 277 274 325 +294 2 2 0 1 68 69 84 +295 2 2 0 1 361 98 378 +296 2 2 0 1 133 136 293 +297 2 2 0 1 296 88 297 +298 2 2 0 1 159 175 263 +299 2 2 0 1 109 110 111 +300 2 2 0 1 249 235 279 +301 2 2 0 1 174 187 188 +302 2 2 0 1 230 207 240 +303 2 2 0 1 233 231 239 +304 2 2 0 1 159 173 174 +305 2 2 0 1 89 44 300 +306 2 2 0 1 194 185 217 +307 2 2 0 1 132 339 342 +308 2 2 0 1 107 109 325 +309 2 2 0 1 248 249 250 +310 2 2 0 1 206 204 207 +311 2 2 0 1 65 78 286 +312 2 2 0 1 37 38 102 +313 2 2 0 1 179 223 224 +314 2 2 0 1 111 112 113 +315 2 2 0 1 185 184 281 +316 2 2 0 1 156 323 326 +317 2 2 0 1 166 168 211 +318 2 2 0 1 300 75 301 +319 2 2 0 1 41 42 79 +320 2 2 0 1 294 211 295 +321 2 2 0 1 154 158 159 +322 2 2 0 1 208 230 231 +323 2 2 0 1 63 88 296 +324 2 2 0 1 113 114 115 +325 2 2 0 1 65 66 78 +326 2 2 0 1 271 269 272 +327 2 2 0 1 69 70 278 +328 2 2 0 1 288 290 291 +329 2 2 0 1 141 164 165 +330 2 2 0 1 188 191 218 +331 2 2 0 1 46 47 303 +332 2 2 0 1 107 106 108 +333 2 2 0 1 136 137 138 +334 2 2 0 1 150 166 294 +335 2 2 0 1 165 171 172 +336 2 2 0 1 302 304 305 +337 2 2 0 1 115 116 144 +338 2 2 0 1 206 207 208 +339 2 2 0 1 214 200 228 +340 2 2 0 1 299 295 336 +341 2 2 0 1 211 168 212 +342 2 2 0 1 116 118 307 +343 2 2 0 1 63 64 88 +344 2 2 0 1 184 222 281 +345 2 2 0 1 128 129 130 +346 2 2 0 1 133 134 136 +347 2 2 0 1 172 184 185 +348 2 2 0 1 171 164 282 +349 2 2 0 1 126 127 128 +350 2 2 0 1 130 131 133 +351 2 2 0 1 166 211 294 +352 2 2 0 1 293 333 338 +353 2 2 0 1 61 62 74 +354 2 2 0 1 294 295 299 +355 2 2 0 1 318 80 319 +356 2 2 0 1 301 75 302 +357 2 2 0 1 207 204 284 +358 2 2 0 1 43 44 89 +359 2 2 0 1 303 47 306 +360 2 2 0 1 118 121 307 +361 2 2 0 1 118 120 121 +362 2 2 0 1 109 108 110 +363 2 2 0 1 200 198 320 +364 2 2 0 1 78 66 290 +365 2 2 0 1 90 309 310 +366 2 2 0 1 310 309 318 +367 2 2 0 1 64 65 286 +368 2 2 0 1 45 75 300 +369 2 2 0 1 309 80 318 +370 2 2 0 1 55 56 312 +371 2 2 0 1 122 124 126 +372 2 2 0 1 333 259 338 +373 2 2 0 1 62 63 296 +374 2 2 0 1 139 153 154 +375 2 2 0 1 241 224 242 +376 2 2 0 1 111 110 112 +377 2 2 0 1 44 45 300 +378 2 2 0 1 166 167 168 +379 2 2 0 1 90 57 309 +380 2 2 0 1 251 250 252 +381 2 2 0 1 221 171 282 +382 2 2 0 1 90 310 311 +383 2 2 0 1 179 155 223 +384 2 2 0 1 45 46 75 +385 2 2 0 1 214 228 229 +386 2 2 0 1 113 112 114 +387 2 2 0 1 227 169 331 +388 2 2 0 1 58 59 80 +389 2 2 0 1 247 235 248 +390 2 2 0 1 302 303 304 +391 2 2 0 1 58 80 309 +392 2 2 0 1 150 152 166 +393 2 2 0 1 313 311 314 +394 2 2 0 1 182 169 227 +395 2 2 0 1 234 233 235 +396 2 2 0 1 232 231 233 +397 2 2 0 1 240 207 284 +398 2 2 0 1 56 90 312 +399 2 2 0 1 75 46 303 +400 2 2 0 1 315 314 316 +401 2 2 0 1 115 114 116 +402 2 2 0 1 56 57 90 +403 2 2 0 1 228 200 320 +404 2 2 0 1 135 259 343 +405 2 2 0 1 57 58 309 +406 2 2 0 1 138 139 308 +407 2 2 0 1 48 49 82 +408 2 2 0 1 276 274 277 +409 2 2 0 1 50 51 91 +410 2 2 0 1 52 53 99 +411 2 2 0 1 339 295 342 +412 2 2 0 1 312 311 313 +413 2 2 0 1 203 194 204 +414 2 2 0 1 121 120 122 +415 2 2 0 1 128 127 129 +416 2 2 0 1 193 185 194 +417 2 2 0 1 288 78 290 +418 2 2 0 1 205 204 206 +419 2 2 0 1 90 311 312 +420 2 2 0 1 313 314 315 +421 2 2 0 1 126 124 127 +422 2 2 0 1 180 162 195 +423 2 2 0 1 291 290 292 +424 2 2 0 1 130 129 131 +425 2 2 0 1 302 75 303 +426 2 2 0 1 133 131 134 +427 2 2 0 1 136 134 137 +428 2 2 0 1 166 152 167 +429 2 2 0 1 195 162 326 +430 2 2 0 1 168 167 178 +431 2 2 0 1 209 332 334 +432 2 2 0 1 229 228 321 +433 2 2 0 1 197 196 198 +434 2 2 0 1 130 133 338 +435 2 2 0 1 215 195 216 +436 2 2 0 1 304 303 306 +437 2 2 0 1 199 198 200 +438 2 2 0 1 138 137 139 +439 2 2 0 1 174 173 187 +440 2 2 0 1 159 158 173 +441 2 2 0 1 146 142 196 +442 2 2 0 1 170 107 325 +443 2 2 0 1 293 165 333 +444 2 2 0 1 228 320 346 +445 2 2 0 1 248 250 251 +446 2 2 0 1 143 147 148 +447 2 2 0 1 208 231 232 +448 2 2 0 1 154 153 158 +449 2 2 0 1 104 85 260 +450 2 2 0 1 161 177 201 +451 2 2 0 1 180 195 215 +452 2 2 0 1 198 196 337 +453 2 2 0 1 336 295 339 +454 2 2 0 1 216 195 323 +455 2 2 0 1 213 200 214 +456 2 2 0 1 172 185 193 +457 2 2 0 1 148 149 150 +458 2 2 0 1 21 20 93 +459 2 2 0 1 251 252 253 +460 2 2 0 1 150 151 152 +461 2 2 0 1 17 16 92 +462 2 2 0 1 118 119 120 +463 2 2 0 1 116 117 118 +464 2 2 0 1 234 235 247 +465 2 2 0 1 211 212 342 +466 2 2 0 1 178 181 219 +467 2 2 0 1 224 223 330 +468 2 2 0 1 122 123 124 +469 2 2 0 1 232 233 234 +470 2 2 0 1 19 18 73 +471 2 2 0 1 203 204 205 +472 2 2 0 1 177 162 180 +473 2 2 0 1 177 180 189 +474 2 2 0 1 193 194 203 +475 2 2 0 1 24 23 76 +476 2 2 0 1 205 206 209 +477 2 2 0 1 219 220 255 +478 2 2 0 1 188 187 191 +479 2 2 0 1 161 151 162 +480 2 2 0 1 14 13 77 +481 2 2 0 1 152 151 161 +482 2 2 0 1 150 149 151 +483 2 2 0 1 178 167 181 +484 2 2 0 1 148 147 149 +485 2 2 0 1 161 201 324 +486 2 2 0 1 149 147 156 +487 2 2 0 1 146 196 197 +488 2 2 0 1 26 25 95 +489 2 2 0 1 186 182 236 +490 2 2 0 1 12 11 94 +491 2 2 0 1 196 142 322 +492 2 2 0 1 197 198 199 +493 2 2 0 1 219 181 220 +494 2 2 0 1 161 162 177 +495 2 2 0 1 332 257 334 +496 2 2 0 1 255 220 256 +497 2 2 0 1 155 138 308 +498 2 2 0 1 122 120 123 +499 2 2 0 1 176 169 182 +500 2 2 0 1 108 106 192 +501 2 2 0 1 124 123 125 +502 2 2 0 1 125 123 142 +503 2 2 0 1 143 146 147 +504 2 2 0 1 160 140 169 +505 2 2 0 1 214 229 246 +506 2 2 0 1 218 191 243 +507 2 2 0 1 112 110 145 +508 2 2 0 1 199 200 213 +509 2 2 0 1 116 114 117 +510 2 2 0 1 186 236 237 +511 2 2 0 1 197 199 202 +512 2 2 0 1 125 142 143 +513 2 2 0 1 131 129 132 +514 2 2 0 1 28 27 81 +515 2 2 0 1 118 117 119 +516 2 2 0 1 120 119 322 +517 2 2 0 1 10 9 83 +518 2 2 0 1 143 142 146 +519 2 2 0 1 176 186 190 +520 2 2 0 1 226 258 345 +521 2 2 0 1 167 152 324 +522 2 2 0 1 139 140 153 +523 2 2 0 1 152 161 324 +524 2 2 0 1 119 117 210 +525 2 2 0 1 123 120 322 +526 2 2 0 1 176 182 186 +527 2 2 0 1 125 299 336 +528 2 2 0 1 226 209 258 +529 2 2 0 1 242 224 330 +530 2 2 0 1 209 206 332 +531 2 2 0 1 205 209 226 +532 2 2 0 1 193 203 259 +533 2 2 0 1 160 169 176 +534 2 2 0 1 181 167 324 +535 2 2 0 1 208 232 332 +536 2 2 0 1 206 208 332 +537 2 2 0 1 156 147 327 +538 2 2 0 1 223 155 308 +539 2 2 0 1 139 137 140 +540 2 2 0 1 153 140 160 +541 2 2 0 1 165 172 333 +542 2 2 0 1 232 234 257 +543 2 2 0 1 162 151 326 +544 2 2 0 1 151 149 326 +545 2 2 0 1 149 156 326 +546 2 2 0 1 147 146 327 +547 2 2 0 1 142 123 322 +548 2 2 0 1 114 112 328 +549 2 2 0 1 124 125 336 +550 2 2 0 1 251 253 254 +551 2 2 0 1 173 158 183 +552 2 2 0 1 146 197 327 +553 2 2 0 1 31 30 86 +554 2 2 0 1 197 202 327 +555 2 2 0 1 7 6 87 +556 2 2 0 1 172 193 333 +557 2 2 0 1 127 124 336 +558 2 2 0 1 232 257 332 +559 2 2 0 1 106 104 329 +560 2 2 0 1 117 114 328 +561 2 2 0 1 126 128 344 +562 2 2 0 1 137 134 331 +563 2 2 0 1 169 140 331 +564 2 2 0 1 112 145 328 +565 2 2 0 1 227 331 341 +566 2 2 0 1 258 209 334 +567 2 2 0 1 104 260 329 +568 2 2 0 1 210 117 328 +569 2 2 0 1 132 129 339 +570 2 2 0 1 133 293 338 +571 2 2 0 1 153 160 335 +572 2 2 0 1 192 106 329 +573 2 2 0 1 140 137 331 +574 2 2 0 1 129 127 339 +575 2 2 0 1 193 259 333 +576 2 2 0 1 322 119 337 +577 2 2 0 1 158 153 335 +578 2 2 0 1 191 187 244 +579 2 2 0 1 127 336 339 +580 2 2 0 1 131 132 341 +581 2 2 0 1 119 210 337 +582 2 2 0 1 196 322 337 +583 2 2 0 1 295 211 342 +584 2 2 0 1 183 158 335 +585 2 2 0 1 134 131 341 +586 2 2 0 1 5 4 100 +587 2 2 0 1 33 32 98 +588 2 2 0 1 331 134 341 +589 2 2 0 1 342 212 347 +590 2 2 0 1 187 173 340 +591 2 2 0 1 344 226 345 +592 2 2 0 1 72 101 414 +593 2 2 0 1 173 183 340 +594 2 2 0 1 205 226 343 +595 2 2 0 1 259 203 343 +596 2 2 0 1 244 187 340 +597 2 2 0 1 203 205 343 +598 2 2 0 1 128 135 344 +599 2 2 0 1 347 212 352 +600 2 2 0 1 321 145 348 +601 2 2 0 1 348 192 356 +602 2 2 0 1 145 110 348 +603 2 2 0 1 108 192 348 +604 2 2 0 1 320 210 346 +605 2 2 0 1 126 344 345 +606 2 2 0 1 314 189 399 +607 2 2 0 1 111 113 358 +608 2 2 0 1 236 352 357 +609 2 2 0 1 236 182 352 +610 2 2 0 1 247 248 353 +611 2 2 0 1 42 43 354 +612 2 2 0 1 178 219 360 +613 2 2 0 1 210 328 346 +614 2 2 0 1 257 234 349 +615 2 2 0 1 168 178 357 +616 2 2 0 1 304 306 364 +617 2 2 0 1 192 262 356 +618 2 2 0 1 60 256 408 +619 2 2 0 1 328 145 346 +620 2 2 0 1 354 246 365 +621 2 2 0 1 214 246 355 +622 2 2 0 1 251 254 353 +623 2 2 0 1 350 222 351 +624 2 2 0 1 199 213 359 +625 2 2 0 1 341 132 347 +626 2 2 0 1 110 108 348 +627 2 2 0 1 357 178 360 +628 2 2 0 1 323 156 363 +629 2 2 0 1 132 342 347 +630 2 2 0 1 227 341 347 +631 2 2 0 1 321 348 356 +632 2 2 0 1 227 347 352 +633 2 2 0 1 94 250 369 +634 2 2 0 1 327 202 363 +635 2 2 0 1 280 239 389 +636 2 2 0 1 234 247 349 +637 2 2 0 1 281 222 350 +638 2 2 0 1 323 363 366 +639 2 2 0 1 352 212 357 +640 2 2 0 1 229 321 356 +641 2 2 0 1 271 273 376 +642 2 2 0 1 182 227 352 +643 2 2 0 1 262 79 365 +644 2 2 0 1 248 251 353 +645 2 2 0 1 79 42 354 +646 2 2 0 1 212 168 357 +647 2 2 0 1 113 157 358 +648 2 2 0 1 213 214 355 +649 2 2 0 1 229 356 365 +650 2 2 0 1 202 199 359 +651 2 2 0 1 237 236 360 +652 2 2 0 1 218 243 361 +653 2 2 0 1 297 88 367 +654 2 2 0 1 236 357 360 +655 2 2 0 1 306 82 364 +656 2 2 0 1 286 287 367 +657 2 2 0 1 156 327 363 +658 2 2 0 1 356 262 365 +659 2 2 0 1 305 304 366 +660 2 2 0 1 363 305 366 +661 2 2 0 1 417 34 418 +662 2 2 0 1 96 362 411 +663 2 2 0 1 189 314 409 +664 2 2 0 1 177 189 368 +665 2 2 0 1 361 243 362 +666 2 2 0 1 246 229 365 +667 2 2 0 1 312 313 379 +668 2 2 0 1 54 55 379 +669 2 2 0 1 101 417 418 +670 2 2 0 1 285 217 380 +671 2 2 0 1 275 276 402 +672 2 2 0 1 362 243 411 +673 2 2 0 1 101 72 417 +674 2 2 0 1 268 29 415 +675 2 2 0 1 319 220 381 +676 2 2 0 1 220 181 381 +677 2 2 0 1 350 93 380 +678 2 2 0 1 230 245 382 +679 2 2 0 1 324 201 381 +680 2 2 0 1 239 231 382 +681 2 2 0 1 281 350 380 +682 2 2 0 1 97 270 383 +683 2 2 0 1 252 94 373 +684 2 2 0 1 88 286 367 +685 2 2 0 1 216 323 366 +686 2 2 0 1 238 218 384 +687 2 2 0 1 35 36 383 +688 2 2 0 1 364 82 390 +689 2 2 0 1 201 177 368 +690 2 2 0 1 20 19 370 +691 2 2 0 1 25 24 372 +692 2 2 0 1 27 26 371 +693 2 2 0 1 11 10 373 +694 2 2 0 1 9 8 374 +695 2 2 0 1 93 20 370 +696 2 2 0 1 273 275 376 +697 2 2 0 1 19 73 370 +698 2 2 0 1 8 7 375 +699 2 2 0 1 24 76 372 +700 2 2 0 1 95 25 372 +701 2 2 0 1 26 95 371 +702 2 2 0 1 94 11 373 +703 2 2 0 1 225 241 386 +704 2 2 0 1 81 27 371 +705 2 2 0 1 10 83 373 +706 2 2 0 1 243 191 385 +707 2 2 0 1 83 9 374 +708 2 2 0 1 91 317 400 +709 2 2 0 1 6 5 377 +710 2 2 0 1 32 31 378 +711 2 2 0 1 7 87 375 +712 2 2 0 1 31 86 378 +713 2 2 0 1 81 267 415 +714 2 2 0 1 87 6 377 +715 2 2 0 1 265 238 387 +716 2 2 0 1 5 100 377 +717 2 2 0 1 98 32 378 +718 2 2 0 1 221 283 391 +719 2 2 0 1 361 378 384 +720 2 2 0 1 242 267 392 +721 2 2 0 1 55 312 379 +722 2 2 0 1 267 81 392 +723 2 2 0 1 292 84 395 +724 2 2 0 1 351 222 391 +725 2 2 0 1 267 268 415 +726 2 2 0 1 244 340 395 +727 2 2 0 1 15 280 389 +728 2 2 0 1 237 297 404 +729 2 2 0 1 217 281 380 +730 2 2 0 1 176 190 393 +731 2 2 0 1 316 314 399 +732 2 2 0 1 335 160 393 +733 2 2 0 1 201 368 413 +734 2 2 0 1 91 51 396 +735 2 2 0 1 52 99 396 +736 2 2 0 1 219 255 407 +737 2 2 0 1 289 288 401 +738 2 2 0 1 60 61 397 +739 2 2 0 1 82 49 400 +740 2 2 0 1 50 91 400 +741 2 2 0 1 180 215 399 +742 2 2 0 1 374 8 375 +743 2 2 0 1 181 324 381 +744 2 2 0 1 95 372 386 +745 2 2 0 1 270 271 398 +746 2 2 0 1 245 240 410 +747 2 2 0 1 53 1 406 +748 2 2 0 1 1 54 406 +749 2 2 0 1 71 72 414 +750 2 2 0 1 319 80 408 +751 2 2 0 1 186 237 404 +752 2 2 0 1 59 60 408 +753 2 2 0 1 314 311 409 +754 2 2 0 1 231 230 382 +755 2 2 0 1 16 15 389 +756 2 2 0 1 298 237 407 +757 2 2 0 1 29 268 405 +758 2 2 0 1 17 92 394 +759 2 2 0 1 73 18 394 +760 2 2 0 1 36 97 383 +761 2 2 0 1 215 216 390 +762 2 2 0 1 351 22 403 +763 2 2 0 1 191 244 385 +764 2 2 0 1 21 93 403 +765 2 2 0 1 252 373 388 +766 2 2 0 1 373 83 388 +767 2 2 0 1 218 361 384 +768 2 2 0 1 29 28 415 +769 2 2 0 1 237 360 407 +770 2 2 0 1 2 35 412 +771 2 2 0 1 34 33 418 +772 2 2 0 1 93 350 403 +773 2 2 0 1 276 388 402 +774 2 2 0 1 4 2 412 +775 2 2 0 1 283 225 386 +776 2 2 0 1 378 86 384 +777 2 2 0 1 362 96 414 +778 2 2 0 1 86 30 405 +779 2 2 0 1 386 372 416 +780 2 2 0 1 268 265 387 +781 2 2 0 1 253 252 388 +782 2 2 0 1 92 16 389 +783 2 2 0 1 222 221 391 +784 2 2 0 1 101 362 414 +785 2 2 0 1 241 242 392 +786 2 2 0 1 271 376 398 +787 2 2 0 1 216 364 390 +788 2 2 0 1 74 255 397 +789 2 2 0 1 160 176 393 +790 2 2 0 1 18 17 394 +791 2 2 0 1 291 292 395 +792 2 2 0 1 350 351 403 +793 2 2 0 1 368 189 409 +794 2 2 0 1 35 383 412 +795 2 2 0 1 51 52 396 +796 2 2 0 1 255 256 397 +797 2 2 0 1 61 74 397 +798 2 2 0 1 189 180 399 +799 2 2 0 1 49 50 400 +800 2 2 0 1 288 291 401 +801 2 2 0 1 283 386 416 +802 2 2 0 1 376 87 398 +803 2 2 0 1 297 367 404 +804 2 2 0 1 81 371 392 +805 2 2 0 1 391 283 416 +806 2 2 0 1 318 319 413 +807 2 2 0 1 387 86 405 +808 2 2 0 1 72 3 417 +809 2 2 0 1 310 368 409 +810 2 2 0 1 54 379 406 +811 2 2 0 1 30 29 405 +812 2 2 0 1 3 34 417 +813 2 2 0 1 22 21 403 +814 2 2 0 1 190 186 404 +815 2 2 0 1 360 219 407 +816 2 2 0 1 99 53 406 +817 2 2 0 1 80 59 408 +818 2 2 0 1 278 96 411 +819 2 2 0 1 240 285 410 +820 2 2 0 1 311 310 409 +821 2 2 0 1 87 377 398 +822 2 2 0 1 367 190 404 +823 2 2 0 1 243 385 411 +824 2 2 0 1 96 71 414 +825 2 2 0 1 381 201 413 +826 2 2 0 1 100 4 412 +827 2 2 0 1 28 81 415 +828 2 2 0 1 383 100 412 +829 2 2 0 1 388 83 402 +830 2 2 0 1 33 98 418 +831 2 2 0 1 379 99 406 +832 2 2 0 1 83 374 402 +833 2 2 0 1 319 381 413 +834 2 2 0 1 372 76 416 +835 2 2 0 1 268 387 405 +836 2 2 0 1 76 391 416 +837 2 2 0 1 385 278 411 +$EndElements diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/main.cc b/test/porousmediumflow/2p/implicit/rotationsymmetry/main.cc new file mode 100644 index 0000000000000000000000000000000000000000..f0880978429829af9df3c8351eaceb79f9945576 --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/main.cc @@ -0,0 +1,155 @@ +// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- +// vi: set et ts=4 sw=4 sts=4: +/***************************************************************************** + * See the file COPYING for full copying permissions. * + * * + * This program is free software: you can redistribute it and/or modify * + * it under the terms of the GNU General Public License as published by * + * the Free Software Foundation, either version 3 of the License, or * + * (at your option) any later version. * + * * + * This program is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * + * GNU General Public License for more details. * + * * + * You should have received a copy of the GNU General Public License * + * along with this program. If not, see <http://www.gnu.org/licenses/>. * + *****************************************************************************/ + +#include <config.h> + +#include <iostream> +#include <dune/common/parallel/mpihelper.hh> + +#include <dumux/common/properties.hh> +#include <dumux/common/parameters.hh> + +#include <dumux/linear/seqsolverbackend.hh> +#include <dumux/linear/linearsolvertraits.hh> +#include <dumux/nonlinear/newtonsolver.hh> +#include <dumux/assembly/fvassembler.hh> + +#include <dumux/io/vtkoutputmodule.hh> +#include <dumux/io/grid/gridmanager_alu.hh> + +#include "properties.hh" + +int main(int argc, char** argv) try +{ + using namespace Dumux; + + // initialize MPI, finalize is done automatically on exit + const auto& mpiHelper = Dune::MPIHelper::instance(argc, argv); + + // parse command line arguments and input file + Parameters::init(argc, argv); + + // the type tag (see properties.hh) + using TypeTag = Properties::TTag::TwoPRotationalSymmetryDome; + + // try to create a grid (from the given grid file or the input file) + GridManager<GetPropType<TypeTag, Properties::Grid>> gridManager; + gridManager.init(); + + // we compute on the leaf grid view + const auto& leafGridView = gridManager.grid().leafGridView(); + + // create the grid geometry + using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; + auto gridGeometry = std::make_shared<GridGeometry>(leafGridView); + gridGeometry->update(); + + // the problem (initial and boundary conditions) + using Problem = GetPropType<TypeTag, Properties::Problem>; + auto problem = std::make_shared<Problem>(gridGeometry); + + // the solution vector + using SolutionVector = GetPropType<TypeTag, Properties::SolutionVector>; + SolutionVector x; + problem->applyInitialSolution(x); + auto xOld = x; + + // the grid variables + using GridVariables = GetPropType<TypeTag, Properties::GridVariables>; + auto gridVariables = std::make_shared<GridVariables>(problem, gridGeometry); + gridVariables->init(x); + + // use non-conforming output for the test with interface solver + VtkOutputModule<GridVariables, SolutionVector> vtkWriter(*gridVariables, x, problem->name()); + GetPropType<TypeTag, Properties::IOFields>::initOutputModule(vtkWriter); // Add model specific output fields + vtkWriter.write(0.0); + + // get some time loop parameters + using Scalar = GetPropType<TypeTag, Properties::Scalar>; + const auto tEnd = getParam<Scalar>("TimeLoop.TEnd"); + const auto maxDt = getParam<Scalar>("TimeLoop.MaxTimeStepSize"); + const auto dt = getParam<Scalar>("TimeLoop.DtInitial"); + + // instantiate time loop + auto timeLoop = std::make_shared<TimeLoop<Scalar>>(0.0, dt, tEnd); + timeLoop->setMaxTimeStepSize(maxDt); + + // the assembler with time loop for instationary problem + using Assembler = FVAssembler<TypeTag, DiffMethod::numeric>; + auto assembler = std::make_shared<Assembler>(problem, gridGeometry, gridVariables, timeLoop, xOld); + + // the linear solver + using LinearSolver = ILU0RestartedGMResBackend; + auto linearSolver = std::make_shared<LinearSolver>(); + + // the non-linear solver + using NewtonSolver = Dumux::NewtonSolver<Assembler, LinearSolver>; + NewtonSolver nonLinearSolver(assembler, linearSolver); + + // time loop + const int outputInterval = getParam<int>("Vtk.OutputInterval"); + timeLoop->start(); do + { + // solve the non-linear system with time step control + nonLinearSolver.solve(x, *timeLoop); + + // make the new solution the old solution + xOld = x; + gridVariables->advanceTimeStep(); + + // advance to the time loop to the next step + timeLoop->advanceTimeStep(); + + // write vtk output every timeStepInterval time step + if (timeLoop->timeStepIndex() % outputInterval == 0) + vtkWriter.write(timeLoop->time()); + + // report statistics of this time step + timeLoop->reportTimeStep(); + + // set new dt as suggested by the Newton solver + timeLoop->setTimeStepSize(nonLinearSolver.suggestTimeStepSize(timeLoop->timeStepSize())); + + } while (!timeLoop->finished()); + + // output some Newton statistics + nonLinearSolver.report(); + + timeLoop->finalize(leafGridView.comm()); + + //////////////////////////////////////////////////////////// + // finalize, print dumux message to say goodbye + //////////////////////////////////////////////////////////// + + // print dumux end message + if (mpiHelper.rank() == 0) + Parameters::print(); + + return 0; +} // end main +catch (const Dumux::ParameterException &e) +{ + std::cerr << std::endl << e << " ---> Abort!" << std::endl; + return 1; +} +catch (const Dune::Exception &e) +{ + std::cerr << "Dune reported error: " << e << " ---> Abort!" << std::endl; + return 2; +} diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/params.input b/test/porousmediumflow/2p/implicit/rotationsymmetry/params.input new file mode 100644 index 0000000000000000000000000000000000000000..efe4baef169dea37dfcdef341136fb6ab0dd6fa6 --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/params.input @@ -0,0 +1,14 @@ +[TimeLoop] +DtInitial = 1e2 # [s] +TEnd = 1e7 # [s] +MaxTimeStepSize = 1e4 + +[Grid] +File = dome.msh + +[Problem] +Name = test_2p_rotationsymmetry_dome +EnableGravity = true + +[Vtk] +OutputInterval = 50 diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/problem.hh b/test/porousmediumflow/2p/implicit/rotationsymmetry/problem.hh new file mode 100644 index 0000000000000000000000000000000000000000..d742db3f709a908ebb4a5afb8b8ab152a32911d8 --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/problem.hh @@ -0,0 +1,86 @@ +// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- +// vi: set et ts=4 sw=4 sts=4: +/***************************************************************************** + * See the file COPYING for full copying permissions. * + * * + * This program is free software: you can redistribute it and/or modify * + * it under the terms of the GNU General Public License as published by * + * the Free Software Foundation, either version 3 of the License, or * + * (at your option) any later version. * + * * + * This program is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * + * GNU General Public License for more details. * + * * + * You should have received a copy of the GNU General Public License * + * along with this program. If not, see <http://www.gnu.org/licenses/>. * + *****************************************************************************/ + +#ifndef DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_PROBLEM_HH +#define DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_PROBLEM_HH + +#include <dumux/common/properties.hh> +#include <dumux/common/boundarytypes.hh> +#include <dumux/porousmediumflow/problem.hh> + +namespace Dumux { + +/*! + * \brief A rotational symmetric 2p problem: water and air separate in a density-driven process in a dome shaped domain + */ +template<class TypeTag> +class TwoPRotationalSymmetryProblem : public PorousMediumFlowProblem<TypeTag> +{ + using ParentType = PorousMediumFlowProblem<TypeTag>; + using GridView = typename GetPropType<TypeTag, Properties::GridGeometry>::GridView; + using Element = typename GridView::template Codim<0>::Entity; + using Scalar = GetPropType<TypeTag, Properties::Scalar>; + using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>; + using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; + using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>; + using BoundaryTypes = Dumux::BoundaryTypes<PrimaryVariables::size()>; + using GlobalPosition = typename Element::Geometry::GlobalCoordinate; + +public: + TwoPRotationalSymmetryProblem(std::shared_ptr<const GridGeometry> gridGeometry) + : ParentType(gridGeometry) + { + FluidSystem::init(/*tempMin=*/293.0, /*tempMax=*/294.0, /*numTemp=*/2, + /*pMin=*/1.0e4, /*pMax=*/1.0e6, /*numP=*/200); + } + + /*! + * \brief The boundary types at position globalPos + */ + BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const + { + BoundaryTypes values; + values.setAllNeumann(); + return values; + } + + /*! + * \brief The initial values at position globalPos + */ + PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const + { + PrimaryVariables values; + const Scalar depth = this->gridGeometry().bBoxMax()[1] - globalPos[1]; + // hydrostatic pressure + values[0] = 1e5 - 1000*this->spatialParams().gravity(globalPos)[1]*depth; + // start with saturation 0.5 -> density driven demixing + values[1] = 0.5; + return values; + } + + /*! + * \brief The temperature \f$\mathrm{[K]}\f$ in the domain + */ + Scalar temperature() const + { return 293.15; } +}; + +} // end namespace Dumux + +#endif diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/properties.hh b/test/porousmediumflow/2p/implicit/rotationsymmetry/properties.hh new file mode 100644 index 0000000000000000000000000000000000000000..f6eb95b78efb9268b296ab2b36be535eb7e2798b --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/properties.hh @@ -0,0 +1,99 @@ +// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- +// vi: set et ts=4 sw=4 sts=4: +/***************************************************************************** + * See the file COPYING for full copying permissions. * + * * + * This program is free software: you can redistribute it and/or modify * + * it under the terms of the GNU General Public License as published by * + * the Free Software Foundation, either version 3 of the License, or * + * (at your option) any later version. * + * * + * This program is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * + * GNU General Public License for more details. * + * * + * You should have received a copy of the GNU General Public License * + * along with this program. If not, see <http://www.gnu.org/licenses/>. * + *****************************************************************************/ + +#ifndef DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_PROPERTIES_HH +#define DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_PROPERTIES_HH + +#include <dune/alugrid/grid.hh> + +#include <dumux/common/properties.hh> + +#include <dumux/discretization/box.hh> +#include <dumux/discretization/rotationsymmetricgridgeometrytraits.hh> + +#include <dumux/porousmediumflow/2p/model.hh> + +#include <dumux/material/components/tabulatedcomponent.hh> +#include <dumux/material/components/h2o.hh> +#include <dumux/material/components/air.hh> +#include <dumux/material/fluidsystems/1pliquid.hh> +#include <dumux/material/fluidsystems/1pgas.hh> +#include <dumux/material/fluidsystems/2pimmiscible.hh> + +#include "problem.hh" +#include "spatialparams.hh" + +namespace Dumux::Properties { + +namespace TTag { +struct TwoPRotationalSymmetryDome +{ using InheritsFrom = std::tuple<TwoP, BoxModel>; }; +} // end namespace TTag + +// Set the grid type +template<class TypeTag> +struct Grid<TypeTag, TTag::TwoPRotationalSymmetryDome> +{ using type = Dune::ALUGrid<2, 2, Dune::simplex, Dune::conforming>; }; + +// Set the problem type +template<class TypeTag> +struct Problem<TypeTag, TTag::TwoPRotationalSymmetryDome> +{ using type = TwoPRotationalSymmetryProblem<TypeTag>; }; + +// Set the fluid system +template<class TypeTag> +struct FluidSystem<TypeTag, TTag::TwoPRotationalSymmetryDome> +{ + using Scalar = GetPropType<TypeTag, Properties::Scalar>; + using WettingPhase = FluidSystems::OnePLiquid<Scalar, Components::TabulatedComponent<Components::H2O<Scalar>, false>>; + using NonwettingPhase = FluidSystems::OnePGas<Scalar, Components::TabulatedComponent<Components::Air<Scalar>, false>>; + using type = FluidSystems::TwoPImmiscible<Scalar, WettingPhase, NonwettingPhase>; +}; + +// Set the spatial parameters +template<class TypeTag> +struct SpatialParams<TypeTag, TTag::TwoPRotationalSymmetryDome> +{ + using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; + using Scalar = GetPropType<TypeTag, Properties::Scalar>; + using type = TwoPRotationalSymmetrySpatialParams<GridGeometry, Scalar>; +}; + +// Set the rotational symmetric grid geometry +template<class TypeTag> +struct GridGeometry<TypeTag, TTag::TwoPRotationalSymmetryDome> +{ + static constexpr bool enableCache = getPropValue<TypeTag, Properties::EnableGridGeometryCache>(); + using GridView = typename GetPropType<TypeTag, Properties::Grid>::LeafGridView; + using Scalar = GetPropType<TypeTag, Properties::Scalar>; + using Traits = RotationSymmetricGridGeometryTraits<BoxDefaultGridGeometryTraits<GridView>, RotationPolicy::toroid>; + using type = BoxFVGridGeometry<Scalar, GridView, enableCache, Traits>; +}; + +// Caching options +template<class TypeTag> +struct EnableGridVolumeVariablesCache<TypeTag, TTag::TwoPRotationalSymmetryDome> { static constexpr bool value = true; }; +template<class TypeTag> +struct EnableGridFluxVariablesCache<TypeTag, TTag::TwoPRotationalSymmetryDome> { static constexpr bool value = true; }; +template<class TypeTag> +struct EnableGridGeometryCache<TypeTag, TTag::TwoPRotationalSymmetryDome> { static constexpr bool value = true; }; + +} // end namespace Dumux::Properties + +#endif diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/solution.png b/test/porousmediumflow/2p/implicit/rotationsymmetry/solution.png new file mode 100644 index 0000000000000000000000000000000000000000..979feb3239a6cc08f15da2b4f13ef24917fba960 Binary files /dev/null and b/test/porousmediumflow/2p/implicit/rotationsymmetry/solution.png differ diff --git a/test/porousmediumflow/2p/implicit/rotationsymmetry/spatialparams.hh b/test/porousmediumflow/2p/implicit/rotationsymmetry/spatialparams.hh new file mode 100644 index 0000000000000000000000000000000000000000..da0361a7f85e0c4d6342bb27655bfbc0b8e3c109 --- /dev/null +++ b/test/porousmediumflow/2p/implicit/rotationsymmetry/spatialparams.hh @@ -0,0 +1,77 @@ +// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- +// vi: set et ts=4 sw=4 sts=4: +/***************************************************************************** + * See the file COPYING for full copying permissions. * + * * + * This program is free software: you can redistribute it and/or modify * + * it under the terms of the GNU General Public License as published by * + * the Free Software Foundation, either version 3 of the License, or * + * (at your option) any later version. * + * * + * This program is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * + * GNU General Public License for more details. * + * * + * You should have received a copy of the GNU General Public License * + * along with this program. If not, see <http://www.gnu.org/licenses/>. * + *****************************************************************************/ +/*! + * \file + * \brief The spatial params 2p rotational symmetry test + */ +#ifndef DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_SPATIAL_PARAMS_HH +#define DUMUX_TEST_TWOP_ROTATIONALSYMMETRY_SPATIAL_PARAMS_HH + +#include <dumux/material/spatialparams/fv.hh> +#include <dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh> +#include <dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh> + +namespace Dumux { + +template<class GridGeometry, class Scalar> +class TwoPRotationalSymmetrySpatialParams +: public FVSpatialParams<GridGeometry, Scalar, TwoPRotationalSymmetrySpatialParams<GridGeometry, Scalar>> +{ + using ThisType = TwoPRotationalSymmetrySpatialParams<GridGeometry, Scalar>; + using ParentType = FVSpatialParams<GridGeometry, Scalar, ThisType>; + + using Element = typename GridGeometry::GridView::template Codim<0>::Entity; + using GlobalPosition = typename Element::Geometry::GlobalCoordinate; +public: + using MaterialLaw = EffToAbsLaw<RegularizedBrooksCorey<Scalar>>; + using MaterialLawParams = typename MaterialLaw::Params; + using PermeabilityType = Scalar; + + TwoPRotationalSymmetrySpatialParams(std::shared_ptr<const GridGeometry> gridGeometry) + : ParentType(gridGeometry) + { + // residual saturations + materialLawParams_.setSwr(0.0); + materialLawParams_.setSnr(0.0); + + // parameters for the material law + materialLawParams_.setPe(0.0); + materialLawParams_.setLambda(2.0); + } + + PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const + { return 1e-11; } + + Scalar porosityAtPos(const GlobalPosition& globalPos) const + { return 0.4; } + + const MaterialLawParams& materialLawParamsAtPos(const GlobalPosition& globalPos) const + { return materialLawParams_; } + + template<class FluidSystem> + int wettingPhaseAtPos(const GlobalPosition& globalPos) const + { return FluidSystem::phase0Idx; } + +private: + MaterialLawParams materialLawParams_; +}; + +} // end namespace Dumux + +#endif diff --git a/test/references/test_2p_rotationsymmetry_dome-reference.vtu b/test/references/test_2p_rotationsymmetry_dome-reference.vtu new file mode 100644 index 0000000000000000000000000000000000000000..5d87141da8599f421ea5d5e1004ff9555ff99044 --- /dev/null +++ b/test/references/test_2p_rotationsymmetry_dome-reference.vtu @@ -0,0 +1,883 @@ +<?xml version="1.0"?> +<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian"> + <UnstructuredGrid> + <Piece NumberOfCells="762" NumberOfPoints="418"> + <PointData Scalars="S_liq"> + <DataArray type="Float32" Name="S_liq" NumberOfComponents="1" format="ascii"> + 0.993144 0.987879 0.98962 0.993821 0.989028 0.993496 0.985678 0.985262 0.984165 0.981629 0.978844 0.960212 + 0.981236 0.957571 0.981299 0.957596 0.977927 0.983448 0.985641 0.988836 0.988749 0.993395 0.99338 0.988808 + 0.985671 0.985662 0.988954 0.993467 0.993959 0.99285 0.989443 0.985667 0.987453 0.983843 0.983292 0.981525 + 0.976423 0.948351 0.977192 0.981976 0.955102 0.559271 0.44786 0.431788 0.428104 0.540599 0.441041 0.512866 + 0.423183 0.419487 0.421547 0.414289 0.412521 0.404787 0.410284 0.416999 0.425777 0.423073 0.414116 0.420045 + 0.428508 0.425222 0.416518 0.432759 0.471676 0.493913 0.43126 0.470397 0.768523 0.435025 0.426823 0.421965 + 0.417387 0.412457 0.407807 0.402611 0.406961 0.410817 0.404488 0.40088 0.396786 0.39021 0.394102 0.382953 + 0.385778 0.392082 0.397678 0.386686 0.374917 0.380644 0.37386 0.368414 0.361108 0.36231 0.35527 0.348575 + 0.340809 0.347282 0.352916 0.366303 0.338278 0.343781 0.357953 0.370952 0.378698 0.375018 0.386674 0.378497 + 0.36621 0.362406 0.352857 0.348666 0.333572 0.338232 0.342077 0.356337 0.326347 0.329813 0.345186 0.359119 + 0.369372 0.381381 0.37189 0.383672 0.392505 0.394602 0.389863 0.40043 0.397514 0.407595 0.40286 0.406207 + 0.396164 0.385389 0.386539 0.397185 0.387094 0.387031 0.397642 0.375717 0.375727 0.375079 0.363437 0.363466 + 0.362678 0.37379 0.36124 0.349194 0.347555 0.332498 0.315694 0.317866 0.334371 0.319069 0.300715 0.300915 + 0.335375 0.319249 0.350053 0.335494 0.350137 0.334758 0.3494 0.332972 0.318489 0.316585 0.300061 0.27846 + 0.277666 0.298209 0.272838 0.252884 0.233786 0.29442 0.265906 0.289194 0.313318 0.330084 0.308903 0.345162 + 0.326346 0.347742 0.361203 0.362752 0.375002 0.38632 0.373539 0.358835 0.341801 0.355643 0.371321 0.368198 + 0.382794 0.384926 0.39543 0.396826 0.394646 0.379756 0.410527 0.849754 0.430654 0.966315 0.750728 0.406827 + 0.453694 0.974299 0.904973 0.407318 0.420088 0.40747 0.397528 0.407103 0.416368 0.415572 0.43204 0.936275 + 0.482404 0.975679 0.981283 0.983 0.985558 0.985653 0.983233 0.985723 0.981369 0.988945 0.989483 0.993813 + 0.993353 0.988807 0.988784 0.993361 0.993375 0.988738 0.993362 0.985499 0.984999 0.980415 0.985332 0.980595 + 0.980767 0.988697 0.988779 0.993362 0.993358 0.98881 0.985404 0.985505 0.987737 0.992961 0.993044 0.989107 + 0.99368 0.985297 0.986768 0.983128 0.979229 0.983219 0.961547 0.48413 0.966137 0.979662 0.543506 0.938855 + 0.375582 0.387005 0.363746 0.350992 0.344381 0.357478 0.335098 0.329482 0.320982 0.349095 0.3233 0.350916 + 0.385585 0.310728 0.294615 0.303117 0.293682 0.261159 0.288408 0.2623 0.208593 0.278143 0.239516 0.234186 + 0.271694 0.926752 0.970966 0.982027 0.982952 0.986122 0.985724 0.977954 0.990027 0.989608 0.99416 0.99357 + 0.988604 0.989938 0.987427 0.993532 0.203421 0.258077 0.194078 0.255224 0.299455 0.285888 0.244305 0.186406 + 0.177759 0.240844 0.27155 0.276326 0.298077 0.303914 0.314449 0.321851 0.337224 0.282799 0.260179 0.244681 + 0.177021 0.17138 0.222474 0.160526 0.228947 0.169786 0.218252 0.164513 0.150479 0.157291 0.151361 0.217751 + 0.207864 0.249702 0.254329 0.255023 0.220265 0.280188 0.279456 0.253282 0.219503 0.158704 0.158388 0.21828 + 0.25092 0.215438 0.277373 0.274463 0.299148 0.296642 0.312679 0.293204 0.270341 0.247293 0.242252 0.211027 + 0.206223 0.150987 0.154317 0.15766 0.148362 0.202048 0.147068 0.2043 0.147656 0.207037 0.148023 0.2015 + 0.146006 0.237699 0.230917 0.24437 0.266509 0.275191 0.258315 0.236486 0.264937 0.288619 0.282945 0.303832 + 0.30876 0.322035 0.316745 0.328151 0.310576 0.297622 0.290404 0.282406 0.303591 0.295271 0.272238 0.314828 + 0.321988 0.331973 0.324696 0.333316 0.316304 0.297059 0.306567 0.285728 0.274893 0.261112 0.248962 0.23126 + 0.24102 0.25744 0.229309 0.202347 0.201865 0.14543 0.140641 0.193769 0.214662 0.149263 + </DataArray> + <DataArray type="Float32" Name="p_liq" NumberOfComponents="1" format="ascii"> + 905373 856540 869563 905308 865180 905084 827718 822588 807976 782351 759748 735880 + 780641 734778 776979 733167 752925 792061 822689 863229 862149 904695 904127 861688 + 819832 821013 861833 903387 902495 901453 866745 820183 844832 799447 788205 773291 + 750785 728394 750593 779221 730418 715402 714737 714463 714508 714972 714761 715031 + 714545 714052 714297 714086 713843 713633 713603 713811 714263 714018 713581 713796 + 714198 714044 713592 714425 714352 715223 714424 715096 716491 714659 714096 713845 + 713607 713376 713154 712934 713150 713360 713143 712919 712714 712484 712703 712270 + 712271 712492 712714 712280 711850 712057 711839 711631 711417 711423 711212 711000 + 710789 710996 711206 711626 710788 710999 711417 711842 712051 712064 712494 712288 + 711858 711637 711433 711214 710796 711012 711231 711654 710811 711030 711451 711878 + 712082 712514 712308 712744 712952 713186 712721 713160 712929 713369 713394 713876 + 713423 712978 713214 713664 713454 713696 713910 713242 713005 712771 712565 712798 + 712335 712538 712105 711903 711677 711252 710833 711055 711478 711278 710858 711078 + 711703 711501 712132 711928 712360 712156 712591 712385 711726 711952 711300 710883 + 711119 711528 711336 710996 711085 711754 711568 711982 712180 712615 712408 713056 + 712845 712823 713268 713033 713482 713942 713722 713504 713289 713740 713961 714202 + 714426 714186 714664 714415 714880 714674 715167 719761 715178 734179 716467 714912 + 715140 742589 722005 714636 714843 714374 714160 714122 714586 714334 714797 726015 + 715066 746214 771549 782448 815676 817666 785256 819971 771915 857404 864771 900341 + 899160 856268 855372 897986 896890 854370 895935 814086 811095 767116 811978 767914 + 775676 853075 852676 895158 894557 851763 812049 812762 842905 894089 893705 852146 + 893379 807760 822607 788421 762861 784338 739496 715514 739842 768694 715692 733454 + 714942 715189 714446 713976 714210 714683 714435 713752 713979 714903 714640 715062 + 715395 714829 715010 714171 713767 714125 714447 713839 713940 714688 714790 714355 + 715430 731556 748612 776670 796669 813357 813732 761912 858143 854869 893118 892953 + 841289 863686 844779 892897 713531 713468 713131 713102 713493 713023 712718 712740 + 712359 712359 712672 712441 712865 712636 713305 713076 713521 712203 711791 712003 + 711991 711635 711618 711292 711376 710964 710691 710650 710352 709806 710071 710232 + 710412 710692 710490 710233 710016 710664 710440 710015 709779 709558 709329 709579 + 709799 709363 710221 710004 710637 710417 710614 710201 709789 709590 709379 709168 + 708996 708752 708925 709118 708599 708858 708465 708728 708351 708656 708258 708567 + 708185 708929 709046 709084 709343 709551 709359 709202 709581 709984 709769 710180 + 710397 710595 710379 710582 710167 709966 709752 709542 709956 709746 709334 710161 + 710371 710578 710369 710579 710161 709745 709952 709539 709330 709130 708917 708723 + 708870 709129 708796 708467 708488 708133 708102 708385 708504 708092 + </DataArray> + <DataArray type="Float32" Name="rho_liq" NumberOfComponents="1" format="ascii"> + 998.573 998.551 998.557 998.573 998.555 998.573 998.538 998.535 998.529 998.517 998.507 998.496 + 998.516 998.495 998.515 998.495 998.504 998.521 998.535 998.554 998.553 998.573 998.573 998.553 + 998.534 998.535 998.553 998.572 998.572 998.571 998.556 998.534 998.546 998.525 998.52 998.513 + 998.503 998.492 998.502 998.516 998.493 998.486 998.486 998.486 998.486 998.486 998.486 998.486 + 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 + 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.487 998.486 998.486 998.486 + 998.486 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 + 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.484 998.484 + 998.484 998.484 998.484 998.485 998.484 998.484 998.485 998.485 998.485 998.485 998.485 998.485 + 998.485 998.485 998.485 998.484 998.484 998.484 998.484 998.485 998.484 998.484 998.485 998.485 + 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.486 + 998.485 998.485 998.485 998.486 998.485 998.486 998.486 998.485 998.485 998.485 998.485 998.485 + 998.485 998.485 998.485 998.485 998.485 998.484 998.484 998.484 998.485 998.484 998.484 998.484 + 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.484 998.484 + 998.484 998.485 998.485 998.484 998.484 998.485 998.485 998.485 998.485 998.485 998.485 998.485 + 998.485 998.485 998.485 998.485 998.486 998.486 998.486 998.486 998.485 998.486 998.486 998.486 + 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.488 998.486 998.495 998.487 998.486 + 998.486 998.499 998.489 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.491 + 998.486 998.5 998.512 998.517 998.532 998.533 998.518 998.534 998.512 998.551 998.555 998.571 + 998.57 998.551 998.55 998.57 998.569 998.55 998.569 998.531 998.53 998.51 998.531 998.51 + 998.514 998.549 998.549 998.568 998.568 998.549 998.531 998.531 998.545 998.568 998.568 998.549 + 998.568 998.529 998.535 998.52 998.508 998.518 998.497 998.486 998.498 998.511 998.487 998.495 + 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 + 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 998.486 + 998.486 998.494 998.502 998.514 998.523 998.531 998.531 998.508 998.552 998.55 998.568 998.568 + 998.544 998.554 998.545 998.567 998.486 998.485 998.485 998.485 998.486 998.485 998.485 998.485 + 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.485 998.486 998.485 998.485 998.485 + 998.485 998.485 998.485 998.484 998.485 998.484 998.484 998.484 998.484 998.484 998.484 998.484 + 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 + 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 + 998.483 998.483 998.483 998.484 998.483 998.483 998.483 998.483 998.483 998.483 998.483 998.483 + 998.483 998.483 998.483 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 + 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 + 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.484 998.483 998.483 + 998.483 998.484 998.483 998.483 998.483 998.483 998.483 998.483 998.483 998.483 + </DataArray> + <DataArray type="Float32" Name="mob_liq" NumberOfComponents="1" format="ascii"> + 971.486 951.033 957.759 974.136 955.469 972.864 942.579 940.988 936.798 927.172 916.688 848.855 + 925.689 839.554 925.925 839.64 913.257 934.068 942.435 954.725 954.389 972.468 972.41 954.617 + 942.547 942.514 955.182 972.75 974.677 970.335 957.072 942.534 949.391 935.571 933.472 926.777 + 907.651 807.68 910.513 928.482 830.928 97.69 40.1728 34.7093 33.5398 85.2832 37.7815 69.0841 + 32.024 30.9197 31.5317 29.4154 28.9165 26.8081 28.2943 30.1926 32.8163 31.9907 29.3662 31.0846 + 33.6664 32.6455 30.0537 35.0226 49.4241 59.4245 34.5398 48.8902 348.329 35.7618 33.1402 31.6568 + 30.3051 28.8986 27.6174 26.2363 27.3889 28.4417 26.7292 25.7882 24.7506 23.1503 24.0878 21.4754 + 22.1163 23.5976 24.9741 22.3252 19.7288 20.9622 19.5072 18.3953 16.9789 17.206 15.9073 14.7416 + 13.4712 14.5241 15.4899 17.9773 13.0755 13.9473 16.3933 18.9074 20.5369 19.7502 22.3224 20.4934 + 17.9591 17.2244 15.4795 14.7571 12.3629 13.0683 13.6729 16.0993 11.3262 11.815 14.1766 16.608 + 18.5875 21.1251 19.0995 21.6373 23.6998 24.2103 23.068 25.6725 24.9328 27.5599 26.3013 27.1864 + 24.596 22.0271 22.2913 24.8505 22.4197 22.4051 24.965 19.8979 19.8999 19.763 17.4213 17.4268 + 17.2762 19.4927 17.0038 14.8466 14.5699 12.2045 9.9181 10.1938 12.4818 10.3491 8.16554 8.18728 + 12.6324 10.3724 14.9933 12.6503 15.0077 12.5396 14.8817 12.2742 10.2739 10.0306 8.09465 6.00365 + 5.93544 7.89672 5.53327 4.08365 2.9829 7.50289 4.99202 6.98425 9.6228 11.8539 9.09183 14.1727 + 11.326 14.6013 16.9969 17.2903 19.7467 22.2409 19.4405 16.5555 13.6286 15.9742 18.9827 18.3521 + 21.4398 21.9214 24.414 24.7607 24.2211 20.7674 28.3614 520.637 34.346 870.643 317.171 27.3527 + 42.3073 899.776 669.739 27.4849 31.0975 27.5262 24.9364 27.4271 30.0104 29.7817 34.7903 767.323 + 54.0762 904.888 925.864 932.364 942.115 942.481 933.25 942.75 926.19 955.144 957.227 974.103 + 972.3 954.613 954.522 972.333 972.387 954.346 972.337 941.888 939.98 922.591 941.252 923.267 + 923.92 954.185 954.503 972.335 972.319 954.623 941.525 941.912 950.482 970.765 971.089 955.769 + 973.581 941.116 946.753 932.852 918.134 933.196 853.586 54.8544 870.004 919.761 87.1325 775.814 + 19.8692 22.399 17.4805 15.1548 14.0449 16.3064 12.5907 11.7677 10.5994 14.8299 10.9089 15.1417 + 22.0722 9.30858 7.52286 8.42948 7.428 4.64495 6.90865 4.72667 1.89043 5.97632 3.28626 3.00334 + 5.44107 736.577 887.528 928.674 932.186 944.274 942.752 913.362 959.333 957.71 975.465 973.149 + 953.826 958.989 949.291 973 1.70979 4.42955 1.41666 4.23688 8.02954 6.67029 3.55703 1.20559 + 0.996988 3.35972 5.42949 5.82169 7.88269 8.51854 9.7625 10.7147 12.9132 6.38668 4.57562 3.57899 + 0.980528 0.861386 2.44612 0.66304 2.74349 0.829783 2.26565 0.731412 0.511994 0.611187 0.5241 2.24494 + 1.86414 3.88197 4.17775 4.22358 2.3504 6.15405 6.08995 4.10943 2.31806 0.633452 0.628418 2.26681 + 3.95827 2.15105 5.91045 5.66626 7.99661 7.73198 9.5446 7.37975 5.33347 3.73431 3.439 1.98022 + 1.80598 0.518944 0.566259 0.616945 0.483783 1.66411 0.467124 1.73953 0.474635 1.83465 0.479373 1.64612 + 0.453784 3.18766 2.83912 3.56082 5.03743 5.72661 4.4459 3.12307 4.91964 6.92884 6.39986 8.50934 + 9.07498 10.7393 10.0508 11.5787 9.29042 7.83465 7.10186 6.3512 8.48235 7.59005 5.48472 9.80964 + 10.733 12.1276 11.0987 12.325 9.99489 7.77555 8.81989 6.65539 5.70187 4.64159 3.83613 2.85604 + 3.36954 4.386 2.76088 1.67397 1.65807 0.446656 0.390672 1.40767 2.12023 0.495644 + </DataArray> + <DataArray type="Float32" Name="S_gas" NumberOfComponents="1" format="ascii"> + 0.00685594 0.0121212 0.0103801 0.00617915 0.0109718 0.00650372 0.0143218 0.0147376 0.015835 0.0183711 0.0211562 0.039788 + 0.0187639 0.042429 0.018701 0.0424043 0.022073 0.0165518 0.014359 0.0111644 0.0112512 0.00660494 0.00661966 0.0111921 + 0.0143295 0.0143383 0.011046 0.00653271 0.00604113 0.00714993 0.0105575 0.0143329 0.0125469 0.0161568 0.0167084 0.0184751 + 0.023577 0.0516494 0.0228083 0.0180245 0.0448978 0.440729 0.55214 0.568212 0.571896 0.459401 0.558959 0.487134 + 0.576817 0.580513 0.578453 0.585711 0.587479 0.595213 0.589716 0.583001 0.574223 0.576927 0.585884 0.579955 + 0.571492 0.574778 0.583482 0.567241 0.528324 0.506087 0.56874 0.529603 0.231477 0.564975 0.573177 0.578035 + 0.582613 0.587543 0.592193 0.597389 0.593039 0.589183 0.595512 0.59912 0.603214 0.60979 0.605898 0.617047 + 0.614222 0.607918 0.602322 0.613314 0.625083 0.619356 0.62614 0.631586 0.638892 0.63769 0.64473 0.651425 + 0.659191 0.652718 0.647084 0.633697 0.661722 0.656219 0.642047 0.629048 0.621302 0.624982 0.613326 0.621503 + 0.633789 0.637594 0.647143 0.651334 0.666428 0.661768 0.657923 0.643663 0.673653 0.670187 0.654814 0.640881 + 0.630628 0.618619 0.62811 0.616328 0.607495 0.605398 0.610137 0.59957 0.602486 0.592405 0.59714 0.593793 + 0.603836 0.614611 0.613461 0.602815 0.612906 0.612969 0.602358 0.624283 0.624273 0.624921 0.636563 0.636534 + 0.637322 0.62621 0.63876 0.650806 0.652445 0.667502 0.684306 0.682134 0.665629 0.680931 0.699284 0.699085 + 0.664625 0.680751 0.649947 0.664506 0.649863 0.665242 0.6506 0.667028 0.681512 0.683415 0.699939 0.72154 + 0.722334 0.701791 0.727162 0.747116 0.766214 0.70558 0.734094 0.710806 0.686682 0.669916 0.691097 0.654838 + 0.673654 0.652258 0.638797 0.637248 0.624998 0.61368 0.626461 0.641165 0.658199 0.644357 0.628679 0.631802 + 0.617206 0.615074 0.60457 0.603174 0.605354 0.620244 0.589473 0.150246 0.569346 0.0336846 0.249272 0.593173 + 0.546306 0.0257014 0.0950271 0.592682 0.579912 0.59253 0.602472 0.592897 0.583632 0.584428 0.56796 0.0637246 + 0.517596 0.0243205 0.0187168 0.0169997 0.0144421 0.0143466 0.0167667 0.0142766 0.0186305 0.0110555 0.0105172 0.00618737 + 0.00664736 0.0111928 0.0112163 0.00663896 0.00662496 0.0112619 0.00663781 0.0145015 0.0150007 0.0195847 0.0146676 0.0194054 + 0.0192325 0.0113034 0.0112212 0.00663823 0.00664227 0.0111899 0.0145962 0.014495 0.0122632 0.00703937 0.00695647 0.0108933 + 0.00631998 0.0147029 0.0132321 0.0168716 0.0207707 0.0167808 0.038453 0.51587 0.0338626 0.0203376 0.456494 0.0611455 + 0.624418 0.612995 0.636254 0.649008 0.655619 0.642522 0.664902 0.670518 0.679018 0.650905 0.6767 0.649084 + 0.614415 0.689272 0.705385 0.696883 0.706318 0.738841 0.711592 0.7377 0.791407 0.721857 0.760484 0.765814 + 0.728306 0.0732482 0.0290345 0.0179735 0.0170477 0.0138778 0.0142756 0.0220455 0.00997299 0.0103918 0.0058396 0.00643007 + 0.0113956 0.010062 0.0125729 0.00646816 0.796579 0.741923 0.805922 0.744776 0.700545 0.714112 0.755695 0.813594 + 0.822241 0.759156 0.72845 0.723674 0.701923 0.696086 0.685551 0.678149 0.662776 0.717201 0.739821 0.755319 + 0.822979 0.82862 0.777526 0.839474 0.771053 0.830214 0.781748 0.835487 0.849521 0.842709 0.848639 0.782249 + 0.792136 0.750298 0.745671 0.744977 0.779735 0.719812 0.720544 0.746718 0.780497 0.841296 0.841612 0.78172 + 0.74908 0.784562 0.722627 0.725537 0.700852 0.703358 0.687321 0.706796 0.729659 0.752707 0.757748 0.788973 + 0.793777 0.849013 0.845683 0.84234 0.851638 0.797952 0.852932 0.7957 0.852345 0.792963 0.851977 0.7985 + 0.853994 0.762301 0.769083 0.75563 0.733491 0.724809 0.741685 0.763514 0.735063 0.711381 0.717055 0.696168 + 0.69124 0.677965 0.683255 0.671849 0.689424 0.702378 0.709596 0.717594 0.696409 0.704729 0.727762 0.685172 + 0.678012 0.668027 0.675304 0.666684 0.683696 0.702941 0.693433 0.714272 0.725107 0.738888 0.751038 0.76874 + 0.75898 0.74256 0.770691 0.797653 0.798135 0.85457 0.859359 0.806231 0.785338 0.850737 + </DataArray> + <DataArray type="Float32" Name="p_gas" NumberOfComponents="1" format="ascii"> + 905373 856540 869563 905308 865180 905084 827718 822588 807976 782351 759748 735880 + 780641 734778 776979 733167 752925 792061 822689 863229 862149 904695 904127 861688 + 819832 821013 861833 903387 902495 901453 866745 820183 844832 799447 788205 773291 + 750785 728394 750593 779221 730418 715402 714737 714463 714508 714972 714761 715031 + 714545 714052 714297 714086 713843 713633 713603 713811 714263 714018 713581 713796 + 714198 714044 713592 714425 714352 715223 714424 715096 716491 714659 714096 713845 + 713607 713376 713154 712934 713150 713360 713143 712919 712714 712484 712703 712270 + 712271 712492 712714 712280 711850 712057 711839 711631 711417 711423 711212 711000 + 710789 710996 711206 711626 710788 710999 711417 711842 712051 712064 712494 712288 + 711858 711637 711433 711214 710796 711012 711231 711654 710811 711030 711451 711878 + 712082 712514 712308 712744 712952 713186 712721 713160 712929 713369 713394 713876 + 713423 712978 713214 713664 713454 713696 713910 713242 713005 712771 712565 712798 + 712335 712538 712105 711903 711677 711252 710833 711055 711478 711278 710858 711078 + 711703 711501 712132 711928 712360 712156 712591 712385 711726 711952 711300 710883 + 711119 711528 711336 710996 711085 711754 711568 711982 712180 712615 712408 713056 + 712845 712823 713268 713033 713482 713942 713722 713504 713289 713740 713961 714202 + 714426 714186 714664 714415 714880 714674 715167 719761 715178 734179 716467 714912 + 715140 742589 722005 714636 714843 714374 714160 714122 714586 714334 714797 726015 + 715066 746214 771549 782448 815676 817666 785256 819971 771915 857404 864771 900341 + 899160 856268 855372 897986 896890 854370 895935 814086 811095 767116 811978 767914 + 775676 853075 852676 895158 894557 851763 812049 812762 842905 894089 893705 852146 + 893379 807760 822607 788421 762861 784338 739496 715514 739842 768694 715692 733454 + 714942 715189 714446 713976 714210 714683 714435 713752 713979 714903 714640 715062 + 715395 714829 715010 714171 713767 714125 714447 713839 713940 714688 714790 714355 + 715430 731556 748612 776670 796669 813357 813732 761912 858143 854869 893118 892953 + 841289 863686 844779 892897 713531 713468 713131 713102 713493 713023 712718 712740 + 712359 712359 712672 712441 712865 712636 713305 713076 713521 712203 711791 712003 + 711991 711635 711618 711292 711376 710964 710691 710650 710352 709806 710071 710232 + 710412 710692 710490 710233 710016 710664 710440 710015 709779 709558 709329 709579 + 709799 709363 710221 710004 710637 710417 710614 710201 709789 709590 709379 709168 + 708996 708752 708925 709118 708599 708858 708465 708728 708351 708656 708258 708567 + 708185 708929 709046 709084 709343 709551 709359 709202 709581 709984 709769 710180 + 710397 710595 710379 710582 710167 709966 709752 709542 709956 709746 709334 710161 + 710371 710578 710369 710579 710161 709745 709952 709539 709330 709130 708917 708723 + 708870 709129 708796 708467 708488 708133 708102 708385 708504 708092 + </DataArray> + <DataArray type="Float32" Name="rho_gas" NumberOfComponents="1" format="ascii"> + 10.7573 10.1771 10.3318 10.7565 10.2797 10.7538 9.83461 9.77365 9.60004 9.29557 9.02701 8.74343 + 9.27526 8.73034 9.23175 8.7112 8.94595 9.41095 9.77486 10.2565 10.2437 10.7492 10.7425 10.2382 + 9.74091 9.75494 10.24 10.7337 10.7231 10.7107 10.2983 9.74508 10.0379 9.49871 9.36513 9.18793 + 8.92052 8.65448 8.91824 9.25839 8.67853 8.50012 8.49221 8.48896 8.48949 8.495 8.4925 8.4957 + 8.48993 8.48408 8.48699 8.48448 8.48159 8.4791 8.47874 8.48121 8.48658 8.48367 8.47848 8.48103 + 8.48581 8.48398 8.47861 8.48851 8.48764 8.49799 8.4885 8.49648 8.51305 8.49129 8.48459 8.48162 + 8.47878 8.47604 8.47341 8.47079 8.47335 8.47585 8.47327 8.47062 8.46818 8.46545 8.46805 8.4629 + 8.46292 8.46554 8.46818 8.46302 8.45791 8.46037 8.45778 8.45531 8.45276 8.45284 8.45033 8.44781 + 8.4453 8.44777 8.45026 8.45524 8.4453 8.4478 8.45277 8.45781 8.4603 8.46045 8.46556 8.46311 + 8.45801 8.45538 8.45296 8.45036 8.44539 8.44796 8.45055 8.45559 8.44557 8.44816 8.45317 8.45824 + 8.46067 8.46581 8.46336 8.46854 8.471 8.47378 8.46826 8.47347 8.47073 8.47596 8.47626 8.48198 + 8.4766 8.47131 8.47411 8.47947 8.47696 8.47985 8.48239 8.47446 8.47163 8.46885 8.46641 8.46917 + 8.46367 8.46609 8.46094 8.45854 8.45585 8.4508 8.44583 8.44846 8.45349 8.45111 8.44612 8.44873 + 8.45616 8.45376 8.46126 8.45884 8.46398 8.46155 8.46672 8.46427 8.45644 8.45913 8.45138 8.44642 + 8.44922 8.45409 8.4518 8.44776 8.44882 8.45677 8.45456 8.45948 8.46183 8.467 8.46454 8.47223 + 8.46973 8.46947 8.47476 8.47196 8.4773 8.48276 8.48015 8.47756 8.47501 8.48037 8.483 8.48586 + 8.48851 8.48567 8.49134 8.48839 8.49392 8.49147 8.49732 8.55191 8.49745 8.72321 8.51277 8.49429 + 8.49701 8.82314 8.57856 8.49102 8.49347 8.4879 8.48536 8.4849 8.49041 8.48743 8.49292 8.62621 + 8.49613 8.86621 9.16723 9.29673 9.69153 9.71517 9.3301 9.74256 9.17158 10.1873 10.2749 10.6975 + 10.6834 10.1738 10.1632 10.6695 10.6565 10.1513 10.6451 9.67263 9.6371 9.11456 9.64758 9.12404 + 9.21626 10.1359 10.1311 10.6359 10.6288 10.1203 9.64843 9.65691 10.0151 10.6232 10.6186 10.1249 + 10.6148 9.59747 9.77388 9.3677 9.064 9.31918 8.78639 8.50145 8.7905 9.1333 8.50356 8.71461 + 8.49464 8.49758 8.48876 8.48318 8.48595 8.49158 8.48862 8.48051 8.48321 8.49419 8.49106 8.49608 + 8.50003 8.49331 8.49546 8.48549 8.48069 8.48494 8.48877 8.48154 8.48274 8.49163 8.49284 8.48767 + 8.50045 8.69205 8.8947 9.22808 9.46569 9.66398 9.66843 9.05273 10.1961 10.1572 10.6117 10.6097 + 9.99585 10.262 10.0373 10.609 8.47789 8.47714 8.47313 8.47278 8.47744 8.47185 8.46823 8.46848 + 8.46396 8.46395 8.46767 8.46493 8.46997 8.46725 8.4752 8.47247 8.47777 8.46211 8.45721 8.45973 + 8.45958 8.45535 8.45516 8.45128 8.45229 8.44738 8.44415 8.44365 8.44011 8.43362 8.43677 8.43868 + 8.44083 8.44415 8.44175 8.4387 8.43612 8.44381 8.44115 8.43611 8.4333 8.43068 8.42795 8.43093 + 8.43354 8.42836 8.43855 8.43598 8.4435 8.44089 8.44323 8.43831 8.43343 8.43106 8.42856 8.42605 + 8.424 8.4211 8.42316 8.42545 8.41928 8.42236 8.41769 8.42082 8.41634 8.41996 8.41523 8.41891 + 8.41437 8.4232 8.42459 8.42504 8.42813 8.4306 8.42831 8.42644 8.43096 8.43575 8.43319 8.43807 + 8.44064 8.443 8.44044 8.44285 8.43791 8.43552 8.43299 8.43048 8.4354 8.43291 8.42801 8.43784 + 8.44034 8.44279 8.44031 8.44281 8.43784 8.4329 8.43536 8.43045 8.42797 8.4256 8.42306 8.42075 + 8.4225 8.42558 8.42163 8.41772 8.41797 8.41375 8.41338 8.41674 8.41816 8.41325 + </DataArray> + <DataArray type="Float32" Name="mob_gas" NumberOfComponents="1" format="ascii"> + 0.0350659 0.193363 0.121526 0.0256814 0.14348 0.0299397 0.318691 0.347206 0.430515 0.671563 1.02443 6.75175 + 0.715445 8.17652 0.708318 8.16247 1.16299 0.491564 0.321194 0.151157 0.154704 0.031358 0.0315681 0.152286 + 0.319231 0.319813 0.146409 0.0303419 0.0240009 0.0397688 0.127856 0.319455 0.21444 0.457262 0.505623 0.683061 + 1.41626 14.6809 1.28269 0.634396 9.67669 7300.98 13329.7 14366.7 14610.3 8169.79 13764.6 9565.13 + 14939 15188.6 15049.2 15543.1 15664.7 16202.4 15819.3 15357.8 14765.3 14946.5 15555.1 15150.8 + 14583.5 14802.4 15390.6 14302.9 11870.3 10591.2 14401.5 11946.2 1199.71 14154.5 14695.5 15021.1 + 15331.4 15669.2 15991.4 16355.5 16050.3 15782.4 16223.4 16477.7 16768.7 17241.6 16960.9 17771.3 + 17564.1 17106.3 16705.1 17497.8 18367.5 17941.6 18446.6 18857.1 19414.9 19322.6 19866.5 20390.7 + 21007 20492.7 20050 19017.4 21209.8 20770 19658.3 18665.2 18085.7 18359.8 17498.7 18100.6 + 19024.4 19315.2 20054.6 20383.4 21589.4 21213.4 20905.6 19783.5 22178.3 21894.8 20658.4 19568.1 + 18784.4 17887.1 18594.4 17718.4 17075.7 16924.9 17266.7 16509.5 16716.7 16006.1 16337.9 16102.9 + 16813 17592.5 17508.4 16740 17467.8 17472.4 16707.4 18307.4 18306.7 18355.1 19236 19233.7 + 19294.2 18451.7 19404.6 20341.7 20470.9 21676.3 23060.5 22879.3 21524.5 22779.1 24328 24310.9 + 21443.3 22764.2 20274.1 21433.7 20267.5 21493.1 20325.4 21637.6 22827.3 22985.8 24384.1 26268.2 + 26338.6 24543 26768.7 28578.4 30356.5 24869.9 27391.5 25323.8 23259.2 21872.3 23630.9 20660 + 22178 20455.9 19407.2 19288.4 18360.8 17524.2 18470.3 19589.7 20927.3 19837 18637 18873 + 17782.7 17626.2 16865.4 16765.6 16921.5 18006.8 15802.2 343.123 14441.4 4.10973 1483.13 16059.5 + 12963.6 1.83279 89.402 16025.3 15147.7 16014.7 16715.5 16040.3 15400.7 15455.1 14350.1 27.4024 + 11243.1 1.55399 0.710154 0.532487 0.326809 0.320381 0.510937 0.315713 0.700401 0.146792 0.126403 0.0257851 + 0.0319671 0.15232 0.15328 0.0318466 0.0316461 0.155156 0.0318306 0.330851 0.366132 0.813271 0.342333 0.791196 + 0.770248 0.156879 0.153485 0.0318369 0.0318952 0.152211 0.337371 0.330415 0.200254 0.0379571 0.0366335 0.140446 + 0.0274788 0.344816 0.251491 0.520544 0.969602 0.512229 6.09863 11143.9 4.1746 0.910349 8030.75 24.2386 + 18317.2 17474 19212 20200.1 20721.7 19694.5 21465.2 21921.2 22619.9 20348.9 22428.2 20205.8 + 17577.7 23476.4 24852.2 24122 24933.3 27821 25391.8 27717.4 32766.5 26295.4 29817.3 30317.9 + 26870 41.4087 2.63771 0.629062 0.536924 0.29007 0.315669 1.15858 0.107815 0.121953 0.0216824 0.0289385 + 0.160758 0.110717 0.215774 0.0294553 33270.3 28101.9 34187.3 28363 24435.5 25612.8 29371.2 34947.2 + 35810.5 29693.9 26883.7 26457.4 24554.1 24054.3 23164.2 22548.2 21294 25884.4 27910.7 29336.5 + 35884.7 36452.3 31430 37552.4 30813.9 36613.4 31834.8 37147.3 38580.4 37882.9 38489.9 31883 + 32838.4 28871.6 28445.8 28382.2 31641.9 26115.3 26180.1 28542 31714.9 37738.7 37771.1 31832.4 + 28759.5 32105.8 26364.9 26624 24462.5 24678.3 23313.2 24975.6 26992.8 29094.7 29563.3 32532 + 32998.4 38528.8 38187.3 37845.5 38798.7 33405.9 38931.9 33186 38871.4 32919.4 38833.7 33459.6 + 39041.4 29989.2 30628 29366.2 27337.7 26559.2 28081.3 30103.1 27479.6 25374.5 25872.2 24061.9 + 23643.5 22533.5 22972.9 22030.6 23490.2 24594 25218.9 25919.8 24082.5 24796.8 26822.9 23133.1 + 22537.4 21719.1 22314.1 21610.1 23009.8 24642.5 23829.3 25627.6 26585.8 27826.6 28940.5 30595.7 + 29678.5 28161.2 30780.3 33376.8 33424 39100.9 39595.6 34219.4 32180.9 38706.2 + </DataArray> + <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii"> + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 + </DataArray> + <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii"> + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 + </DataArray> + </PointData> + <CellData Scalars="process rank"> + <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii"> + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 + </DataArray> + </CellData> + <Points> + <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii"> + 0 0 0 0 5 0 3.66025 3.66025 0 5 0 0 + 7.65424 4.09664 0 10 0 0 4.52996 7.95934 0 9.97182 8.46996 0 + 0 10 0 3.30799 12.6456 0 0 15 0 4.4451 17.4535 0 + 7.56491 12.8149 0 9.8537 17.5683 0 13.0835 13.1884 0 14.1111 17.8459 0 + 17.5159 15.6527 0 17.5017 11.5713 0 15.0021 8.4238 0 12.5422 4.26744 0 + 17.5008 4.33008 0 15 0 0 20 0 0 22.4663 4.31026 0 + 20.0006 8.65927 0 24.7673 8.46252 0 27.2977 4.21326 0 25 0 0 + 30 0 0 35 0 0 31.2901 3.63161 0 30.0255 8.44728 0 + 34.814 5.80374 0 34.6484 10.4842 0 27.0465 11.8187 0 30.5113 13.2923 0 + 34.6615 15.5301 0 30.3354 17.9917 0 26.223 15.7578 0 22.1065 12.8385 0 + 21.6904 17.947 0 17.5219 20.4027 0 21.658 22.9035 0 17.3595 25.3866 0 + 25.9686 25.4082 0 25.9837 20.425 0 30.2915 22.9318 0 34.6153 20.4554 0 + 34.5983 25.4372 0 25.9536 30.3892 0 30.2755 27.9129 0 34.5804 30.4163 0 + 30.2589 32.8916 0 34.5585 35.3945 0 25.9383 35.367 0 21.6455 32.8677 0 + 21.6577 27.8792 0 17.3611 30.3198 0 17.3493 35.2799 0 13.0252 32.6839 0 + 13.0307 27.7417 0 8.56562 30 0 8.66262 35.0464 0 9.24601 25.7382 0 + 12.6416 22.3257 0 7.32207 21.8499 0 4.66999 26.7738 0 3.26565 22.5 0 + 0 20 0 0 25 0 0 30 0 4.37181 32.3888 0 + 0 35 0 4.33261 37.5128 0 0 40 0 4.32961 42.5012 0 + 8.67087 40.0586 0 13.0331 37.7003 0 17.3283 40.2668 0 13.0264 42.7357 0 + 8.66389 45.0362 0 12.9945 47.7275 0 17.3029 45.2932 0 17.2793 50.2731 0 + 8.65995 50.1372 0 4.33838 47.5575 0 0 45 0 0 50 0 + 0 55 0 4.32887 52.5923 0 8.66055 55.1634 0 4.33443 57.5621 0 + 8.66491 60.11 0 0 60 0 4.33013 62.5 0 0 65 0 + 4.33243 67.5116 0 8.66546 65.072 0 12.9773 62.6395 0 12.9684 57.7066 0 + 12.9847 67.5822 0 17.2915 65.1441 0 17.2666 60.2155 0 17.2757 55.2447 0 + 12.9735 52.7462 0 21.599 52.7555 0 21.5979 47.7905 0 25.9052 50.2766 0 + 25.9015 55.2386 0 21.5885 57.7202 0 25.8959 60.1933 0 21.5926 62.6736 0 + 21.5855 67.6174 0 25.8834 65.1419 0 30.1818 62.6733 0 30.1917 57.7203 0 + 30.1769 67.6156 0 34.4968 65.1486 0 34.4761 60.2102 0 34.491 55.2533 0 + 30.2034 52.7629 0 30.2149 47.801 0 34.5055 50.2939 0 34.5214 45.3307 0 + 30.2248 42.8374 0 34.5373 40.3664 0 25.9143 45.3109 0 25.9251 40.3418 0 + 21.6097 42.8173 0 21.6402 37.8345 0 30.2384 37.8688 0 38.8708 32.9245 0 + 38.844 37.9048 0 38.8282 42.8635 0 43.1315 40.3986 0 43.1522 35.4378 0 + 47.4364 37.9334 0 51.74 35.4689 0 47.4553 32.973 0 51.7252 40.4237 0 + 47.4225 42.89 0 43.1304 45.3472 0 47.4368 47.823 0 51.7324 45.3647 0 + 43.1502 50.2875 0 38.8238 47.8179 0 38.8205 52.7689 0 43.1451 55.2312 0 + 38.8217 57.7129 0 38.8177 62.6652 0 38.8534 67.5806 0 43.1688 65.0989 0 + 43.167 60.1566 0 47.452 62.627 0 47.4431 67.548 0 51.7234 65.1046 0 + 47.4631 57.6914 0 51.734 60.1711 0 47.4739 52.7512 0 51.7456 55.2364 0 + 51.7399 50.2992 0 56.027 52.7793 0 56.0385 47.8405 0 60.3062 50.3212 0 + 56.0125 57.7149 0 60.2927 55.2567 0 55.9991 62.6492 0 55.9514 67.5167 0 + 60.1365 64.8989 0 60.2533 60.1439 0 64.4773 62.5175 0 63.858 66.4584 0 + 68.665 65.5733 0 64.5308 57.6818 0 68.6798 59.9648 0 68.7727 55.2042 0 + 64.573 52.7983 0 64.5853 47.8662 0 68.8514 50.3437 0 64.5972 42.9291 0 + 68.8662 45.4113 0 60.3192 45.389 0 60.3143 40.4455 0 56.0221 42.9022 0 + 56.0274 37.958 0 56.0423 33.0048 0 60.3268 35.4969 0 64.6077 37.9849 0 + 68.8853 40.471 0 68.9011 35.5274 0 64.625 33.0355 0 68.9265 30.5724 0 + 64.6446 28.0795 0 60.3432 30.5409 0 60.3585 25.582 0 56.0598 28.0487 0 + 64.6644 23.1176 0 68.9514 25.6152 0 68.9712 20.6537 0 64.749 17.6559 0 + 60.3026 20.7413 0 60.041 16.5132 0 55.9939 18.2629 0 56.0528 23.1067 0 + 51.7588 20.6161 0 51.8701 15.8256 0 47.4957 18.1143 0 51.7564 25.563 0 + 47.4989 23.0617 0 47.4809 28.021 0 51.7569 30.5165 0 43.1785 30.469 0 + 43.2072 25.4935 0 38.8979 27.9439 0 38.9144 22.9728 0 38.9263 18.0051 0 + 43.2107 20.5425 0 43.1033 15.7719 0 47.5008 13.0027 0 52.43 11.74 0 + 49.9896 8.41003 0 44.9436 8.3423 0 42.6391 11.7418 0 39.5314 8.247 0 + 38.8475 13.2251 0 42.5 4.33013 0 38.3691 3.67617 0 40 0 0 + 45 0 0 47.5006 4.32206 0 52.5013 4.29319 0 50 0 0 + 55 0 0 57.5002 4.28716 0 60 0 0 55.0027 8.44751 0 + 59.8845 8.64593 0 56.8219 13.1842 0 64.9714 8.47247 0 62.4941 12.9726 0 + 67.4684 12.126 0 62.4972 4.32888 0 67.4877 4.30025 0 65 0 0 + 70 0 0 72.4774 4.34119 0 69.9869 8.40998 0 74.1075 8.30391 0 + 77.6329 5.20206 0 75 0 0 80 0 0 82.7066 4.22277 0 + 85 0 0 81.6442 8.75924 0 85.4962 7.218 0 77.6258 10.762 0 + 81.841 13.3415 0 85.9173 11.1179 0 86.1017 15.6984 0 81.8637 18.2245 0 + 77.5815 15.7288 0 72.8145 12.8043 0 73.2779 18.1914 0 69.5432 16.4254 0 + 73.2586 23.1491 0 77.561 20.6871 0 73.2399 28.1065 0 73.1969 33.0702 0 + 77.4942 30.6086 0 77.5336 25.6472 0 81.7981 28.1488 0 77.4574 35.5623 0 + 81.7392 33.1168 0 81.8393 23.1808 0 86.0566 25.763 0 86.1469 20.7148 0 + 90.1803 18.5403 0 90.0527 23.6258 0 93.8814 21.2607 0 86.5732 30.9962 0 + 85.2219 35.4887 0 91.336 31.5552 0 90.4811 27.9119 0 89.4855 34.7476 0 + 94.1544 33.689 0 93.4225 25.225 0 97.0031 24.298 0 95.694 29.0285 0 + 98.0785 19.509 0 94.5045 16.4049 0 98.9176 14.673 0 95.0156 11.8371 0 + 99.5185 9.80171 0 96.2185 8.10777 0 90.6864 8.09048 0 90.3794 13.3637 0 + 91.1721 3.56345 0 87.3687 3.91634 0 90 0 0 95 0 0 + 94.2846 5.27024 0 96.9079 2.98226 0 99.8795 4.90677 0 100 0 0 + 92.388 38.2683 0 87.7188 38.8831 0 90.3989 42.7555 0 85.4749 42.9681 0 + 82.2905 38.499 0 80.9295 43.7414 0 83.601 47.2642 0 88.1921 47.1397 0 + 85.7729 51.4103 0 81.0689 51.2883 0 79.9289 47.6924 0 76.79 50.2197 0 + 76.867 45.4147 0 73.0081 47.8895 0 77.3336 40.4985 0 73.1141 42.9525 0 + 73.1705 38.0153 0 72.9608 52.8098 0 72.8901 57.5321 0 77.638 55.2532 0 + 83.147 55.557 0 80.3208 59.5699 0 76.3695 59.6387 0 77.301 63.4393 0 + 73.173 62.3325 0 74.0951 67.1559 0 65.7476 70.0774 0 70.7107 70.7107 0 + 67.1559 74.0951 0 59.5699 80.3208 0 63.4393 77.301 0 59.4242 75.2886 0 + 63.2129 73.276 0 60.5575 69.9092 0 55.7628 72.1368 0 51.6184 75.0432 0 + 55.5708 77.7264 0 51.6722 69.9615 0 47.4317 72.4757 0 47.374 77.5005 0 + 51.4922 80.4223 0 55.557 83.147 0 51.4103 85.7729 0 47.2678 82.6897 0 + 43.071 79.966 0 42.9551 85.1667 0 43.1646 74.9449 0 38.898 77.4165 0 + 43.1774 70.0189 0 38.8981 72.4986 0 34.5616 70.0532 0 34.5727 74.9675 0 + 34.6239 79.8753 0 38.8566 82.3503 0 34.6174 84.7886 0 38.6403 87.4035 0 + 34.3259 89.3881 0 38.2683 92.388 0 42.7555 90.3989 0 47.1397 88.1921 0 + 33.689 94.1544 0 30.0322 90.9701 0 29.0285 95.694 0 25.5912 92.4747 0 + 24.298 97.0031 0 20.8932 93.2834 0 19.509 98.0785 0 16.3643 94.31 0 + 14.673 98.9176 0 17.179 89.9456 0 26.8255 88.6589 0 22.1999 88.1345 0 + 17.2692 84.9403 0 21.5415 82.4943 0 25.9666 84.8737 0 30.3376 86.8415 0 + 30.3241 82.278 0 30.253 77.4497 0 25.9202 79.9424 0 25.8941 75.0212 0 + 30.2227 72.5326 0 25.878 70.0884 0 21.567 72.5726 0 17.281 70.086 0 + 17.2628 75.0467 0 21.5726 77.5159 0 17.2499 80.0182 0 12.9434 82.5077 0 + 12.9554 77.5286 0 8.64059 80.0119 0 8.62681 84.9845 0 8.65387 75.0276 0 + 12.9693 72.5528 0 8.66476 70.0488 0 4.33291 72.5176 0 0 70 0 + 0 75 0 0 80 0 4.32928 77.5092 0 4.31611 82.4894 0 + 0 85 0 4.27921 87.4196 0 0 90 0 3.97375 92.3567 0 + 8.87157 90.5902 0 12.9386 87.4845 0 13.4198 91.5118 0 7.04154 95.4599 0 + 11.7453 95.2288 0 9.80171 99.5185 0 4.90677 99.8795 0 3.44133 96.4463 0 + 0 95 0 0 100 0 + </DataArray> + </Points> + <Cells> + <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii"> + 0 1 2 2 3 0 3 4 5 3 2 4 + 4 2 6 7 4 6 6 2 1 1 8 6 + 8 9 6 10 9 8 11 9 10 12 9 11 + 11 13 12 6 9 12 12 7 6 7 12 14 + 12 13 14 13 15 14 14 15 16 16 17 14 + 14 17 18 14 18 7 19 18 20 18 19 7 + 7 19 4 5 4 19 5 19 21 21 19 20 + 20 22 21 20 23 22 24 20 18 23 20 24 + 24 25 23 26 23 25 23 26 27 22 23 27 + 28 27 26 29 28 30 26 30 28 31 30 26 + 32 30 31 31 33 32 31 26 25 25 34 31 + 31 35 33 31 34 35 33 35 36 36 35 37 + 38 37 35 35 34 38 38 34 39 40 38 39 + 39 34 25 25 24 39 39 24 17 18 17 24 + 39 17 16 16 40 39 40 16 41 40 41 42 + 43 42 41 44 45 42 45 40 42 40 45 38 + 45 46 37 37 38 45 37 47 36 47 37 46 + 46 48 47 45 44 46 44 49 50 50 46 44 + 46 50 48 50 51 48 51 50 52 52 53 51 + 49 54 52 49 52 50 49 55 54 55 49 56 + 44 56 49 42 56 44 43 57 56 42 43 56 + 56 57 55 57 58 55 58 57 59 60 59 57 + 61 62 59 59 60 61 60 63 61 64 63 60 + 64 60 43 57 43 60 43 41 64 41 15 64 + 16 15 41 64 15 13 13 65 64 65 63 64 + 66 63 65 65 67 66 11 65 13 11 67 65 + 11 68 67 11 10 68 69 67 68 66 67 69 + 66 69 70 61 63 66 71 61 66 61 71 62 + 70 71 66 71 70 72 71 72 73 72 74 73 + 75 73 74 73 75 76 76 62 73 73 62 71 + 77 59 62 59 77 58 78 58 77 79 78 77 + 79 77 76 62 76 77 76 80 79 81 79 80 + 78 79 82 81 82 79 83 82 81 84 81 80 + 80 85 84 80 76 75 85 80 75 75 86 85 + 74 86 75 86 87 85 87 88 89 89 85 87 + 89 84 85 90 84 89 89 91 90 92 90 91 + 91 89 88 88 93 91 94 91 93 95 94 93 + 94 95 96 94 96 97 97 92 94 91 94 92 + 98 99 92 92 97 98 100 98 97 98 100 101 + 101 102 98 99 98 102 102 103 99 104 99 103 + 90 92 99 99 104 90 84 90 104 104 81 84 + 104 83 81 103 83 104 105 83 103 105 106 83 + 106 105 107 105 108 107 108 105 109 109 110 108 + 103 109 105 109 103 102 102 111 109 102 101 111 + 101 112 111 113 111 112 111 113 110 110 109 111 + 110 114 115 113 116 114 114 110 113 117 114 116 + 117 118 114 115 114 118 119 115 118 120 115 119 + 115 108 110 108 115 120 107 108 120 107 120 121 + 120 119 122 121 120 122 121 122 123 124 121 123 + 124 123 125 126 121 124 127 126 124 121 126 107 + 106 107 126 106 82 83 128 82 106 82 128 78 + 126 128 106 128 126 127 128 127 129 54 129 127 + 128 129 78 58 78 129 129 55 58 55 129 54 + 54 130 52 127 124 130 54 127 130 130 124 125 + 130 125 53 53 52 130 131 53 132 133 132 125 + 125 132 53 132 133 134 132 135 131 136 137 138 + 138 135 136 139 136 140 140 136 134 136 135 134 + 132 134 135 141 134 133 134 141 140 140 142 143 + 142 140 141 141 144 142 144 141 145 133 145 141 + 145 133 123 125 123 133 145 123 122 146 122 119 + 122 146 145 144 145 146 147 146 148 148 119 118 + 119 148 146 118 149 148 150 149 117 118 117 149 + 149 150 151 152 148 149 149 151 152 152 151 153 + 151 154 153 153 154 155 152 153 156 156 153 157 + 158 156 159 147 156 158 147 152 156 152 147 148 + 158 144 147 146 147 144 158 142 144 142 158 160 + 159 160 158 160 159 161 161 162 160 161 163 162 + 159 164 161 165 161 164 157 164 159 157 159 156 + 155 157 153 155 166 157 155 167 166 168 166 167 + 164 157 166 169 164 166 166 168 169 170 169 168 + 168 171 170 170 171 172 169 170 173 174 173 170 + 175 176 173 165 173 176 173 165 169 164 169 165 + 176 163 165 163 161 165 176 177 163 177 176 178 + 179 177 180 179 181 177 182 181 179 181 182 183 + 163 181 162 177 181 163 160 143 142 143 160 162 + 143 139 140 139 143 183 162 183 143 183 162 181 + 184 183 182 183 184 139 139 137 136 184 137 139 + 137 184 185 185 184 186 182 186 184 186 182 187 + 179 187 182 187 179 188 187 189 190 186 187 190 + 191 192 190 192 193 190 186 190 193 185 186 193 + 194 193 192 194 195 193 196 194 192 197 196 192 + 196 198 199 199 200 196 199 201 200 200 201 202 + 200 194 196 194 200 203 202 203 200 203 202 204 + 202 205 204 204 205 206 204 207 203 204 208 207 + 207 209 210 207 210 195 195 203 207 203 195 194 + 210 185 195 193 195 185 185 210 137 137 210 138 + 138 209 211 210 209 138 209 207 208 209 208 212 + 211 209 212 212 213 211 135 138 211 211 131 135 + 131 211 213 53 131 51 213 51 131 51 213 48 + 213 214 48 214 213 212 47 48 214 214 215 47 + 215 36 47 215 214 216 216 217 215 217 216 206 + 212 216 214 216 212 208 208 204 206 208 206 216 + 218 206 205 217 206 218 218 219 220 218 220 221 + 221 222 218 223 222 221 218 222 217 224 215 217 + 217 222 224 215 224 36 36 224 33 33 224 223 + 224 222 223 225 226 223 223 226 32 223 32 33 + 32 226 29 29 30 32 29 226 227 227 226 225 + 227 225 228 229 228 225 221 225 223 225 221 229 + 220 229 221 229 220 230 228 229 231 230 231 229 + 231 230 232 232 233 234 233 232 230 230 235 233 + 233 235 236 235 230 220 220 219 235 205 219 218 + 237 219 205 205 202 237 202 201 237 235 219 237 + 236 235 237 238 236 239 239 236 237 237 201 239 + 239 201 199 240 239 199 239 240 238 241 238 242 + 238 241 236 241 233 236 234 233 241 243 234 241 + 241 242 243 243 242 244 244 242 245 246 242 238 + 242 246 245 246 247 245 245 247 248 249 244 245 + 249 245 248 248 250 249 251 252 250 248 251 250 + 253 251 248 253 254 251 255 253 248 248 247 255 + 253 256 257 256 253 255 258 256 259 259 256 260 + 260 256 255 260 255 261 261 262 260 261 263 262 + 255 247 261 261 247 246 261 246 240 238 240 246 + 240 263 261 199 263 240 198 263 199 262 263 198 + 264 265 262 198 264 262 264 198 197 196 197 198 + 266 197 191 192 191 197 190 189 191 189 267 191 + 191 267 266 267 268 266 197 266 264 266 269 264 + 268 270 269 269 266 268 268 267 271 271 272 268 + 270 268 272 273 269 270 270 274 273 273 265 269 + 265 264 269 260 262 265 265 259 260 265 273 259 + 275 259 273 258 259 275 275 276 258 276 275 277 + 277 278 276 274 277 275 275 273 274 279 274 270 + 270 272 279 272 280 279 279 281 282 279 283 281 + 284 281 283 279 282 274 282 277 274 282 285 277 + 286 285 287 287 285 282 282 281 287 281 284 287 + 288 278 286 278 285 286 277 285 278 278 288 289 + 289 276 278 290 289 288 291 289 290 292 293 291 + 291 290 292 291 293 294 294 295 291 289 291 295 + 276 289 295 295 258 276 258 257 256 257 258 295 + 257 254 253 294 254 257 257 295 294 294 296 297 + 297 254 294 251 254 297 251 297 252 252 297 298 + 297 296 298 299 298 296 300 301 299 299 296 300 + 300 296 294 293 300 294 301 300 293 293 302 301 + 302 293 292 303 301 302 299 301 303 284 283 304 + 304 305 306 306 305 307 304 283 305 283 280 305 + 279 280 283 305 280 308 307 305 308 307 309 310 + 307 310 311 311 306 307 311 310 312 313 312 310 + 310 314 313 309 314 310 313 314 315 316 317 315 + 315 314 316 318 319 316 316 309 318 316 314 309 + 308 309 307 308 318 309 308 280 272 308 272 271 + 320 318 271 271 318 308 319 318 320 188 319 320 + 320 271 267 189 320 267 188 189 187 188 320 189 + 179 180 188 180 319 188 317 316 319 319 180 317 + 178 317 180 177 178 180 175 321 178 176 175 178 + 173 174 175 322 175 174 175 322 321 321 322 323 + 317 178 321 321 315 317 323 315 321 313 315 323 + 312 313 324 324 313 323 323 325 324 323 326 325 + 322 326 323 327 326 328 325 326 327 329 328 172 + 329 327 328 328 326 322 328 322 174 174 170 172 + 172 328 174 172 330 331 330 332 331 172 331 329 + 333 334 335 334 336 335 332 336 334 337 336 330 + 330 336 332 172 171 330 330 171 337 337 171 168 + 337 168 167 337 338 335 335 336 337 338 339 340 + 337 167 338 341 338 167 167 155 341 154 341 155 + 341 342 339 339 338 341 339 343 344 344 340 339 + 340 335 338 335 340 333 340 344 345 345 333 340 + 346 345 344 344 347 346 347 344 343 343 348 347 + 349 347 348 348 343 350 350 351 348 343 339 342 + 342 350 343 341 154 342 352 342 154 154 151 352 + 150 352 151 353 350 352 342 352 350 350 353 351 + 150 353 352 353 150 354 117 354 150 353 354 355 + 355 351 353 355 356 351 357 348 351 348 357 349 + 358 357 356 351 356 357 357 358 359 358 360 359 + 359 360 361 359 361 362 349 359 362 359 349 357 + 347 349 363 363 346 347 362 363 349 360 364 361 + 364 365 366 365 367 366 366 367 368 369 370 368 + 370 371 372 369 371 370 373 371 369 368 367 369 + 367 374 375 375 369 367 373 369 375 376 373 375 + 375 377 376 375 374 378 375 378 377 378 379 380 + 378 374 379 365 374 367 379 374 365 364 360 365 + 379 365 360 379 360 358 380 379 358 358 356 380 + 381 380 356 380 382 378 382 380 381 381 383 382 + 383 381 384 356 355 381 384 381 355 384 355 354 + 116 384 354 354 117 116 384 116 385 116 113 385 + 384 385 383 386 383 385 385 112 386 112 385 113 + 387 386 112 112 101 387 386 387 388 389 388 390 + 389 386 388 383 386 389 389 382 383 382 389 377 + 377 378 382 390 377 389 390 376 377 391 376 390 + 391 390 392 391 393 394 392 393 391 393 392 395 + 396 395 392 392 388 396 388 392 390 396 388 387 + 387 101 100 100 396 387 396 100 397 97 397 100 + 397 395 396 395 397 398 397 97 96 96 398 397 + 398 96 399 95 399 96 399 400 398 400 401 402 + 402 398 400 398 402 395 402 393 395 402 403 393 + 403 394 393 401 403 402 404 403 401 404 405 403 + 406 405 404 406 407 405 405 407 408 394 405 408 + 394 403 405 409 391 394 376 391 409 409 373 376 + 409 410 373 408 410 409 408 409 394 408 411 412 + 412 410 408 371 410 412 373 410 371 371 412 372 + 372 412 413 413 412 411 414 413 411 411 415 414 + 407 411 408 407 415 411 416 415 407 416 407 406 + 417 415 416 417 414 415 + </DataArray> + <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii"> + 3 6 9 12 15 18 21 24 27 30 33 36 + 39 42 45 48 51 54 57 60 63 66 69 72 + 75 78 81 84 87 90 93 96 99 102 105 108 + 111 114 117 120 123 126 129 132 135 138 141 144 + 147 150 153 156 159 162 165 168 171 174 177 180 + 183 186 189 192 195 198 201 204 207 210 213 216 + 219 222 225 228 231 234 237 240 243 246 249 252 + 255 258 261 264 267 270 273 276 279 282 285 288 + 291 294 297 300 303 306 309 312 315 318 321 324 + 327 330 333 336 339 342 345 348 351 354 357 360 + 363 366 369 372 375 378 381 384 387 390 393 396 + 399 402 405 408 411 414 417 420 423 426 429 432 + 435 438 441 444 447 450 453 456 459 462 465 468 + 471 474 477 480 483 486 489 492 495 498 501 504 + 507 510 513 516 519 522 525 528 531 534 537 540 + 543 546 549 552 555 558 561 564 567 570 573 576 + 579 582 585 588 591 594 597 600 603 606 609 612 + 615 618 621 624 627 630 633 636 639 642 645 648 + 651 654 657 660 663 666 669 672 675 678 681 684 + 687 690 693 696 699 702 705 708 711 714 717 720 + 723 726 729 732 735 738 741 744 747 750 753 756 + 759 762 765 768 771 774 777 780 783 786 789 792 + 795 798 801 804 807 810 813 816 819 822 825 828 + 831 834 837 840 843 846 849 852 855 858 861 864 + 867 870 873 876 879 882 885 888 891 894 897 900 + 903 906 909 912 915 918 921 924 927 930 933 936 + 939 942 945 948 951 954 957 960 963 966 969 972 + 975 978 981 984 987 990 993 996 999 1002 1005 1008 + 1011 1014 1017 1020 1023 1026 1029 1032 1035 1038 1041 1044 + 1047 1050 1053 1056 1059 1062 1065 1068 1071 1074 1077 1080 + 1083 1086 1089 1092 1095 1098 1101 1104 1107 1110 1113 1116 + 1119 1122 1125 1128 1131 1134 1137 1140 1143 1146 1149 1152 + 1155 1158 1161 1164 1167 1170 1173 1176 1179 1182 1185 1188 + 1191 1194 1197 1200 1203 1206 1209 1212 1215 1218 1221 1224 + 1227 1230 1233 1236 1239 1242 1245 1248 1251 1254 1257 1260 + 1263 1266 1269 1272 1275 1278 1281 1284 1287 1290 1293 1296 + 1299 1302 1305 1308 1311 1314 1317 1320 1323 1326 1329 1332 + 1335 1338 1341 1344 1347 1350 1353 1356 1359 1362 1365 1368 + 1371 1374 1377 1380 1383 1386 1389 1392 1395 1398 1401 1404 + 1407 1410 1413 1416 1419 1422 1425 1428 1431 1434 1437 1440 + 1443 1446 1449 1452 1455 1458 1461 1464 1467 1470 1473 1476 + 1479 1482 1485 1488 1491 1494 1497 1500 1503 1506 1509 1512 + 1515 1518 1521 1524 1527 1530 1533 1536 1539 1542 1545 1548 + 1551 1554 1557 1560 1563 1566 1569 1572 1575 1578 1581 1584 + 1587 1590 1593 1596 1599 1602 1605 1608 1611 1614 1617 1620 + 1623 1626 1629 1632 1635 1638 1641 1644 1647 1650 1653 1656 + 1659 1662 1665 1668 1671 1674 1677 1680 1683 1686 1689 1692 + 1695 1698 1701 1704 1707 1710 1713 1716 1719 1722 1725 1728 + 1731 1734 1737 1740 1743 1746 1749 1752 1755 1758 1761 1764 + 1767 1770 1773 1776 1779 1782 1785 1788 1791 1794 1797 1800 + 1803 1806 1809 1812 1815 1818 1821 1824 1827 1830 1833 1836 + 1839 1842 1845 1848 1851 1854 1857 1860 1863 1866 1869 1872 + 1875 1878 1881 1884 1887 1890 1893 1896 1899 1902 1905 1908 + 1911 1914 1917 1920 1923 1926 1929 1932 1935 1938 1941 1944 + 1947 1950 1953 1956 1959 1962 1965 1968 1971 1974 1977 1980 + 1983 1986 1989 1992 1995 1998 2001 2004 2007 2010 2013 2016 + 2019 2022 2025 2028 2031 2034 2037 2040 2043 2046 2049 2052 + 2055 2058 2061 2064 2067 2070 2073 2076 2079 2082 2085 2088 + 2091 2094 2097 2100 2103 2106 2109 2112 2115 2118 2121 2124 + 2127 2130 2133 2136 2139 2142 2145 2148 2151 2154 2157 2160 + 2163 2166 2169 2172 2175 2178 2181 2184 2187 2190 2193 2196 + 2199 2202 2205 2208 2211 2214 2217 2220 2223 2226 2229 2232 + 2235 2238 2241 2244 2247 2250 2253 2256 2259 2262 2265 2268 + 2271 2274 2277 2280 2283 2286 + </DataArray> + <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii"> + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 + </DataArray> + </Cells> + </Piece> + </UnstructuredGrid> +</VTKFile>