From 7aa71af5b545fcee5bdbf62f9f03dfbda6f14b87 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?J=C3=B6rg=20Buchwald?= <joerg.buchwald@ufz.de>
Date: Mon, 17 Oct 2022 21:58:42 +0200
Subject: [PATCH] [T/TRM] removing body force from displacement test
---
.../ThermoRichardsMechanics/Tests.cmake | 17 +
..._force_for_deformation_ts_0_t_0.000000.vtu | 48 +++
.../BodyForce/square.prj | 349 ++++++++++++++++++
.../BodyForce/square_1x1.gml | 29 ++
.../BodyForce/square_1x1_quad8_1e2.vtu | 33 ++
5 files changed, 476 insertions(+)
create mode 100644 Tests/Data/ThermoRichardsMechanics/BodyForce/expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu
create mode 100644 Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
create mode 100644 Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1.gml
create mode 100644 Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1_quad8_1e2.vtu
diff --git a/ProcessLib/ThermoRichardsMechanics/Tests.cmake b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
index 07c6a62321b..4b972582d32 100644
--- a/ProcessLib/ThermoRichardsMechanics/Tests.cmake
+++ b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
@@ -278,3 +278,20 @@ AddTest(
expected_Column_ts_68_t_7200000.000000.vtu Column_ts_68_t_7200000.000000.vtu epsilon epsilon 1e-5 1e-5
expected_Column_ts_68_t_7200000.000000.vtu Column_ts_68_t_7200000.000000.vtu sigma sigma 1e-5 1e-5
)
+# ThermoRichardsMechanics; test for removing body force from displacement equation
+AddTest(
+ NAME ThermoRichardsMechanics_dont_apply_body_force_for_deformation
+ PATH ThermoRichardsMechanics/BodyForce
+ RUNTIME 1
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS square.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu displacement displacement 1e-10 1e-10
+ expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu pressure pressure 5e-2 1e-8
+ expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu temperature temperature 1e-10 1e-10
+ expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu epsilon epsilon 1e-10 1e-10
+ expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu sigma sigma 5e-2 1e-8
+)
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu b/Tests/Data/ThermoRichardsMechanics/BodyForce/expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu
new file mode 100644
index 00000000000..74610fac729
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="27" format="appended" RangeMin="45" RangeMax="121" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="900" format="appended" RangeMin="0" RangeMax="0" offset="324" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="900" format="appended" RangeMin="0.13" RangeMax="0.13" offset="436" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="900" format="appended" RangeMin="1" RangeMax="1" offset="544" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="900" format="appended" RangeMin="4743582.8085" RangeMax="4758015.9156" offset="644" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="900" format="appended" RangeMin="0" RangeMax="0" offset="1388" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="900" format="appended" RangeMin="0.13" RangeMax="0.13" offset="1500" />
+ </FieldData>
+ <Piece NumberOfPoints="341" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="1608" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="1688" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1768" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1852" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1936" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="0" RangeMax="2000000" offset="2024" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="1990000" RangeMax="2000000" offset="2168" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="2416" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="4743416.4903" RangeMax="4758182.4261" offset="2584" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="288.15" offset="3800" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="288.15" RangeMax="288.15" offset="3892" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="999.97851095" RangeMax="999.98269323" offset="3984" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.13" RangeMax="0.13" offset="4144" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="4220" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="4744154.3767" RangeMax="4757443.7201" offset="4292" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4556" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="4632" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6088" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7828" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="8096" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABsAAAAAAAAAIwAAAAAAAAA=eF4z0zPRM9A1tjAz1k1PTUo1sTAxTDTTS8ksKqkEAFqFB3o=AQAAAAAAAAAAgAAAAAAAAIBwAAAAAAAAMwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAdwNwgAABAQAAAAAAAAAAgAAAAAAAACAcAAAAAAAALwAAAAAAAAA=eF7txTENADAIALAfTwuW8ARC5gI787G0T6dObN8c27Zt27Zt27Zt29//ADz+dhQ=AQAAAAAAAAAAgAAAAAAAACAcAAAAAAAAKAAAAAAAAAA=eF7txTENAAAIA7A5w78bJCzYIO3T5OzEtm3btm3btm3b9vsLNU4peQ==AQAAAAAAAAAAgAAAAAAAAIBwAAAAAAAADQIAAAAAAAA=eF7tzL9L1HEcx/GwxRA0yh/npd+vnebioLk4RA3eJlFOkYgOIgQNIQQNkUlq56nX+fMS+gEZKIENGQaBKH2jSaK+g0MkQoggOB04BYrpvQZ7/wsfnq/lxYsXPMK6sC3/fDwY6HzyoznhBeHR/jx4Izj1X9aO/uNuTXm5fvj05G9Yr58fLY4HvbOPL3dMesHxvpI4+UN8fHx8fHx8fHx8x/z9YvlbGfnvE9ZfLZGffSn/3pD1+0rlH8zKr0/i4+Pj4+Pj4+Pju+tnS+X3z8l/k7T+xzL54wvyu4atfz8i//Wi/OoRfHx8fHx8fHx8fHf97Yj82Cf50yPWny+X37As/9ao9e9E5V8N5Jel8PHx8fHx8fHx8d31f0Xlf/0mP5my/osL8sM1+S3PrN9eIX8zlF+QxsfHx8fHx8fHx3fX/14hv2td/qO09dOV8nt+y782Zv2bnvzeP/IPx/Dx8fHx8fHx8fHd9b948vO25feMW7/fl392V37jhPXjVfIrs/L3JvDx8fHx8fHx8fHd9Zeq5L/dk989af0HF+V/+Cu/dsr6TTH5K4fyd6bw8fHx8fHx8fHx3fXfxeQ3n/Zz/u1p69+tlt96xs91NGP9uhr5HUV+zt/I4OPj4+Pj4+Pj47vrv6qRv3VO/vXn1u+8JD8bkV84Y32/Vv6BJ//nDD4+Pj4+Pj4+Pr6b/j+L1isdAQAAAAAAAAAAgAAAAAAAAIBwAAAAAAAAMwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAdwNwgAABAQAAAAAAAAAAgAAAAAAAACAcAAAAAAAALwAAAAAAAAA=eF7txTENADAIALAfTwuW8ARC5gI787G0T6dObN8c27Zt27Zt27Zt29//ADz+dhQ=AQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAGgAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAODPAAqoAAE=AQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAGgAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAODPAAqoAAE=AQAAAAAAAAAAgAAAAAAAAFAVAAAAAAAAHAAAAAAAAAA=eF7twYEAAAAAw6D5U1/hAFUBAAAAAPAYFVAAAQ==AQAAAAAAAAAAgAAAAAAAAFAVAAAAAAAAHAAAAAAAAAA=eF7twYEAAAAAw6D5U1/hAFUBAAAAAPAYFVAAAQ==AQAAAAAAAAAAgAAAAAAAAKAqAAAAAAAAIQAAAAAAAAA=eF7twTEBAAAAwqD1T20JT6AAAAAAAAAAAAAA4GcqoAABAQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAASQAAAAAAAAA=eF5jYGBgaGixc2SgAT2jgTb0hhra0CcqaEM/KKEN/aOANrRALm1ojUza0A6ptKEjEmlDF8TShh4Fo2AUjIJRMAqGGwAAFNGvZw==AQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAmQAAAAAAAAA=eF5jYGBgaGixc2SgAT2jgTb0hhra0CcqaEM/KKEN/aOANrRALm1ojUza0A6ptKEjEmlDF8TShoblu56mwcmH0bB8vKRucPJhNKxc2FM1OPkwGlbOXCkbnHwYDSu33hQNTj6MhpWDLPmDkw+jYeWqTPbg5MNoWDltkj44+TAaVu77JA9OPoyG1SMp8YOTD6MJ1Vek0AAttuljAQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAWwAAAAAAAAA=eF7tjksKwCAMRPVk3v80tdCfUNIiTjZZD4gw2TwegeHtX7+tHCTe4AleJD6gdge1O7jabgvOYk79akmgu4X/LH9B74u08J/l6uS6OrmuTq6rk+veaWT+8hHy9g==AQAAAAAAAAAAgAAAAAAAAKAqAAAAAAAAbwMAAAAAAAA=eF7tmclP1VAUxqkujAsHBOExvT7GjSYOG904JHWjC2GjLpyHxCFRSBwShwQ0CqhMynNeQEJ0g4kDC02MQ8GNbwNhyUbDkg2QKB3erc97ek5ibymkf8A5yc3J7S/f15PzbtvXWyuXy42uMcxZmZtvJ80ZP+818yhc4nPEfzNXuCDuEndCPCvnI4WGacusp5O+z8VWVQ/cIW4zX5S7IQ7zD0WG+UfmrwPY39Ptav+B28Rt5goXxJ0FOPR3MIH9PTqE/T14bz63iFvMFS6Iu8SdEId5X6nh31dyw0m/z/Udqj/wOeI2c4UL4g5xN8TheG856vvGUW90qf7AbeIOc4VnibuL8LYk9nfHJD4ft/QszC3mChfEbeJuiEO/r6dQ/3MW+7vuofr7AbeJW8wVLog7xJ0IfWMV3n+bNd3vs56ezy3iFnOFC+IOcTfEYf2frMH1q+djfwseq3rgNnGbeSR3iLshLqe5A3X0/y2lm9DvZU//87/EbeIuc4WLEHdCXJNjTL7fLZW5Rb7faX5WeabA8I83dCTNJTJfuxNfL0+X+yH1kOulHuq9ynrWU3jEPeJuiEN8l+/HkFPy/RjypVb1/F/WGv66nHmB6/NCW3w91PVZ6qGuaamH+XnWs57CI54lLiL8Pxah/7cB9D8T2F/Q5HhfjOuzZxCvjxN34+uhnnfFeP10D+L1c5z1rKfwiAviIsL/dQL9jw2h/6HA/owmx6sSXJ8bP+H63H8/vh7qelmCHPRQ3z7Ws57CIw55g+TZCP/+UtTnjaB/Q2B/C9bl8zJcn2MZfL7v6YyvzyM91DWawfv7btazngKOPyvD5ztwEeGfLkf//nH039Wl8q4KfD9qmsD75/bu+Ho4f2cF1tc4gfVtYz3rKTzigrgX4d+eRP+dk+i/NbC/qslxS8f1uXoK1+fmB/H1UM9NHf9/rJrC+/sm1rOewiPuEXcj/G+k0P/XLPqvD+xPa3JcqcTn+1sHn+91vfH1UNflSrx+3jh4/dSynvUUHvEscRHh31SF/i2a7udUWuXnqnF9NizX/ftn6aP4eqjrbDXtf0k9zEtYz3oKj7ggno3wP1WD/ql89C8MfB/R5DhSS/ufCd3fh135JL4e6jpcixz0UN8K1rOewiMOeTqB35eCHPTwfQk4fF+CefD7EvPF+T9Nw69RAQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAJAAAAAAAAAA=eF5LSwMCpiKHtFF6lB6lhyzNMApGwSgYBaNgFAxDAACO6kYyAQAAAAAAAAAAgAAAAAAAAKgKAAAAAAAAIwAAAAAAAAA=eF5LSwMCpiKHtFF6lB6lR+lRepQepUfpUXqUHkQ0AFzpl0g=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAVwAAAAAAAAA=eF7b8s+l7459v8MWKtEGIYtm3aYivWbl/2W3qEhr/4/ZfJOK9PKQXftvUJFWWSV+5joV6QX/S25coyItF3rpyVUq0tNX6X+8QkValKH372Uq0gAvVaCWAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAGAAAAAAAAAA=eF5bUmDLdX3xAfslo/QoPUpj0ADiun7gAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAFQAAAAAAAAA=eF5jYACBD/YMo/QoPUpj0AD4RXZdAQAAAAAAAAAAgAAAAAAAAIAMAAAAAAAAowAAAAAAAAA=eF5jYGBgOCXsfBBIMQT0yIHpqlY/MA0Do/L45feJQuQ/zIbI57WPypMiv0kcIj9hNUQ+qXNUnhT5ZZIQeYPdEPmw7lF5UuRnSUPkL5yCyHv1jsqTIt8nC5EvuAWRt+sflSdFvkkeIi/wCiJvNHFUnhT5MkWI/IafEHm1yaPypMhnKUPkAzjlwbTU1FF5UuTjVCHyHyQg8nzTR+WJlQcAX+swFQ==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAGAAAAAAAAAA=eF77s/LjJd+kAPs/o/QoPUpj0AAu19DoAQAAAAAAAAAAgAAAAAAAAPgfAAAAAAAAIwQAAAAAAAA=eF51mEFqHTEQRH2nZO/b5E7JDf4Rsv4rgQMDHwyCLGQUwogJg44Q/EUxdFVXr6zyo9WqHs186eUljx/fP+Pna64X0b8842H4KvqKZvJ00b8946/JP0T/9fYZ/8y80+i6boxzvYgOH3K+ir6imTxddPiQ5x+iw4d83ml07XP0gfUievSB+Sr6imbydNGjD5x/iB594Hmn0fW5hg+5XkTHOOer6CuaydNFhw95/iE6fMjnnUbXfRyfB9aL6Mib81X0Fc3k6aLH54HzD9Hj88DzTqMjGuXH/mG9iB79Zb6KHsecp4se35Ocf4ge35M87zS6vqfjvmC9iA4fcr6KvqKZPF30uC84/xA97guedxpdv0vwIdeL6Cvvw/BV9BXN5OmiY5znH6LDh3zeaXT9Dsf3A+tFdPiQ81X0Fc3k6aKj7jz/ED2+H3jeaXT93QEfcr2IDh9yvoq+opk8XXT4kOcfomOczzuNfv0ff8X3JOtFdPiQ81X0Fc3k6aLH55jzD9HjOnneafLcjC/6/6/POu/Cx33E+TbDF+HXet+Fj/vu4td6f5v8Vfg1/jDrbcI/07/9MfV0s95d+Liv2c/D1D+MP6fw8T3A9Uzjw+01Pu9XcJ8xZh7j/LnYDF+ER9+ZR9+ZR9/z/FV49D1fbxMefc/r6Wa9u/Doe+7nYeofxp9TePQ9r2caH9B3PUes2KjOu/CxDt7fm+GL8LHvF798fggf+875q/Cx77zeJnzsO9fThV+xCx/f2+znYeofwse+X3x8z3M9U3T0Pfp6BfqAMfrOPPrO/IrNnjeZx5h5jJlH3/P8VXj0PV9vEx59z+vpZr27Pecyj77n9Q/jz2nPxXk90/iA/a7nRP5+xv2u5+j8e7sZvgiP+pjHOpmP+13P48zH/a7ndObR97yebta7C4++534epv5h/DmFj79vuZ5pfMD3Xs/F8Anj+J3QewPmV2yGL8LH51LvGZiP69H7B+bjWO8lmEffmUff8/XuwqPvuZ+H8Oh77s8pfPxdx/VM4wPe83oPENfN73m9J2F+xWb4Ijz6zjz6znx8z+t9C/NYf77eJjzGeT3drHcXHn1nPr7n9T4n9+cUHn3P65nGh5u991ixU513c9+C33V8jtnsfRHz6Dvza50P4dH3PH8VHn3P19uEh195PV14jJmPY/bzsPdXzC/utPdazK+YoqPvcZ4r+HwZv++D8uAcx+fRzfBF+LWed+HRd+bj953zV+HR93y9TXj0Pa+nm/XuwqNPuZ+HqX8Yf07h4zmO65nGh5u918I8GKPvzKPvzK/Y7H0g8+g788vnh/Doe56/Co++5+ttwqPveT3drHe395DMox95/cP4c9p7y7yeaXzA77pJ9dxFX7GJvvK+ix5/d/E94Yfo8Jn1Fbvo8Tnmek7R/wO2lrzJAQAAAAAAAAAAgAAAAAAAAAAZAAAAAAAA+AQAAAAAAAA=eF5d11f86HMdx3F/e+89Dw7H3nv+bcdelZUOKSO0yF7tQcMoyjhlFA1SlJFQCJFVshVFRmXvceH1vDi/383z8bv4vK9eN9/JJnv/G8kZc4Y8Mo/Jb+ZR6W7ynCntHJ3H5bfTnrspcua0c2yekKekPXdT5ixp5/j8fJ6WJ6S7qXLWtHNifjG/m/bcTZ2zpZ0v5JfzjLTnbpqcPe18Kb+a30977qbNOdLOV/LreVbaczddzpl2vpYn5Tlpz930OVfa+UaenBPTns70Mk/OnXr7Tp6X30p3epk37ejt1Lwg7bnTy3xpR2+n54/Tnju9zJ929Pa9vCjtudPLAmlHb2fmT9OeO70smHb09oP8edpzp5eF0o7ezs5L0p47vSycdvR2bv4y7bnTyyJpR28/zF+lPXe6G5N2JuaP8vK0pzO9LJaLpt4uzN/m+elOL4unHb39JK9Ke+70Mjbt6O3ivCbtudPLEmlHbz/La9OeO70smXb09ou8Lu2508u4tKO3S/OGtOdOL0ulHb1dln9Me+70snTa0duv86a0504vy6QdvV2Rf0p77sbksmlHb7/JW9OezvSyfC6Xers6b88r051eVkg7evtd/iXtudPLimlHb7/Pu9KeO72slHb0dn3ek/bc6WXltHNd/iH/mvbc6WWVtKO3G/O+tOdOL6umHb3dnPenPXd6WS3t6O2WfDDtudPL6mlHb7flw2nPne7WSDt6+3M+mvZ0ppe1cs3U2535z7wj3ell7bSjt7vzibTnTi/rpB293Zv/Tnvu9LJu2tHb3/KptOdOL+ulHb39PZ9Oe+70sn7a0dsD+Wzac6eXDdKO3h7K/6Y9d3rZMO3o7ZH8f9pzp5fRtKO3x/KFtOdOdxvlaOrtH/lSPpY608smuXHq7V/5aj6e7vSyadrR25P5etpzp5fN0o7e/pNvpj13etk87ejtmXw77bnTyxZpR2/P5btpz51etkw7evtfjvRgs+dOL+PTjt6ezym6t+dOL1ulHb29mFN1b8/daG6ddvT2ck7Tvb3R1N02aUdvr+R03dvTmV62y21Tb2/kjN2/1r87vWyfdvT2Vs7cvT13etkh7ejtnZy1e3vu9LJj2tGbh/rsac+dXnZKO3qbvLs5PfjTnV52Tjt6mzLnTnvu9LJL2tHb1DnvYM+dXj6QdvQ2bc4/2HOnlw+mHb1NnwsO9tzp7kNpR28z5MKDPZ3pZbfcNfU2Sy6aM41MeqeX3dOO3mbLxQd77vSyR9rR2xy5xGDPnV72TDt6myvHDfbc6eXDaUdv8+TSgz13etkr7ehtvlx2sOdOLx9JO3pbIJcf7LnTy4S0o7eFcsXBnju97J0TUm+L5MqDPXe62yft6G1MrjrY05le9s2Ppt7G5hq52Mikd3r5WNrR25K51mDPnV4+nnb0tlSuM9hzp5f90o7elsn1Bnvu9LJ/2tHbcrnBYM+dXg5IO3pbIUcHe+70cmDa0dtKufFgz92E/ETa0dsquelgz93eeVDa0dtquflgz53uDk47els9txzs6Uwvn8xDUm9r59a55sikd3r5VNrR27q5zWDPnV4+nXb0tn5uO9hzp5fPpB29bZjbDfbc6eWzaUdvG+X2gz13ejk07ehtk9xhsOdOL4elHb1tljsO9tzp5XNpR29b5E6DPXd6OTzt6G187jzYc6e7I9KO3rbKXQZ77wEQhC2+AQAAAAAAAAAAgAAAAAAAACADAAAAAAAApgAAAAAAAAA=eF4txRFghAAAAMC2SRiGYRg+huFjGIZhGIZh+BiGYRiGYfgYhoPu5MLgETl24tSZX85d+O3SlWs3bt259+DRH0+evXj15t2HT399+Xbw8xQ6cuzEqTO/nLvw26Ur127cunPvwaM/njx78erNuw+f/vry7eD3KXTk2IlTZ345d+G3S1eu3bh1596DR388efbi1Zt3Hz799eXbwd9T6MixE/8DLbMvQA==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF4TF6c9AADGKQj9
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
new file mode 100644
index 00000000000..09b4c1c342e
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
@@ -0,0 +1,349 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <geometry>square_1x1.gml</geometry>
+ <mesh>square_1x1_quad8_1e2.vtu</mesh>
+ <processes>
+ <process>
+ <name>BodyForceTest</name>
+ <type>THERMO_RICHARDS_MECHANICS</type>
+ <mass_lumping>false</mass_lumping>
+ <integration_order>3</integration_order>
+ <constitutive_relation id="0">
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>YoungsModulusClay</youngs_modulus>
+ <poissons_ratio>PoissonsRatioClay</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 -9.81</specific_body_force>
+ <apply_body_force_for_deformation>false</apply_body_force_for_deformation><!-- default is true //-->
+ <initial_stress>Initial_stress</initial_stress>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>bulk_modulus</name>
+ <type>Constant</type>
+ <value>2.15e9</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>4181.3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Linear</type>
+ <reference_value>999.1</reference_value>
+ <independent_variable>
+ <variable_name>phase_pressure</variable_name>
+ <reference_condition>1.0e5</reference_condition>
+ <slope>4.6511627906976743e-10</slope>
+ </independent_variable>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>0.001</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>4e-4</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>995</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>1.5e-5</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2689.65517241379</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.13</value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Parameter</type>
+ <parameter_name>IntrinsicPermClay</parameter_name>
+ </property>
+ <property>
+ <name>relative_permeability</name>
+ <type>RelativePermeabilityVanGenuchten</type>
+ <residual_liquid_saturation>0
+ </residual_liquid_saturation>
+ <residual_gas_saturation>0
+ </residual_gas_saturation>
+ <exponent>0.6</exponent>
+ <minimum_relative_permeability_liquid>1e-6
+ </minimum_relative_permeability_liquid>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>SaturationVanGenuchten</type>
+ <residual_liquid_saturation>0.0
+ </residual_liquid_saturation>
+ <residual_gas_saturation>0
+ </residual_gas_saturation>
+ <exponent>0.6</exponent>
+ <p_b>20000000</p_b>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>1.0</cutoff_value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="BodyForceTest">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-12</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>864000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1826</repeat>
+ <delta_t>86400</delta_t>
+ </pair>
+ <pair>
+ <repeat> 1</repeat>
+ <delta_t>21600</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>dont_apply_body_force_for_deformation</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <timesteps>
+ <pair>
+ <repeat>1000</repeat>
+ <each_steps>10</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>temperature</variable>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>b0</basis_vector_0>
+ <basis_vector_1>b1</basis_vector_1>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>b0</name>
+ <type>Constant</type>
+ <values> 0.829037572555042 0.559192903470747</values>
+ </parameter>
+ <parameter>
+ <name>b1</name>
+ <type>Constant</type>
+ <values>-0.559192903470747 0.829037572555042</values>
+ </parameter>
+ <parameter>
+ <name>ThermalConductivityClay</name>
+ <type>Constant</type>
+ <values>2.75862068965517 0 0 1.49425287356322</values>
+ <use_local_coordinate_system>true
+ </use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>IntrinsicPermClay</name>
+ <type>Constant</type>
+ <values>1e-20 0 0 1e-20</values>
+ <use_local_coordinate_system>false</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>YoungsModulusClay</name>
+ <type>Constant</type>
+ <value>8e9</value>
+ </parameter>
+ <parameter>
+ <name>PoissonsRatioClay</name>
+ <type>Constant</type>
+ <value>0.35</value>
+ </parameter>
+ <parameter>
+ <mesh>square_1x1_quad8_1e2</mesh>
+ <name>Initial_stress</name>
+ <type>Function</type>
+ <expression>-2.5e6+-1.0e4*y</expression><!-- apply_body_force_for_deformation=false //-->
+ <expression>-0.5e6+-1.0e4*y</expression><!-- allows for constant total stress //-->
+ <expression>-4.0e6+-1.0e4*y</expression><!-- despite nonzero body force //-->
+ <expression>0.0</expression>
+ </parameter>
+ <parameter>
+ <name>displacement_ic</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Function</type>
+ <expression>2.e6+-1.0e4*y</expression>
+ </parameter>
+ <parameter>
+ <name>pressure_high</name>
+ <type>Constant</type>
+ <value>2.0e6</value>
+ </parameter>
+ <parameter>
+ <name>pressure_low</name>
+ <type>Constant</type>
+ <value>1.99e6</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <values>288.15</values>
+ </parameter>
+ <parameter>
+ <name>dirichlet</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pressure_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_low</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_high</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>100</max_iter>
+ <linear_solver>linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1.gml b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1.gml
new file mode 100644
index 00000000000..fe082f12d74
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1.gml
@@ -0,0 +1,29 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
+<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
+ <name>square_1x1_geometry</name>
+ <points>
+ <point id="0" x="0" y="0" z="0" name="origin"/>
+ <point id="1" x="0" y="1" z="0"/>
+ <point id="2" x="1" y="0" z="0"/>
+ <point id="3" x="1" y="1" z="0"/>
+ </points>
+ <polylines>
+ <polyline id="0" name="left">
+ <pnt>0</pnt>
+ <pnt>1</pnt>
+ </polyline>
+ <polyline id="1" name="right">
+ <pnt>2</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ <polyline id="2" name="bottom">
+ <pnt>0</pnt>
+ <pnt>2</pnt>
+ </polyline>
+ <polyline id="3" name="top">
+ <pnt>1</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ </polylines>
+</OpenGeoSysGLI>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1_quad8_1e2.vtu b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1_quad8_1e2.vtu
new file mode 100644
index 00000000000..32a0c26bd0d
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_1x1_quad8_1e2.vtu
@@ -0,0 +1,33 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="341" NumberOfCells="100">
+ <PointData>
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.4142135624">
+ AQAAAACAAAD4HwAAAwQAAA==eJx1WEuK3TAQnDsl+9wmd0pu8I6Q9awEEzAEBgRZvEEhWDgYHSHMM81QP+1ULlqtakvq7qcnP75/ex8/vni8Cf7pMX4Ffhf8GvdgZwj+9TH+BvtT8J8v7+NfWHcFXPddc483wUsHz++CX+Me7AzBSwdvfwpeOvh1V8A1zqgD401w1IH5XXDUge0MwVEHtj8FRx143RVw/a9LB483wWvu+V3w0sHbGYKXDt7+FLx08OuugOs5xv+B8SY4/g/M74KXDt7OEBz/B7Y/Bcf/gdddAf/wi3XweBMc9WV+FxznbGcIjvck25+C4z3J666A6z2N54LxJjieC+Z3wcsPb2cIjueC7U/B8Vzwuivg+i6VDh5vgpcOnt8FLx28nSF4zb39KXjp4NddAdd3GO8HxpvgeD8wvwteOng7Q3C8H9j+FBzvB153BVzzjtLB403w0sHzu+Clg7czBC8dvP0peM39uivgH99ZB483wfGeZH4XvHTwdobg+B+z/Sk47pPXXcHOLeii3z8//HwO56rOEdvbYt7G/Gu/rzG/Yf61398xD2D+NX+L7yXzH+Zf/sR3xe93j/ev1/OI95TX54zn2fuzgg63mFdznGue6g7/X2yxHmF+xT3VKcyvuKf6hfkV91TXML/inuodv9891kFezyPWR16fM9ZN3p8VdKi4ax1xjS3EXess5tc81V/Mx7hrXcZ8jLvWa8zHuGsdx3yMu9Z3fr+78PHeZj2P4P8M+pzCx3ue/VlBh1usm/jerbgzH/NGvqe3WG8yv+bMxzzzVeLu7XfhV9z9fu/Cr7h7f0bY7x7rXK/nEetfr88Z62Lvzwo61HnXOpHfTzzvWkf793YL/CZ8PO9adzMfz7vW48zH8651OvPxvGv97ve7Cx/zVdbzCP7PoM8pfMxv2Z8VdKj3XutizoPwndC+gc+btsBvwsf/UvsMzMf9aP+B+TjXvgTzMa/TfoXf7y58zOtYzyP4P4M+p/Axr2N/VtDhFvsAnM/iPa99Ep//brF/wny857Wvwny857Xfwny857UPw3y857U/4/e7x76N1/OI/Ryvzxn7PN6fFXS4xb5H7QP9fA79lsrruI7ZYr+I+RV35uN7tkvcvf0u/Iq73+9d+KWX92eE/e6xT8X8irv3fwZ9ztjX8v6soEO979rn4foS33ftg/l6dAv8Jnx837Vvxnx837Wfxnx837XPxnx837X/5ve7C7/i5PU8gv8z6HMKH+s49mcFHW6xr4XvyUfcU9+P+dfYYj+Q+RX31CdkfsU99Q+ZX3FPfUXmV9xTv9Hvd499SK/nEfuTXp8z9i29PyvoUHndkvgyfo1N8IoX45h3cZ/wTXB8h3jdPfh5BH9Owf8Dtpa8yQ==
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="340">
+ AQAAAACAAAAAGQAAuwQAAA==eJxdl0XUkFUYBvno7s6f7u7+aUFsVGxssRW7JezuFluxMcDG7ha7G7u7Fj6z+OZuZvXMas57zq1S5f9XhPXDeuHB4WHhKeEh2lUNG8hzaHhEeJp87KqFDeU5PDwqPEM+dtXDRvIcGS4Mz5KPXY2wsTxHh4vDc+RjVzNsIs+i8JjwPPnY1QqbyrMkPC68QD52tcNm8hwbnhBeJB+7OmFzeY4PTwovkY9d3bCFPCeGJ4dL5aMzemkVtgzp7fTwivBU7eiltTz0dmZ4lXzs6KWNPPR2dniNfOzopa089HZueJ187OilnTz0dn54vXzs6KW9PPR2YXijfOzopYM89HZxeLN87Oilozz0dml4q3zs6KWTPPR2WXibfOzorkKepeHl4R3y0Rm9dAk7h/R2dbgyvFI7eukqD71dG94tHzt66SYPvS0L75WPHb10l4febgjvl48dvfSQh95uClfJx45eespDb7eED8nHjl56yUNvy8NH5GNHL73lobfbw8fkY0cvfeShtzvDJ+RjVxH2lYfeVoRPyUdn9NI/7BfS2z3hs+Fd2tHLAHno7b7wefnY0ctAeejtgfBF+djRyyB56O3B8GX52NHLYHlWhQ+Hq+VjRy9D5KG3R8PX5GNHL0PlobfHwzfkY0cvw+ShtyfDt+RjRy/D5aG3p8N35GNHdyPkobdnwvfkozN6GRWODOnthfDD8Dnt6GW0PPT2UvixfOzoZYw89PZK+Kl87OhlrDz09mr4uXzs6GWcPPT2eviFfOzoZbw89PZm+JV87Ohlgjz09nb4jXzs6GWiPPT2bvidfOzopVIeens//EE+dnQ3SR56+yD8ST46o5cp4eSQ3j4Jfw0/0o5epspDb5+Fv8vHjl6myUNva8I/5WNHL9Plobcvw7/lY0cvM+Sht6/Df+VjRy9ryUNv34ZFUfaxo5eZ8tDb92G1ouxjRy+z5KG3H8MaRdnHrjJcWx56+zmsVZR97Ohutjz09ktYpyj76Ixe1g3XCentj7B+9r9pRy/ryUNvf4UNi7KPHb2sLw+9/RM2Lso+dvSygTz0xke9aVH2saOXDeWht6rZNefDX5R39LKRPPRWPWxZlH3s6GWOPPRWM2wtHzt62VgeeqsdtpWPHb1sIg+91Q3by8eO7jaVh97qhR3lozN62SycG9Jbo7Bz2KAo7+hlc3norUnYVT529LKFPPTWLOwuHzt62VIeemsR9pSPHb1sJQ+9tQp7y8eOXraWh97ahH3lY0cv28hDb+3C/vKxo5d58tBbh3CgfOzoZVt56K1TOFg+dnS3nTz0VhEOlY/O6GWHcPuQ3rqFI8IuRXlHLzvKQ289wlHysaOXneSht17hGPnY0cvO8tBbn3CcfOzoZRd56K1fOEE+dvQyXx56GxBWyseOXnaVh94GhZPlYzcv3E0eehsSTpWPHb3sLg+9DQuny8eO7vaQh96Gc+/lozN62SvcM6S30dzbcGRR3tHL3vLQ29hwtnzs6GUfeehtPPdaPnb0sq889DaRey0fO3pZIA+9TeJey8eOXvaTh96mcK/lY0cv+8tDb9O41/Kxo5cD5KG3Gdxr+djRy4Hy0NtM7rV87OjuIHnobRb3Wr7/ABCELb4=
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="800">
+ AQAAAACAAAAgAwAApgAAAA==eJwtxRFghAAAAMC2SRiGYRg+huFjGIZhGIZh+BiGYRiGYfgYhoPu5MLgETl24tSZX85d+O3SlWs3bt259+DRH0+evXj15t2HT399+Xbw8xQ6cuzEqTO/nLvw26Ur127cunPvwaM/njx78erNuw+f/vry7eD3KXTk2IlTZ345d+G3S1eu3bh1596DR388efbi1Zt3Hz799eXbwd9T6MixE/8DLbMvQA==
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="23" RangeMax="23">
+ AQAAAACAAABkAAAADAAAAA==eJwTF6c9AADGKQj9
+ </DataArray>
+ <DataArray type="Int64" Name="faces" format="binary" RangeMin="0" RangeMax="0">
+ AQAAAACAAAAgAwAAEAAAAA==eJxjYBgFo2AU4AIAAyAAAQ==
+ </DataArray>
+ <DataArray type="Int64" Name="faceoffsets" format="binary" RangeMin="1" RangeMax="100">
+ AQAAAACAAAAgAwAArwAAAA==eJwtxdciAgAAAEAjkYRskazskpUtFJGRUbYI//8NHrp7ueamhha3OuA2B93uDofc6bC7HHG3e9zrqPvc7wEPesjDHvGoYx7zuOOecMKTnvK0ZzzrpOc87wUvesnLXnHKaa864zWve8Ob3nLW297xrve87wMf+sg5H/vEp8674DOfu+gLX7rkK1/7xrcu+873fvCjK676yc9+8avf/O4Pf7rmL3+77h//+s//+ewTuw==
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab