diff --git a/CHANGELOG.md b/CHANGELOG.md
index c4d29c5139a5304e5bab72d5bb926c231063278e..d211f2d63c941ea0a55e06cc56776ad789892307 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -15,6 +15,8 @@ build can be downloaded from the job artifacts.
- __Projection__: In addition to the L2-projector projecting between different grids added a helper that
computes the L2-projection of analytic functions in to discrete FEM spaces (requires `dune-functions`).
+- __Poromechanics__: Fixed a bug in `PoroElasticLocalResidual`, where the average density between fluid and solid was computed incorrectly, potentially leading to unphysical body forces.
+
### Immediate interface changes not allowing/requiring a deprecation period:
- __Assembler/Newton/PDE/Solver__: We now distinguish between `SolutionVector` and `ResidualType`/`ResidualVector`. The former
diff --git a/dumux/geomechanics/poroelastic/localresidual.hh b/dumux/geomechanics/poroelastic/localresidual.hh
index 4a7f6544f8cdc53b51872ab425d8aa187c4bd9d6..6e0b8e1c27aa44a2ef0eceabd93f8fec81e1ae8b 100644
--- a/dumux/geomechanics/poroelastic/localresidual.hh
+++ b/dumux/geomechanics/poroelastic/localresidual.hh
@@ -92,7 +92,7 @@ public:
const auto& vv = elemVolVars[scv];
const auto phi = vv.porosity();
const auto rhoFluid = problem.spatialParams().effectiveFluidDensity(element, scv);
- const auto rhoAverage = phi*rhoFluid + (1.0 - phi*vv.solidDensity());
+ const auto rhoAverage = phi*rhoFluid + (1.0 - phi)*vv.solidDensity();
// add body force
const auto& g = problem.spatialParams().gravity(scv.center());
diff --git a/test/multidomain/poromechanics/el2p/CMakeLists.txt b/test/multidomain/poromechanics/el2p/CMakeLists.txt
index 635de4e58bbd683926dbae35379a5b6ffc16fd82..252a6f6ef3eabf8941398d72681c7cd9fc38be29 100644
--- a/test/multidomain/poromechanics/el2p/CMakeLists.txt
+++ b/test/multidomain/poromechanics/el2p/CMakeLists.txt
@@ -13,3 +13,17 @@ dumux_add_test(NAME test_md_poromechanics_el2p
${CMAKE_CURRENT_BINARY_DIR}/test_md_poromechanics_el2p_poroelastic-00010.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/test_md_poromechanics_el2p params.input
-Vtk.OutputName test_md_poromechanics_el2p")
+
+dumux_add_test(NAME test_md_poromechanics_el2p_gravity
+ LABELS multidomain poromechanics 2p poroelastic
+ TARGET test_md_poromechanics_el2p
+ CMAKE_GUARD HAVE_UMFPACK
+ TIMEOUT 3600
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/test_md_poromechanics_el2p_gravity_2p-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_poromechanics_el2p_gravity_twop-00010.vtu
+ ${CMAKE_SOURCE_DIR}/test/references/test_md_poromechanics_el2p_gravity_poroelastic-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_poromechanics_el2p_gravity_poroelastic-00010.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_md_poromechanics_el2p params.input
+ -Vtk.OutputName test_md_poromechanics_el2p_gravity -Problem.EnableGravity true -Newton.MaxRelativeShift 1e-8")
diff --git a/test/references/test_md_poromechanics_el2p_gravity_2p-reference.vtu b/test/references/test_md_poromechanics_el2p_gravity_2p-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d92817936501ca1175b9406faf27093869b4a6b7
--- /dev/null
+++ b/test/references/test_md_poromechanics_el2p_gravity_2p-reference.vtu
@@ -0,0 +1,196 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="64" NumberOfPoints="125">
+ <CellData Scalars="S_liq">
+ <DataArray type="Float32" Name="S_liq" NumberOfComponents="1" format="ascii">
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 0.999999 0.999999 1
+ 1 0.999999 0.999999 1 1 1 1 1 1 1 1 1
+ 1 0.999999 0.999999 1 1 0.999999 0.999999 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1
+ </DataArray>
+ <DataArray type="Float32" Name="p_liq" NumberOfComponents="1" format="ascii">
+ 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07
+ 1.50003e+07 1.50003e+07 1.50003e+07 1.50003e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.50016e+07 1.50016e+07 1.5e+07
+ 1.5e+07 1.50016e+07 1.50016e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07
+ 1.5e+07 1.50016e+07 1.50016e+07 1.5e+07 1.5e+07 1.50016e+07 1.50016e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07
+ 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07 1.5e+07
+ 1.5e+07 1.5e+07 1.5e+07 1.5e+07
+ </DataArray>
+ <DataArray type="Float32" Name="rho_liq" NumberOfComponents="1" format="ascii">
+ 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64
+ 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64
+ 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64
+ 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64
+ 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64 1077.64
+ 1077.64 1077.64 1077.64 1077.64
+ </DataArray>
+ <DataArray type="Float32" Name="mob_liq" NumberOfComponents="1" format="ascii">
+ 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437
+ 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437
+ 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437
+ 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437
+ 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437 692.437
+ 692.437 692.437 692.437 692.437
+ </DataArray>
+ <DataArray type="Float32" Name="S_gas" 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 6.7843e-07 6.7843e-07 0
+ 0 6.7843e-07 6.7843e-07 0 0 0 0 0 0 0 0 0
+ 0 6.79485e-07 6.79485e-07 0 0 6.79485e-07 6.79485e-07 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="p_gas" NumberOfComponents="1" format="ascii">
+ 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07
+ 1.50179e+07 1.50179e+07 1.50179e+07 1.50179e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50192e+07 1.50192e+07 1.50176e+07
+ 1.50176e+07 1.50192e+07 1.50192e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07
+ 1.50176e+07 1.50192e+07 1.50192e+07 1.50176e+07 1.50176e+07 1.50192e+07 1.50192e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07
+ 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07
+ 1.50176e+07 1.50176e+07 1.50176e+07 1.50176e+07
+ </DataArray>
+ <DataArray type="Float32" Name="rho_gas" NumberOfComponents="1" format="ascii">
+ 920.521 920.521 920.521 920.521 920.521 920.521 920.521 920.521 920.521 920.521 920.521 920.521
+ 920.521 920.521 920.521 920.521 920.518 920.518 920.518 920.518 920.518 920.529 920.529 920.518
+ 920.518 920.529 920.529 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518
+ 920.518 920.529 920.529 920.518 920.518 920.529 920.529 920.518 920.518 920.518 920.518 920.518
+ 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518 920.518
+ 920.518 920.518 920.518 920.518
+ </DataArray>
+ <DataArray type="Float32" Name="mob_gas" 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
+ </DataArray>
+ <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+ 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+ 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+ 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+ 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+ 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+ 17569.9 17569.9 17569.9 17569.9
+ </DataArray>
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.20019 0.200355 0.200355 0.20019 0.200355 0.200655 0.200655 0.200355 0.200355 0.200655 0.200655 0.200355
+ 0.20019 0.200355 0.200355 0.20019 0.200073 0.200114 0.200114 0.200073 0.200114 0.200156 0.200156 0.200114
+ 0.200114 0.200156 0.200156 0.200114 0.200073 0.200114 0.200114 0.200073 0.199927 0.199886 0.199886 0.199927
+ 0.199886 0.199845 0.199845 0.199886 0.199886 0.199845 0.199845 0.199886 0.199927 0.199886 0.199886 0.199927
+ 0.19981 0.199645 0.199645 0.19981 0.199645 0.199343 0.199343 0.199645 0.199645 0.199343 0.199343 0.199645
+ 0.19981 0.199645 0.199645 0.19981
+ </DataArray>
+ <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
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 250 0 0 0 250 0 250 250 0
+ 0 0 250 250 0 250 0 250 250 250 250 250
+ 500 0 0 500 250 0 500 0 250 500 250 250
+ 750 0 0 750 250 0 750 0 250 750 250 250
+ 1000 0 0 1000 250 0 1000 0 250 1000 250 250
+ 0 500 0 250 500 0 0 500 250 250 500 250
+ 500 500 0 500 500 250 750 500 0 750 500 250
+ 1000 500 0 1000 500 250 0 750 0 250 750 0
+ 0 750 250 250 750 250 500 750 0 500 750 250
+ 750 750 0 750 750 250 1000 750 0 1000 750 250
+ 0 1000 0 250 1000 0 0 1000 250 250 1000 250
+ 500 1000 0 500 1000 250 750 1000 0 750 1000 250
+ 1000 1000 0 1000 1000 250 0 0 500 250 0 500
+ 0 250 500 250 250 500 500 0 500 500 250 500
+ 750 0 500 750 250 500 1000 0 500 1000 250 500
+ 0 500 500 250 500 500 500 500 500 750 500 500
+ 1000 500 500 0 750 500 250 750 500 500 750 500
+ 750 750 500 1000 750 500 0 1000 500 250 1000 500
+ 500 1000 500 750 1000 500 1000 1000 500 0 0 750
+ 250 0 750 0 250 750 250 250 750 500 0 750
+ 500 250 750 750 0 750 750 250 750 1000 0 750
+ 1000 250 750 0 500 750 250 500 750 500 500 750
+ 750 500 750 1000 500 750 0 750 750 250 750 750
+ 500 750 750 750 750 750 1000 750 750 0 1000 750
+ 250 1000 750 500 1000 750 750 1000 750 1000 1000 750
+ 0 0 1000 250 0 1000 0 250 1000 250 250 1000
+ 500 0 1000 500 250 1000 750 0 1000 750 250 1000
+ 1000 0 1000 1000 250 1000 0 500 1000 250 500 1000
+ 500 500 1000 750 500 1000 1000 500 1000 0 750 1000
+ 250 750 1000 500 750 1000 750 750 1000 1000 750 1000
+ 0 1000 1000 250 1000 1000 500 1000 1000 750 1000 1000
+ 1000 1000 1000
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 4 5 7 6 1 8 9 3
+ 5 10 11 7 8 12 13 9 10 14 15 11
+ 12 16 17 13 14 18 19 15 2 3 21 20
+ 6 7 23 22 3 9 24 21 7 11 25 23
+ 9 13 26 24 11 15 27 25 13 17 28 26
+ 15 19 29 27 20 21 31 30 22 23 33 32
+ 21 24 34 31 23 25 35 33 24 26 36 34
+ 25 27 37 35 26 28 38 36 27 29 39 37
+ 30 31 41 40 32 33 43 42 31 34 44 41
+ 33 35 45 43 34 36 46 44 35 37 47 45
+ 36 38 48 46 37 39 49 47 4 5 7 6
+ 50 51 53 52 5 10 11 7 51 54 55 53
+ 10 14 15 11 54 56 57 55 14 18 19 15
+ 56 58 59 57 6 7 23 22 52 53 61 60
+ 7 11 25 23 53 55 62 61 11 15 27 25
+ 55 57 63 62 15 19 29 27 57 59 64 63
+ 22 23 33 32 60 61 66 65 23 25 35 33
+ 61 62 67 66 25 27 37 35 62 63 68 67
+ 27 29 39 37 63 64 69 68 32 33 43 42
+ 65 66 71 70 33 35 45 43 66 67 72 71
+ 35 37 47 45 67 68 73 72 37 39 49 47
+ 68 69 74 73 50 51 53 52 75 76 78 77
+ 51 54 55 53 76 79 80 78 54 56 57 55
+ 79 81 82 80 56 58 59 57 81 83 84 82
+ 52 53 61 60 77 78 86 85 53 55 62 61
+ 78 80 87 86 55 57 63 62 80 82 88 87
+ 57 59 64 63 82 84 89 88 60 61 66 65
+ 85 86 91 90 61 62 67 66 86 87 92 91
+ 62 63 68 67 87 88 93 92 63 64 69 68
+ 88 89 94 93 65 66 71 70 90 91 96 95
+ 66 67 72 71 91 92 97 96 67 68 73 72
+ 92 93 98 97 68 69 74 73 93 94 99 98
+ 75 76 78 77 100 101 103 102 76 79 80 78
+ 101 104 105 103 79 81 82 80 104 106 107 105
+ 81 83 84 82 106 108 109 107 77 78 86 85
+ 102 103 111 110 78 80 87 86 103 105 112 111
+ 80 82 88 87 105 107 113 112 82 84 89 88
+ 107 109 114 113 85 86 91 90 110 111 116 115
+ 86 87 92 91 111 112 117 116 87 88 93 92
+ 112 113 118 117 88 89 94 93 113 114 119 118
+ 90 91 96 95 115 116 121 120 91 92 97 96
+ 116 117 122 121 92 93 98 97 117 118 123 122
+ 93 94 99 98 118 119 124 123
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 8 16 24 32 40 48 56 64 72 80 88 96
+ 104 112 120 128 136 144 152 160 168 176 184 192
+ 200 208 216 224 232 240 248 256 264 272 280 288
+ 296 304 312 320 328 336 344 352 360 368 376 384
+ 392 400 408 416 424 432 440 448 456 464 472 480
+ 488 496 504 512
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/test/references/test_md_poromechanics_el2p_gravity_poroelastic-reference.vtu b/test/references/test_md_poromechanics_el2p_gravity_poroelastic-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..fb72dcd8d38b3835d99b47602179b9425ea59165
--- /dev/null
+++ b/test/references/test_md_poromechanics_el2p_gravity_poroelastic-reference.vtu
@@ -0,0 +1,165 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="64" NumberOfPoints="125">
+ <PointData Scalars="porosity" Vectors="u">
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.200047 0.200175 0.200175 0.20064 0.200035 0.200083 0.200083 0.200091 0.200187 0.200688 0.200093 0.200119
+ 0.200134 0.200492 0.200079 0.200146 0.200047 0.200175 0.200035 0.200083 0.200187 0.200688 0.200093 0.200119
+ 0.200745 0.200152 0.200526 0.200172 0.200187 0.200093 0.200134 0.200492 0.200079 0.200146 0.200526 0.200172
+ 0.200378 0.200165 0.200134 0.200079 0.200047 0.200175 0.200035 0.200083 0.200187 0.200093 0.200134 0.200079
+ 0.200047 0.200035 0.199982 0.199943 0.199943 0.199842 0.199935 0.199811 0.199953 0.199863 0.199982 0.199943
+ 0.199935 0.199811 0.199776 0.19984 0.199935 0.199953 0.199863 0.19984 0.199884 0.199953 0.199982 0.199943
+ 0.199935 0.199953 0.199982 0.199936 0.199799 0.199799 0.199426 0.199784 0.199381 0.199833 0.199498 0.199936
+ 0.199799 0.199784 0.199381 0.199326 0.199461 0.199784 0.199833 0.199498 0.199461 0.199573 0.199833 0.199936
+ 0.199799 0.199784 0.199833 0.199936 0.199953 0.199825 0.199825 0.199358 0.199812 0.199311 0.199866 0.199507
+ 0.199953 0.199825 0.199812 0.199311 0.199254 0.199473 0.199812 0.199866 0.199507 0.199473 0.199622 0.199866
+ 0.199953 0.199825 0.199812 0.199866 0.199953
+ </DataArray>
+ <DataArray type="Float32" Name="u" NumberOfComponents="3" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0.0203967 0.0203967 0.196367
+ 0 0 0 0 0 0 0 0 0 -5.14796e-12 0.0232633 0.223806
+ 0 0 0 0 0 0 0 0 0 -0.0203967 0.0203967 0.196367
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0.0232633 -1.12797e-13 0.223806
+ 0 0 0 -6.92946e-12 5.21607e-13 0.263098 0 0 0 -0.0232633 5.73196e-13 0.223806
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0.0203967 -0.0203967 0.196367 0 0 0 -4.19502e-12 -0.0232633 0.223806
+ 0 0 0 -0.0203967 -0.0203967 0.196367 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 -9.5166e-06 -9.5166e-06 0.247008 0 0 0 -8.22174e-12 -1.94599e-05 0.292708
+ 0 0 0 9.51661e-06 -9.5166e-06 0.247008 0 0 0 0 0 0
+ 0 0 0 -1.94599e-05 1.06135e-12 0.292708 -8.84654e-12 4.37292e-12 0.356343 1.94599e-05 3.83762e-12 0.292708
+ 0 0 0 0 0 0 -9.5166e-06 9.5166e-06 0.247008 -5.53219e-12 1.94599e-05 0.292708
+ 9.5166e-06 9.5166e-06 0.247008 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 -0.0204115 -0.0204115 0.196421 0 0 0
+ -5.43843e-12 -0.0232888 0.223877 0 0 0 0.0204115 -0.0204115 0.196421 0 0 0
+ 0 0 0 0 0 0 -0.0232888 1.35531e-12 0.223877 -5.67541e-12 4.58532e-12 0.263199
+ 0.0232888 3.16176e-12 0.223877 0 0 0 0 0 0 -0.0204115 0.0204115 0.196421
+ -4.04598e-12 0.0232888 0.223877 0.0204115 0.0204115 0.196421 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 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>
+ </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
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 250 0 0 0 250 0 250 250 0
+ 0 0 250 250 0 250 0 250 250 250 250 250
+ 500 0 0 500 250 0 500 0 250 500 250 250
+ 750 0 0 750 250 0 750 0 250 750 250 250
+ 1000 0 0 1000 250 0 1000 0 250 1000 250 250
+ 0 500 0 250 500 0 0 500 250 250 500 250
+ 500 500 0 500 500 250 750 500 0 750 500 250
+ 1000 500 0 1000 500 250 0 750 0 250 750 0
+ 0 750 250 250 750 250 500 750 0 500 750 250
+ 750 750 0 750 750 250 1000 750 0 1000 750 250
+ 0 1000 0 250 1000 0 0 1000 250 250 1000 250
+ 500 1000 0 500 1000 250 750 1000 0 750 1000 250
+ 1000 1000 0 1000 1000 250 0 0 500 250 0 500
+ 0 250 500 250 250 500 500 0 500 500 250 500
+ 750 0 500 750 250 500 1000 0 500 1000 250 500
+ 0 500 500 250 500 500 500 500 500 750 500 500
+ 1000 500 500 0 750 500 250 750 500 500 750 500
+ 750 750 500 1000 750 500 0 1000 500 250 1000 500
+ 500 1000 500 750 1000 500 1000 1000 500 0 0 750
+ 250 0 750 0 250 750 250 250 750 500 0 750
+ 500 250 750 750 0 750 750 250 750 1000 0 750
+ 1000 250 750 0 500 750 250 500 750 500 500 750
+ 750 500 750 1000 500 750 0 750 750 250 750 750
+ 500 750 750 750 750 750 1000 750 750 0 1000 750
+ 250 1000 750 500 1000 750 750 1000 750 1000 1000 750
+ 0 0 1000 250 0 1000 0 250 1000 250 250 1000
+ 500 0 1000 500 250 1000 750 0 1000 750 250 1000
+ 1000 0 1000 1000 250 1000 0 500 1000 250 500 1000
+ 500 500 1000 750 500 1000 1000 500 1000 0 750 1000
+ 250 750 1000 500 750 1000 750 750 1000 1000 750 1000
+ 0 1000 1000 250 1000 1000 500 1000 1000 750 1000 1000
+ 1000 1000 1000
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 4 5 7 6 1 8 9 3
+ 5 10 11 7 8 12 13 9 10 14 15 11
+ 12 16 17 13 14 18 19 15 2 3 21 20
+ 6 7 23 22 3 9 24 21 7 11 25 23
+ 9 13 26 24 11 15 27 25 13 17 28 26
+ 15 19 29 27 20 21 31 30 22 23 33 32
+ 21 24 34 31 23 25 35 33 24 26 36 34
+ 25 27 37 35 26 28 38 36 27 29 39 37
+ 30 31 41 40 32 33 43 42 31 34 44 41
+ 33 35 45 43 34 36 46 44 35 37 47 45
+ 36 38 48 46 37 39 49 47 4 5 7 6
+ 50 51 53 52 5 10 11 7 51 54 55 53
+ 10 14 15 11 54 56 57 55 14 18 19 15
+ 56 58 59 57 6 7 23 22 52 53 61 60
+ 7 11 25 23 53 55 62 61 11 15 27 25
+ 55 57 63 62 15 19 29 27 57 59 64 63
+ 22 23 33 32 60 61 66 65 23 25 35 33
+ 61 62 67 66 25 27 37 35 62 63 68 67
+ 27 29 39 37 63 64 69 68 32 33 43 42
+ 65 66 71 70 33 35 45 43 66 67 72 71
+ 35 37 47 45 67 68 73 72 37 39 49 47
+ 68 69 74 73 50 51 53 52 75 76 78 77
+ 51 54 55 53 76 79 80 78 54 56 57 55
+ 79 81 82 80 56 58 59 57 81 83 84 82
+ 52 53 61 60 77 78 86 85 53 55 62 61
+ 78 80 87 86 55 57 63 62 80 82 88 87
+ 57 59 64 63 82 84 89 88 60 61 66 65
+ 85 86 91 90 61 62 67 66 86 87 92 91
+ 62 63 68 67 87 88 93 92 63 64 69 68
+ 88 89 94 93 65 66 71 70 90 91 96 95
+ 66 67 72 71 91 92 97 96 67 68 73 72
+ 92 93 98 97 68 69 74 73 93 94 99 98
+ 75 76 78 77 100 101 103 102 76 79 80 78
+ 101 104 105 103 79 81 82 80 104 106 107 105
+ 81 83 84 82 106 108 109 107 77 78 86 85
+ 102 103 111 110 78 80 87 86 103 105 112 111
+ 80 82 88 87 105 107 113 112 82 84 89 88
+ 107 109 114 113 85 86 91 90 110 111 116 115
+ 86 87 92 91 111 112 117 116 87 88 93 92
+ 112 113 118 117 88 89 94 93 113 114 119 118
+ 90 91 96 95 115 116 121 120 91 92 97 96
+ 116 117 122 121 92 93 98 97 117 118 123 122
+ 93 94 99 98 118 119 124 123
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 8 16 24 32 40 48 56 64 72 80 88 96
+ 104 112 120 128 136 144 152 160 168 176 184 192
+ 200 208 216 224 232 240 248 256 264 272 280 288
+ 296 304 312 320 328 336 344 352 360 368 376 384
+ 392 400 408 416 424 432 440 448 456 464 472 480
+ 488 496 504 512
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12 12 12 12 12 12 12 12 12
+ 12 12 12 12
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>