From ac38a707afb089382e135dc915d39d520a28d4b8 Mon Sep 17 00:00:00 2001
From: renchao_lu <renchao.lu@gmail.com>
Date: Sun, 2 May 2021 20:16:27 +0200
Subject: [PATCH] [PL/RCT] create media map.

---
 .../CreateRichardsComponentTransportProcess.cpp            | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp b/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
index ea8c8cafe6f..18b82257304 100644
--- a/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
+++ b/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
@@ -9,9 +9,11 @@
  */
 
 #include "CreateRichardsComponentTransportProcess.h"
-
 #include "CreatePorousMediaProperties.h"
+
 #include "MaterialLib/Fluid/FluidProperties/CreateFluidProperties.h"
+#include "MaterialLib/MPL/CreateMaterialSpatialDistributionMap.h"
+#include "MaterialLib/MPL/MaterialSpatialDistributionMap.h"
 #include "ParameterLib/ConstantParameter.h"
 #include "ParameterLib/Utils.h"
 #include "ProcessLib/Output/CreateSecondaryVariables.h"
@@ -151,6 +153,9 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
         std::copy_n(b.data(), b.size(), specific_body_force.data());
     }
 
+    auto media_map =
+        MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh);
+
     RichardsComponentTransportProcessData process_data{
         std::move(porous_media_properties),
         fluid_reference_density,
-- 
GitLab