From 7606f8d87e083f2149e6a77e492443e31eacaf2f Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Fri, 1 Dec 2017 13:45:53 +0100
Subject: [PATCH] [test][richards] Cleanup and add cc parallel test
---
.../richards/implicit/CMakeLists.txt | 40 +-
.../richards/implicit/test_richardslens.input | 4 -
.../richards/implicit/test_richardslens_fv.cc | 4 +-
.../richardslenscc-reference-parallel.vtu | 407 ++++++++++++++++++
4 files changed, 434 insertions(+), 21 deletions(-)
create mode 100644 test/references/richardslenscc-reference-parallel.vtu
diff --git a/test/porousmediumflow/richards/implicit/CMakeLists.txt b/test/porousmediumflow/richards/implicit/CMakeLists.txt
index 36aed5fe8a..617f3250d6 100644
--- a/test/porousmediumflow/richards/implicit/CMakeLists.txt
+++ b/test/porousmediumflow/richards/implicit/CMakeLists.txt
@@ -1,33 +1,43 @@
add_input_file_links()
+add_executable(test_richards_tpfa EXCLUDE_FROM_ALL test_richardslens_fv.cc)
+target_compile_definitions(test_richards_tpfa PUBLIC TYPETAG=RichardsLensCCProblem)
+
+add_executable(test_richards_box EXCLUDE_FROM_ALL test_richardslens_fv.cc)
+target_compile_definitions(test_richards_box PUBLIC TYPETAG=RichardsLensBoxProblem)
+
if(MPI_FOUND)
-dune_add_test(SOURCES test_richardslens_fv.cc
- NAME test_boxrichards_parallel
- COMPILE_DEFINITIONS TYPETAG=RichardsLensBoxProblem
+dune_add_test(NAME test_richards_box_parallel
+ TARGET test_richards_box
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
CMD_ARGS --script fuzzy
--files ${CMAKE_SOURCE_DIR}/test/references/richardslensbox-reference-parallel.vtu
- ${CMAKE_CURRENT_BINARY_DIR}/s0002-p0000-richardslensbox-00008.vtu
- --command "${MPIEXEC} -np 2 ${CMAKE_CURRENT_BINARY_DIR}/test_boxrichards test_richardslens.input -Problem.Name test_boxrichards")
+ ${CMAKE_CURRENT_BINARY_DIR}/s0002-p0000-test_richards_box_parallel-00008.vtu
+ --command "${MPIEXEC} -np 2 ${CMAKE_CURRENT_BINARY_DIR}/test_richards_box test_richardslens.input -Problem.Name test_richards_box_parallel")
+
+
+dune_add_test(NAME test_richards_tpfa_parallel
+ TARGET test_richards_tpfa
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/richardslenscc-reference-parallel.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/s0002-p0000-test_richards_tpfa_parallel-00007.vtu
+ --command "${MPIEXEC} -np 2 ${CMAKE_CURRENT_BINARY_DIR}/test_richards_tpfa test_richardslens.input -Problem.Name test_richards_tpfa_parallel")
endif()
# isothermal tests
-dune_add_test(SOURCES test_richardslens_fv.cc
- NAME test_boxrichards
- COMPILE_DEFINITIONS TYPETAG=RichardsLensBoxProblem
+dune_add_test(TARGET test_richards_box
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
CMD_ARGS --script fuzzy
--files ${CMAKE_SOURCE_DIR}/test/references/richardslensbox-reference.vtu
- ${CMAKE_CURRENT_BINARY_DIR}/test_boxrichards-00008.vtu
- --command "${CMAKE_CURRENT_BINARY_DIR}/test_boxrichards test_richardslens.input -Problem.Name test_boxrichards")
+ ${CMAKE_CURRENT_BINARY_DIR}/test_richards_box-00008.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_richards_box test_richardslens.input -Problem.Name test_richards_box")
-dune_add_test(SOURCES test_richardslens_fv.cc
- NAME test_ccrichards
- COMPILE_DEFINITIONS TYPETAG=RichardsLensCCProblem
+dune_add_test(TARGET test_richards_tpfa
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
CMD_ARGS --script fuzzy
--files ${CMAKE_SOURCE_DIR}/test/references/richardslenscc-reference.vtu
- ${CMAKE_CURRENT_BINARY_DIR}/test_ccrichards-00007.vtu
- --command "${CMAKE_CURRENT_BINARY_DIR}/test_ccrichards test_richardslens.input -Problem.Name test_ccrichards")
+ ${CMAKE_CURRENT_BINARY_DIR}/test_richards_tpfa-00007.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_richards_tpfa test_richardslens.input -Problem.Name test_richards_tpfa")
# comparison to analytical solution - only with cc
dune_add_test(SOURCES test_ccrichardsanalytical.cc
diff --git a/test/porousmediumflow/richards/implicit/test_richardslens.input b/test/porousmediumflow/richards/implicit/test_richardslens.input
index c602ef1cf9..e0ea32dc43 100644
--- a/test/porousmediumflow/richards/implicit/test_richardslens.input
+++ b/test/porousmediumflow/richards/implicit/test_richardslens.input
@@ -10,10 +10,6 @@ Cells = 24 16
Name = richardslens
EnableGravity = 1 # enable gravity
-[Implicit]
-EnablePartialReassemble = true # enable partial reassembly of the Jacobian matrix
-EnableJacobianRecycling = true #
-
[Newton]
TargetSteps = 18 # set the "desireable" number of newton iterations of a time step
EnableChop = true # chop for better convergence
diff --git a/test/porousmediumflow/richards/implicit/test_richardslens_fv.cc b/test/porousmediumflow/richards/implicit/test_richardslens_fv.cc
index dd3ebbdd69..dd81e7f93f 100644
--- a/test/porousmediumflow/richards/implicit/test_richardslens_fv.cc
+++ b/test/porousmediumflow/richards/implicit/test_richardslens_fv.cc
@@ -153,8 +153,8 @@ int main(int argc, char** argv) try
auto assembler = std::make_shared<Assembler>(problem, fvGridGeometry, gridVariables, timeLoop);
// the linear solver
- using LinearSolver = Dumux::ILU0BiCGSTABBackend<TypeTag>;
- auto linearSolver = std::make_shared<LinearSolver>();
+ using LinearSolver = Dumux::AMGBackend<TypeTag>;
+ auto linearSolver = std::make_shared<LinearSolver>(leafGridView, fvGridGeometry->dofMapper());
// the non-linear solver
using NewtonController = Dumux::RichardsNewtonController<TypeTag>;
diff --git a/test/references/richardslenscc-reference-parallel.vtu b/test/references/richardslenscc-reference-parallel.vtu
new file mode 100644
index 0000000000..e41093ec1e
--- /dev/null
+++ b/test/references/richardslenscc-reference-parallel.vtu
@@ -0,0 +1,407 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="192" NumberOfPoints="221">
+ <CellData Scalars="Sw">
+ <DataArray type="Float32" Name="Sw" NumberOfComponents="1" format="ascii">
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17 -4.36236e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 0.00714573 0.0116268 0.0116268 0.00714573
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 0.0033326 0.296918 0.339898 0.339898 0.296918
+ 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 9.35838e-17 2.85357e-06 0.0790732 0.814737 0.847266 0.847266 0.814737
+ </DataArray>
+ <DataArray type="Float32" Name="Sn" 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 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 0.992854 0.988373 0.988373 0.992854
+ 1 1 1 1 1 1 1 0.996667 0.703082 0.660102 0.660102 0.703082
+ 1 1 1 1 1 1 0.999997 0.920927 0.185263 0.152734 0.152734 0.185263
+ </DataArray>
+ <DataArray type="Float32" Name="pw" NumberOfComponents="1" format="ascii">
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+ 97470.4 97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+ 97470.4 97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+ 97470.4 97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97661.6 97781.4 97781.4 97661.6
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97559.5 99627.1 99646.8 99646.8 99627.1
+ 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.5 99308.3 99788.3 99798.9 99798.9 99788.3
+ </DataArray>
+ <DataArray type="Float32" Name="pn" NumberOfComponents="1" format="ascii">
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+ 2529.62 2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+ 2529.62 2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+ 2529.62 2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2338.44 2218.55 2218.55 2338.44
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2440.46 372.863 353.176 353.176 372.863
+ 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.54 691.737 211.699 201.127 201.127 211.699
+ </DataArray>
+ <DataArray type="Float32" Name="density" NumberOfComponents="1" format="ascii">
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ </DataArray>
+ <DataArray type="Float32" Name="mobility" 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 10.7357 17.6661 17.6661 10.7357
+ 0 0 0 0 0 0 0 0.0154609 407.511 473.279 473.279 407.511
+ </DataArray>
+ <DataArray type="Float32" Name="kr" 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.0107357 0.0176661 0.0176661 0.0107357
+ 0 0 0 0 0 0 0 1.54609e-05 0.407511 0.473279 0.473279 0.407511
+ </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
+ </DataArray>
+ <DataArray type="Float32" Name="pressure head" NumberOfComponents="1" format="ascii">
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149
+ -25.7861 -25.7861 -25.7861 -25.7861 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149
+ -25.7861 -25.7861 -25.7861 -25.7861 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149
+ -25.7861 -25.7861 -25.7861 -25.7861 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149 -219.149
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -23.8373 -22.6152 -22.6152 -23.8373
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -24.8773 -3.80085 -3.60016 -3.60016 -3.80085
+ -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7861 -25.7854 -7.05135 -2.15799 -2.05022 -2.05022 -2.15799
+ </DataArray>
+ <DataArray type="Float32" Name="water content" NumberOfComponents="1" format="ascii">
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17 -1.74495e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 0.00285829 0.00465071 0.00465071 0.00285829
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 0.00133304 0.118767 0.135959 0.135959 0.118767
+ 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 3.74335e-17 1.14143e-06 0.0316293 0.325895 0.338906 0.338906 0.325895
+ </DataArray>
+ <DataArray type="Float32" Name="phasePresence" NumberOfComponents="1" format="ascii">
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ 3 3 3 3 3 3 3 3 3 3 3 3
+ </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 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 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.25 0 0 0 0.25 0 0.25 0.25 0
+ 0.5 0 0 0.5 0.25 0 0.75 0 0 0.75 0.25 0
+ 1 0 0 1 0.25 0 1.25 0 0 1.25 0.25 0
+ 1.5 0 0 1.5 0.25 0 1.75 0 0 1.75 0.25 0
+ 2 0 0 2 0.25 0 2.25 0 0 2.25 0.25 0
+ 2.5 0 0 2.5 0.25 0 2.75 0 0 2.75 0.25 0
+ 3 0 0 3 0.25 0 0 0.5 0 0.25 0.5 0
+ 0.5 0.5 0 0.75 0.5 0 1 0.5 0 1.25 0.5 0
+ 1.5 0.5 0 1.75 0.5 0 2 0.5 0 2.25 0.5 0
+ 2.5 0.5 0 2.75 0.5 0 3 0.5 0 0 0.75 0
+ 0.25 0.75 0 0.5 0.75 0 0.75 0.75 0 1 0.75 0
+ 1.25 0.75 0 1.5 0.75 0 1.75 0.75 0 2 0.75 0
+ 2.25 0.75 0 2.5 0.75 0 2.75 0.75 0 3 0.75 0
+ 0 1 0 0.25 1 0 0.5 1 0 0.75 1 0
+ 1 1 0 1.25 1 0 1.5 1 0 1.75 1 0
+ 2 1 0 2.25 1 0 2.5 1 0 2.75 1 0
+ 3 1 0 0 1.25 0 0.25 1.25 0 0.5 1.25 0
+ 0.75 1.25 0 1 1.25 0 1.25 1.25 0 1.5 1.25 0
+ 1.75 1.25 0 2 1.25 0 2.25 1.25 0 2.5 1.25 0
+ 2.75 1.25 0 3 1.25 0 0 1.5 0 0.25 1.5 0
+ 0.5 1.5 0 0.75 1.5 0 1 1.5 0 1.25 1.5 0
+ 1.5 1.5 0 1.75 1.5 0 2 1.5 0 2.25 1.5 0
+ 2.5 1.5 0 2.75 1.5 0 3 1.5 0 0 1.75 0
+ 0.25 1.75 0 0.5 1.75 0 0.75 1.75 0 1 1.75 0
+ 1.25 1.75 0 1.5 1.75 0 1.75 1.75 0 2 1.75 0
+ 2.25 1.75 0 2.5 1.75 0 2.75 1.75 0 3 1.75 0
+ 0 2 0 0.25 2 0 0.5 2 0 0.75 2 0
+ 1 2 0 1.25 2 0 1.5 2 0 1.75 2 0
+ 2 2 0 2.25 2 0 2.5 2 0 2.75 2 0
+ 3 2 0 0 2.25 0 0.25 2.25 0 0.5 2.25 0
+ 0.75 2.25 0 1 2.25 0 1.25 2.25 0 1.5 2.25 0
+ 1.75 2.25 0 2 2.25 0 2.25 2.25 0 2.5 2.25 0
+ 2.75 2.25 0 3 2.25 0 0 2.5 0 0.25 2.5 0
+ 0.5 2.5 0 0.75 2.5 0 1 2.5 0 1.25 2.5 0
+ 1.5 2.5 0 1.75 2.5 0 2 2.5 0 2.25 2.5 0
+ 2.5 2.5 0 2.75 2.5 0 3 2.5 0 0 2.75 0
+ 0.25 2.75 0 0.5 2.75 0 0.75 2.75 0 1 2.75 0
+ 1.25 2.75 0 1.5 2.75 0 1.75 2.75 0 2 2.75 0
+ 2.25 2.75 0 2.5 2.75 0 2.75 2.75 0 3 2.75 0
+ 0 3 0 0.25 3 0 0.5 3 0 0.75 3 0
+ 1 3 0 1.25 3 0 1.5 3 0 1.75 3 0
+ 2 3 0 2.25 3 0 2.5 3 0 2.75 3 0
+ 3 3 0 0 3.25 0 0.25 3.25 0 0.5 3.25 0
+ 0.75 3.25 0 1 3.25 0 1.25 3.25 0 1.5 3.25 0
+ 1.75 3.25 0 2 3.25 0 2.25 3.25 0 2.5 3.25 0
+ 2.75 3.25 0 3 3.25 0 0 3.5 0 0.25 3.5 0
+ 0.5 3.5 0 0.75 3.5 0 1 3.5 0 1.25 3.5 0
+ 1.5 3.5 0 1.75 3.5 0 2 3.5 0 2.25 3.5 0
+ 2.5 3.5 0 2.75 3.5 0 3 3.5 0 0 3.75 0
+ 0.25 3.75 0 0.5 3.75 0 0.75 3.75 0 1 3.75 0
+ 1.25 3.75 0 1.5 3.75 0 1.75 3.75 0 2 3.75 0
+ 2.25 3.75 0 2.5 3.75 0 2.75 3.75 0 3 3.75 0
+ 0 4 0 0.25 4 0 0.5 4 0 0.75 4 0
+ 1 4 0 1.25 4 0 1.5 4 0 1.75 4 0
+ 2 4 0 2.25 4 0 2.5 4 0 2.75 4 0
+ 3 4 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 2 3 27 26 3 5 28 27 5 7 29 28
+ 7 9 30 29 9 11 31 30 11 13 32 31
+ 13 15 33 32 15 17 34 33 17 19 35 34
+ 19 21 36 35 21 23 37 36 23 25 38 37
+ 26 27 40 39 27 28 41 40 28 29 42 41
+ 29 30 43 42 30 31 44 43 31 32 45 44
+ 32 33 46 45 33 34 47 46 34 35 48 47
+ 35 36 49 48 36 37 50 49 37 38 51 50
+ 39 40 53 52 40 41 54 53 41 42 55 54
+ 42 43 56 55 43 44 57 56 44 45 58 57
+ 45 46 59 58 46 47 60 59 47 48 61 60
+ 48 49 62 61 49 50 63 62 50 51 64 63
+ 52 53 66 65 53 54 67 66 54 55 68 67
+ 55 56 69 68 56 57 70 69 57 58 71 70
+ 58 59 72 71 59 60 73 72 60 61 74 73
+ 61 62 75 74 62 63 76 75 63 64 77 76
+ 65 66 79 78 66 67 80 79 67 68 81 80
+ 68 69 82 81 69 70 83 82 70 71 84 83
+ 71 72 85 84 72 73 86 85 73 74 87 86
+ 74 75 88 87 75 76 89 88 76 77 90 89
+ 78 79 92 91 79 80 93 92 80 81 94 93
+ 81 82 95 94 82 83 96 95 83 84 97 96
+ 84 85 98 97 85 86 99 98 86 87 100 99
+ 87 88 101 100 88 89 102 101 89 90 103 102
+ 91 92 105 104 92 93 106 105 93 94 107 106
+ 94 95 108 107 95 96 109 108 96 97 110 109
+ 97 98 111 110 98 99 112 111 99 100 113 112
+ 100 101 114 113 101 102 115 114 102 103 116 115
+ 104 105 118 117 105 106 119 118 106 107 120 119
+ 107 108 121 120 108 109 122 121 109 110 123 122
+ 110 111 124 123 111 112 125 124 112 113 126 125
+ 113 114 127 126 114 115 128 127 115 116 129 128
+ 117 118 131 130 118 119 132 131 119 120 133 132
+ 120 121 134 133 121 122 135 134 122 123 136 135
+ 123 124 137 136 124 125 138 137 125 126 139 138
+ 126 127 140 139 127 128 141 140 128 129 142 141
+ 130 131 144 143 131 132 145 144 132 133 146 145
+ 133 134 147 146 134 135 148 147 135 136 149 148
+ 136 137 150 149 137 138 151 150 138 139 152 151
+ 139 140 153 152 140 141 154 153 141 142 155 154
+ 143 144 157 156 144 145 158 157 145 146 159 158
+ 146 147 160 159 147 148 161 160 148 149 162 161
+ 149 150 163 162 150 151 164 163 151 152 165 164
+ 152 153 166 165 153 154 167 166 154 155 168 167
+ 156 157 170 169 157 158 171 170 158 159 172 171
+ 159 160 173 172 160 161 174 173 161 162 175 174
+ 162 163 176 175 163 164 177 176 164 165 178 177
+ 165 166 179 178 166 167 180 179 167 168 181 180
+ 169 170 183 182 170 171 184 183 171 172 185 184
+ 172 173 186 185 173 174 187 186 174 175 188 187
+ 175 176 189 188 176 177 190 189 177 178 191 190
+ 178 179 192 191 179 180 193 192 180 181 194 193
+ 182 183 196 195 183 184 197 196 184 185 198 197
+ 185 186 199 198 186 187 200 199 187 188 201 200
+ 188 189 202 201 189 190 203 202 190 191 204 203
+ 191 192 205 204 192 193 206 205 193 194 207 206
+ 195 196 209 208 196 197 210 209 197 198 211 210
+ 198 199 212 211 199 200 213 212 200 201 214 213
+ 201 202 215 214 202 203 216 215 203 204 217 216
+ 204 205 218 217 205 206 219 218 206 207 220 219
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320 324 328 332 336
+ 340 344 348 352 356 360 364 368 372 376 380 384
+ 388 392 396 400 404 408 412 416 420 424 428 432
+ 436 440 444 448 452 456 460 464 468 472 476 480
+ 484 488 492 496 500 504 508 512 516 520 524 528
+ 532 536 540 544 548 552 556 560 564 568 572 576
+ 580 584 588 592 596 600 604 608 612 616 620 624
+ 628 632 636 640 644 648 652 656 660 664 668 672
+ 676 680 684 688 692 696 700 704 708 712 716 720
+ 724 728 732 736 740 744 748 752 756 760 764 768
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab