diff --git a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp b/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
index eb8de233e12a3da6c7f9db08142e7b5b3417cfd8..aba644015886ca73ea48984d348f597c01e0579c 100644
--- a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
+++ b/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
@@ -114,19 +114,19 @@ std::unique_ptr<Process> createComponentTransportProcess(
         //! \ogs_file_param_special{prj__processes__process__ComponentTransport__fluid_reference_density}
         "fluid_reference_density", parameters, 1);
     DBUG("Use '%s' as fluid_reference_density parameter.",
-         fluid_reference_density.name.c_str());
+         fluid_reference_density.name.c_str(), &mesh);
 
     // Parameter for the retardation factor.
     auto const& retardation_factor = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__ComponentTransport__retardation_factor}
-        "retardation_factor", parameters, 1);
+        "retardation_factor", parameters, 1, &mesh);
 
     // Parameter for the decay rate.
     auto const& decay_rate = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__ComponentTransport__decay_rate}
-        "decay_rate", parameters, 1);
+        "decay_rate", parameters, 1, &mesh);
 
     // Specific body force parameter.
     Eigen::VectorXd specific_body_force;
diff --git a/ProcessLib/GroundwaterFlow/CreateGroundwaterFlowProcess.cpp b/ProcessLib/GroundwaterFlow/CreateGroundwaterFlowProcess.cpp
index 772e003eb281dedcab53dd94ce65b87a38e834e7..816bd67db4c4524a32532bdc60dd8d2a9a3ea7a7 100644
--- a/ProcessLib/GroundwaterFlow/CreateGroundwaterFlowProcess.cpp
+++ b/ProcessLib/GroundwaterFlow/CreateGroundwaterFlowProcess.cpp
@@ -53,7 +53,8 @@ std::unique_ptr<Process> createGroundwaterFlowProcess(
     auto& hydraulic_conductivity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__GROUNDWATER_FLOW__hydraulic_conductivity}
-        "hydraulic_conductivity", parameters, 0 /*arbitrary many components*/);
+        "hydraulic_conductivity", parameters, 0 /*arbitrary many components*/,
+        &mesh);
 
     DBUG("Use '%s' as hydraulic conductivity parameter.",
          hydraulic_conductivity.name.c_str());
diff --git a/ProcessLib/HT/CreateHTProcess.cpp b/ProcessLib/HT/CreateHTProcess.cpp
index 69c8a75f39e077d964a821abe3b94e7f3296a87e..78b8117c712fb80299fe7f7ae6a76ff7889369a4 100644
--- a/ProcessLib/HT/CreateHTProcess.cpp
+++ b/ProcessLib/HT/CreateHTProcess.cpp
@@ -115,12 +115,12 @@ std::unique_ptr<Process> createHTProcess(
     {
         solid_thermal_expansion = &ParameterLib::findParameter<double>(
             //! \ogs_file_param_special{prj__processes__process__HT__solid_thermal_expansion__thermal_expansion}
-            *solid_config, "thermal_expansion", parameters, 1);
+            *solid_config, "thermal_expansion", parameters, 1, &mesh);
         DBUG("Use '%s' as solid thermal expansion.",
              solid_thermal_expansion->name.c_str());
         biot_constant = &ParameterLib::findParameter<double>(
             //! \ogs_file_param_special{prj__processes__process__HT__solid_thermal_expansion__biot_constant}
-            *solid_config, "biot_constant", parameters, 1);
+            *solid_config, "biot_constant", parameters, 1, &mesh);
         DBUG("Use '%s' as Biot's constant.", biot_constant->name.c_str());
     }
 
diff --git a/ProcessLib/HeatConduction/CreateHeatConductionProcess.cpp b/ProcessLib/HeatConduction/CreateHeatConductionProcess.cpp
index 69b0f8c5ec37408d701ffe45ad6a33dec10c7267..94c7520a31d6d2ad55beb1fad4df71111e0aeee7 100644
--- a/ProcessLib/HeatConduction/CreateHeatConductionProcess.cpp
+++ b/ProcessLib/HeatConduction/CreateHeatConductionProcess.cpp
@@ -49,7 +49,7 @@ std::unique_ptr<Process> createHeatConductionProcess(
     auto& thermal_conductivity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_CONDUCTION__thermal_conductivity}
-        "thermal_conductivity", parameters, 1);
+        "thermal_conductivity", parameters, 1, &mesh);
 
     DBUG("Use '%s' as thermal conductivity parameter.",
          thermal_conductivity.name.c_str());
@@ -58,7 +58,7 @@ std::unique_ptr<Process> createHeatConductionProcess(
     auto& heat_capacity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_CONDUCTION__heat_capacity}
-        "heat_capacity", parameters, 1);
+        "heat_capacity", parameters, 1, &mesh);
 
     DBUG("Use '%s' as heat capacity parameter.", heat_capacity.name.c_str());
 
@@ -66,7 +66,7 @@ std::unique_ptr<Process> createHeatConductionProcess(
     auto& density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_CONDUCTION__density}
-        "density", parameters, 1);
+        "density", parameters, 1, &mesh);
 
     DBUG("Use '%s' as density parameter.", density.name.c_str());
 
diff --git a/ProcessLib/HeatTransportBHE/CreateHeatTransportBHEProcess.cpp b/ProcessLib/HeatTransportBHE/CreateHeatTransportBHEProcess.cpp
index 2af6e6ab19758c8fbb84b6181440f1e40915d2c0..fddb992557bd38691323b2065adb0d33316445a9 100644
--- a/ProcessLib/HeatTransportBHE/CreateHeatTransportBHEProcess.cpp
+++ b/ProcessLib/HeatTransportBHE/CreateHeatTransportBHEProcess.cpp
@@ -90,7 +90,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& thermal_conductivity_solid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__thermal_conductivity_solid}
-        "thermal_conductivity_solid", parameters, 1);
+        "thermal_conductivity_solid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as solid phase thermal conductivity parameter.",
          thermal_conductivity_solid.name.c_str());
@@ -99,7 +99,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& thermal_conductivity_fluid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__thermal_conductivity_fluid}
-        "thermal_conductivity_fluid", parameters, 1);
+        "thermal_conductivity_fluid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as fluid phase thermal conductivity parameter.",
          thermal_conductivity_fluid.name.c_str());
@@ -108,7 +108,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& thermal_conductivity_gas = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__thermal_conductivity_gas}
-        "thermal_conductivity_gas", parameters, 1);
+        "thermal_conductivity_gas", parameters, 1, &mesh);
 
     DBUG("Use '%s' as gas phase thermal conductivity parameter.",
          thermal_conductivity_gas.name.c_str());
@@ -117,7 +117,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& heat_capacity_solid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__heat_capacity_solid}
-        "heat_capacity_solid", parameters, 1);
+        "heat_capacity_solid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as solid phase heat capacity parameter.",
          heat_capacity_solid.name.c_str());
@@ -126,7 +126,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& heat_capacity_fluid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__heat_capacity_fluid}
-        "heat_capacity_fluid", parameters, 1);
+        "heat_capacity_fluid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as fluid phase heat capacity parameter.",
          heat_capacity_fluid.name.c_str());
@@ -135,7 +135,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& heat_capacity_gas = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__heat_capacity_gas}
-        "heat_capacity_gas", parameters, 1);
+        "heat_capacity_gas", parameters, 1, &mesh);
 
     DBUG("Use '%s' as gas phase heat capacity parameter.",
          heat_capacity_gas.name.c_str());
@@ -144,7 +144,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& density_solid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__density_solid}
-        "density_solid", parameters, 1);
+        "density_solid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as solid phase density parameter.",
          density_solid.name.c_str());
@@ -153,7 +153,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& density_fluid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__density_fluid}
-        "density_fluid", parameters, 1);
+        "density_fluid", parameters, 1, &mesh);
 
     DBUG("Use '%s' as fluid phase density parameter.",
          density_fluid.name.c_str());
@@ -162,7 +162,7 @@ std::unique_ptr<Process> createHeatTransportBHEProcess(
     auto& density_gas = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HEAT_TRANSPORT_BHE__density_gas}
-        "density_gas", parameters, 1);
+        "density_gas", parameters, 1, &mesh);
 
     DBUG("Use '%s' as gas phase density parameter.", density_gas.name.c_str());
 
diff --git a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
index d42075c4fdb2bad4c707a9bc32514236562cd8dd..ee0c9f538c933615ba2e673e1ca80afd429470d2 100644
--- a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
+++ b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
@@ -111,7 +111,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& intrinsic_permeability = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__intrinsic_permeability}
-        "intrinsic_permeability", parameters, 1);
+        "intrinsic_permeability", parameters, 1, &mesh);
 
     DBUG("Use '%s' as intrinsic conductivity parameter.",
          intrinsic_permeability.name.c_str());
@@ -120,7 +120,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& specific_storage = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__specific_storage}
-        "specific_storage", parameters, 1);
+        "specific_storage", parameters, 1, &mesh);
 
     DBUG("Use '%s' as storage coefficient parameter.",
          specific_storage.name.c_str());
@@ -129,7 +129,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& fluid_viscosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__fluid_viscosity}
-        "fluid_viscosity", parameters, 1);
+        "fluid_viscosity", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid viscosity parameter.",
          fluid_viscosity.name.c_str());
 
@@ -137,14 +137,14 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& fluid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__fluid_density}
-        "fluid_density", parameters, 1);
+        "fluid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid density parameter.", fluid_density.name.c_str());
 
     // Biot coefficient
     auto& biot_coefficient = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__biot_coefficient}
-        "biot_coefficient", parameters, 1);
+        "biot_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as Biot coefficient parameter.",
          biot_coefficient.name.c_str());
 
@@ -152,14 +152,14 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& porosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__porosity}
-        "porosity", parameters, 1);
+        "porosity", parameters, 1, &mesh);
     DBUG("Use '%s' as porosity parameter.", porosity.name.c_str());
 
     // Solid density
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Specific body force
diff --git a/ProcessLib/LIE/HydroMechanics/CreateHydroMechanicsProcess.cpp b/ProcessLib/LIE/HydroMechanics/CreateHydroMechanicsProcess.cpp
index 3c5e291baf872d061c752081e0203d0f51fc5f85..68d6e22a4ed4e063358a170a8e69443155ba1401 100644
--- a/ProcessLib/LIE/HydroMechanics/CreateHydroMechanicsProcess.cpp
+++ b/ProcessLib/LIE/HydroMechanics/CreateHydroMechanicsProcess.cpp
@@ -138,7 +138,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& intrinsic_permeability = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__intrinsic_permeability}
-        "intrinsic_permeability", parameters, 1);
+        "intrinsic_permeability", parameters, 1, &mesh);
 
     DBUG("Use '%s' as intrinsic permeabiltiy parameter.",
          intrinsic_permeability.name.c_str());
@@ -147,7 +147,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& specific_storage = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__specific_storage}
-        "specific_storage", parameters, 1);
+        "specific_storage", parameters, 1, &mesh);
 
     DBUG("Use '%s' as specific storage parameter.",
          specific_storage.name.c_str());
@@ -156,7 +156,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& fluid_viscosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__fluid_viscosity}
-        "fluid_viscosity", parameters, 1);
+        "fluid_viscosity", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid viscosity parameter.",
          fluid_viscosity.name.c_str());
 
@@ -164,14 +164,14 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& fluid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__fluid_density}
-        "fluid_density", parameters, 1);
+        "fluid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid density parameter.", fluid_density.name.c_str());
 
     // Biot coefficient
     auto& biot_coefficient = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__biot_coefficient}
-        "biot_coefficient", parameters, 1);
+        "biot_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as Biot coefficient parameter.",
          biot_coefficient.name.c_str());
 
@@ -179,14 +179,14 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
     auto& porosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__porosity}
-        "porosity", parameters, 1);
+        "porosity", parameters, 1, &mesh);
     DBUG("Use '%s' as porosity parameter.", porosity.name.c_str());
 
     // Solid density
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Specific body force
@@ -264,13 +264,16 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
             fracture_properties_config.getConfigParameter<int>("material_id"),
             ParameterLib::findParameter<double>(
                 //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__fracture_properties__initial_aperture}
-                fracture_properties_config, "initial_aperture", parameters, 1),
+                fracture_properties_config, "initial_aperture", parameters, 1,
+                &mesh),
             ParameterLib::findParameter<double>(
                 //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__fracture_properties__specific_storage}
-                fracture_properties_config, "specific_storage", parameters, 1),
+                fracture_properties_config, "specific_storage", parameters, 1,
+                &mesh),
             ParameterLib::findParameter<double>(
                 //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__fracture_properties__biot_coefficient}
-                fracture_properties_config, "biot_coefficient", parameters, 1));
+                fracture_properties_config, "biot_coefficient", parameters, 1,
+                &mesh));
         if (frac_prop->aperture0.isTimeDependent())
         {
             OGS_FATAL(
@@ -285,7 +288,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__initial_effective_stress}
         "initial_effective_stress", parameters,
-        MathLib::KelvinVector::KelvinVectorDimensions<GlobalDim>::value);
+        MathLib::KelvinVector::KelvinVectorDimensions<GlobalDim>::value, &mesh);
     DBUG("Use '%s' as initial effective stress parameter.",
          initial_effective_stress.name.c_str());
 
@@ -294,7 +297,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
         double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS_WITH_LIE__initial_fracture_effective_stress}
-        "initial_fracture_effective_stress", parameters, GlobalDim);
+        "initial_fracture_effective_stress", parameters, GlobalDim, &mesh);
     DBUG("Use '%s' as initial fracture effective stress parameter.",
          initial_fracture_effective_stress.name.c_str());
 
diff --git a/ProcessLib/LIE/SmallDeformation/CreateSmallDeformationProcess.cpp b/ProcessLib/LIE/SmallDeformation/CreateSmallDeformationProcess.cpp
index 92fbfdac86d4ae824d2cd45b8bc98f1b879bc000..b8dff71f0157b5bbcbc1275359f3fe468e238b72 100644
--- a/ProcessLib/LIE/SmallDeformation/CreateSmallDeformationProcess.cpp
+++ b/ProcessLib/LIE/SmallDeformation/CreateSmallDeformationProcess.cpp
@@ -161,7 +161,8 @@ std::unique_ptr<Process> createSmallDeformationProcess(
             fracture_properties_config.getConfigParameter<int>("material_id"),
             ParameterLib::findParameter<double>(
                 //! \ogs_file_param_special{prj__processes__process__SMALL_DEFORMATION_WITH_LIE__fracture_properties__initial_aperture}
-                fracture_properties_config, "initial_aperture", parameters, 1));
+                fracture_properties_config, "initial_aperture", parameters, 1,
+                &mesh));
     }
 
     if (n_var_du < fracture_properties.size())
diff --git a/ProcessLib/ProcessVariable.cpp b/ProcessLib/ProcessVariable.cpp
index bb1ad5aafe16db70d3407ff775bc7e16b10a9702..d22d5187325769425eb65fd019a6b96f506b4e31 100644
--- a/ProcessLib/ProcessVariable.cpp
+++ b/ProcessLib/ProcessVariable.cpp
@@ -104,7 +104,7 @@ ProcessVariable::ProcessVariable(
       _initial_condition(ParameterLib::findParameter<double>(
           //! \ogs_file_param{prj__process_variables__process_variable__initial_condition}
           config.getConfigParameter<std::string>("initial_condition"),
-          parameters, _n_components))
+          parameters, _n_components, &mesh))
 {
     DBUG("Constructing process variable %s", _name.c_str());
 
diff --git a/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp b/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
index 4b41756284c00404fe8818dce927d60b018ee3bc..716a63eba490dd360abc0c1752447cee25731830 100644
--- a/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
+++ b/ProcessLib/RichardsComponentTransport/CreateRichardsComponentTransportProcess.cpp
@@ -87,7 +87,7 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
     auto& fluid_reference_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__fluid_reference_density}
-        "fluid_reference_density", parameters, 1);
+        "fluid_reference_density", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid_reference_density parameter.",
          fluid_reference_density.name.c_str());
 
@@ -96,7 +96,7 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
         double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__molecular_diffusion_coefficient}
-        "molecular_diffusion_coefficient", parameters, 1);
+        "molecular_diffusion_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as molecular diffusion coefficient parameter.",
          molecular_diffusion_coefficient.name.c_str());
 
@@ -105,7 +105,7 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
         double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__solute_dispersivity_longitudinal}
-        "solute_dispersivity_longitudinal", parameters, 1);
+        "solute_dispersivity_longitudinal", parameters, 1, &mesh);
     DBUG("Use '%s' as longitudinal solute dispersivity parameter.",
          solute_dispersivity_longitudinal.name.c_str());
 
@@ -114,7 +114,7 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
         double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__solute_dispersivity_transverse}
-        "solute_dispersivity_transverse", parameters, 1);
+        "solute_dispersivity_transverse", parameters, 1, &mesh);
     DBUG("Use '%s' as transverse solute dispersivity parameter.",
          solute_dispersivity_transverse.name.c_str());
 
@@ -122,13 +122,13 @@ std::unique_ptr<Process> createRichardsComponentTransportProcess(
     auto const& retardation_factor = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__retardation_factor}
-        "retardation_factor", parameters, 1);
+        "retardation_factor", parameters, 1, &mesh);
 
     // Parameter for the decay rate.
     auto const& decay_rate = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RichardsComponentTransport__decay_rate}
-        "decay_rate", parameters, 1);
+        "decay_rate", parameters, 1, &mesh);
 
     // Specific body force parameter.
     Eigen::VectorXd specific_body_force;
diff --git a/ProcessLib/RichardsFlow/CreateRichardsFlowProcess.cpp b/ProcessLib/RichardsFlow/CreateRichardsFlowProcess.cpp
index cf78a715d565d48a98cdfcf6ac0829f6e731d514..2c5fb72c9f530810e6cd5facc2edbd9a8542b06b 100644
--- a/ProcessLib/RichardsFlow/CreateRichardsFlowProcess.cpp
+++ b/ProcessLib/RichardsFlow/CreateRichardsFlowProcess.cpp
@@ -78,7 +78,7 @@ std::unique_ptr<Process> createRichardsFlowProcess(
     auto& temperature = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_FLOW__temperature}
-        "temperature", parameters, 1);
+        "temperature", parameters, 1, &mesh);
 
     //! \ogs_file_param{prj__processes__process__RICHARDS_FLOW__material_property}
     auto const& mat_config = config.getConfigSubtree("material_property");
diff --git a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
index 2ca3e94b47b9cea25578ec012892bdb3ccbee664..32602b0bc4fc634f79267579cb40f1c70a66a442 100644
--- a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
+++ b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
@@ -112,7 +112,7 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
     auto& fluid_bulk_modulus = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_MECHANICS__fluid_bulk_modulus}
-        "fluid_bulk_modulus", parameters, 1);
+        "fluid_bulk_modulus", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid bulk modulus parameter.",
          fluid_bulk_modulus.name.c_str());
 
@@ -120,7 +120,7 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
     auto& biot_coefficient = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_MECHANICS__biot_coefficient}
-        "biot_coefficient", parameters, 1);
+        "biot_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as Biot coefficient parameter.",
          biot_coefficient.name.c_str());
 
@@ -128,14 +128,14 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_MECHANICS__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Solid bulk modulus
     auto& solid_bulk_modulus = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_MECHANICS__solid_bulk_modulus}
-        "solid_bulk_modulus", parameters, 1);
+        "solid_bulk_modulus", parameters, 1, &mesh);
     DBUG("Use '%s' as solid bulk modulus parameter.",
          solid_bulk_modulus.name.c_str());
 
@@ -161,7 +161,7 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
     auto& temperature = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__RICHARDS_MECHANICS__temperature}
-        "temperature", parameters, 1);
+        "temperature", parameters, 1, &mesh);
 
     auto const& flow_material_config =
         //! \ogs_file_param{prj__processes__process__RICHARDS_MECHANICS__material_property}
diff --git a/ProcessLib/SmallDeformation/CreateSmallDeformationProcess.cpp b/ProcessLib/SmallDeformation/CreateSmallDeformationProcess.cpp
index 18de53ae860a67546a7209b76664f2e5b78ec0d6..f5c9b068c50760ea05f28ea9a2b16b8fa16f0c0c 100644
--- a/ProcessLib/SmallDeformation/CreateSmallDeformationProcess.cpp
+++ b/ProcessLib/SmallDeformation/CreateSmallDeformationProcess.cpp
@@ -73,7 +73,7 @@ std::unique_ptr<Process> createSmallDeformationProcess(
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__SMALL_DEFORMATION__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Specific body force
diff --git a/ProcessLib/SmallDeformationNonlocal/CreateSmallDeformationNonlocalProcess.cpp b/ProcessLib/SmallDeformationNonlocal/CreateSmallDeformationNonlocalProcess.cpp
index e42f887d54b6b35c09fbe830f8807818ac90e371..23e7b42a1712d4d024534adb7656b19c585f3c6c 100644
--- a/ProcessLib/SmallDeformationNonlocal/CreateSmallDeformationNonlocalProcess.cpp
+++ b/ProcessLib/SmallDeformationNonlocal/CreateSmallDeformationNonlocalProcess.cpp
@@ -75,7 +75,7 @@ std::unique_ptr<Process> createSmallDeformationNonlocalProcess(
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__SMALL_DEFORMATION_NONLOCAL__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Specific body force
diff --git a/ProcessLib/ThermalTwoPhaseFlowWithPP/CreateThermalTwoPhaseFlowWithPPProcess.cpp b/ProcessLib/ThermalTwoPhaseFlowWithPP/CreateThermalTwoPhaseFlowWithPPProcess.cpp
index 2b16dd30f0b91badae698904e6cf38a1af798594..dca8ef5122903c556930c9df7ef5ed51176ad7d5 100644
--- a/ProcessLib/ThermalTwoPhaseFlowWithPP/CreateThermalTwoPhaseFlowWithPPProcess.cpp
+++ b/ProcessLib/ThermalTwoPhaseFlowWithPP/CreateThermalTwoPhaseFlowWithPPProcess.cpp
@@ -81,25 +81,25 @@ std::unique_ptr<Process> createThermalTwoPhaseFlowWithPPProcess(
     auto& diff_coeff_b = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_THERMAL__diffusion_coeff_component_b}
-        "diffusion_coeff_component_b", parameters, 1);
+        "diffusion_coeff_component_b", parameters, 1, &mesh);
     auto& diff_coeff_a = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_THERMAL__diffusion_coeff_component_a}
-        "diffusion_coeff_component_a", parameters, 1);
+        "diffusion_coeff_component_a", parameters, 1, &mesh);
 
     // Parameter for the density of the solid.
 
     auto& density_solid = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_THERMAL__density_solid}
-        "density_solid", parameters, 1);
+        "density_solid", parameters, 1, &mesh);
     DBUG("Use '%s' as density_solid parameter.", density_solid.name.c_str());
 
     // Parameter for the latent heat of evaporation.
     auto& latent_heat_evaporation = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_THERMAL__latent_heat_evaporation}
-        "latent_heat_evaporation", parameters, 1);
+        "latent_heat_evaporation", parameters, 1, &mesh);
     DBUG("Use '%s' as latent_heat_evaporation parameter.",
          latent_heat_evaporation.name.c_str());
 
diff --git a/ProcessLib/ThermoHydroMechanics/CreateThermoHydroMechanicsProcess.cpp b/ProcessLib/ThermoHydroMechanics/CreateThermoHydroMechanicsProcess.cpp
index f319d69ca25988543fac89e3f24b68261ddb2a7f..2b30d9659f4c888839ebc46271192966993a46cb 100644
--- a/ProcessLib/ThermoHydroMechanics/CreateThermoHydroMechanicsProcess.cpp
+++ b/ProcessLib/ThermoHydroMechanics/CreateThermoHydroMechanicsProcess.cpp
@@ -130,7 +130,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& intrinsic_permeability = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__intrinsic_permeability}
-        "intrinsic_permeability", parameters, 1);
+        "intrinsic_permeability", parameters, 1, &mesh);
 
     DBUG("Use '%s' as intrinsic conductivity parameter.",
          intrinsic_permeability.name.c_str());
@@ -139,7 +139,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& specific_storage = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__specific_storage}
-        "specific_storage", parameters, 1);
+        "specific_storage", parameters, 1, &mesh);
 
     DBUG("Use '%s' as storage coefficient parameter.",
          specific_storage.name.c_str());
@@ -148,7 +148,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& fluid_viscosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__fluid_viscosity}
-        "fluid_viscosity", parameters, 1);
+        "fluid_viscosity", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid viscosity parameter.",
          fluid_viscosity.name.c_str());
 
@@ -156,14 +156,14 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& fluid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__fluid_density}
-        "fluid_density", parameters, 1);
+        "fluid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid density parameter.", fluid_density.name.c_str());
 
     // Biot coefficient
     auto& biot_coefficient = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__biot_coefficient}
-        "biot_coefficient", parameters, 1);
+        "biot_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as Biot coefficient parameter.",
          biot_coefficient.name.c_str());
 
@@ -171,14 +171,14 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& porosity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__porosity}
-        "porosity", parameters, 1);
+        "porosity", parameters, 1, &mesh);
     DBUG("Use '%s' as porosity parameter.", porosity.name.c_str());
 
     // Solid density
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // linear thermal expansion coefficient for solid
@@ -186,7 +186,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
         ParameterLib::findParameter<double>(
             config,
             //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__solid_linear_thermal_expansion_coefficient}
-            "solid_linear_thermal_expansion_coefficient", parameters, 1);
+            "solid_linear_thermal_expansion_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as solid linear thermal expansion coefficient parameter.",
          solid_linear_thermal_expansion_coefficient.name.c_str());
 
@@ -195,7 +195,8 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
         ParameterLib::findParameter<double>(
             config,
             //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__fluid_volumetric_thermal_expansion_coefficient}
-            "fluid_volumetric_thermal_expansion_coefficient", parameters, 1);
+            "fluid_volumetric_thermal_expansion_coefficient", parameters, 1,
+            &mesh);
     DBUG(
         "Use '%s' as fluid volumetric thermal expansion coefficient "
         "parameter.",
@@ -205,7 +206,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& solid_specific_heat_capacity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__solid_specific_heat_capacity}
-        "solid_specific_heat_capacity", parameters, 1);
+        "solid_specific_heat_capacity", parameters, 1, &mesh);
     DBUG("Use '%s' as solid specific heat capacity parameter.",
          solid_specific_heat_capacity.name.c_str());
 
@@ -213,7 +214,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& fluid_specific_heat_capacity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__fluid_specific_heat_capacity}
-        "fluid_specific_heat_capacity", parameters, 1);
+        "fluid_specific_heat_capacity", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid specific heat capacity parameter.",
          fluid_specific_heat_capacity.name.c_str());
 
@@ -221,7 +222,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& solid_thermal_conductivity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__solid_thermal_conductivity}
-        "solid_thermal_conductivity", parameters, 1);
+        "solid_thermal_conductivity", parameters, 1, &mesh);
     DBUG("Use '%s' as solid thermal conductivity parameter.",
          solid_thermal_conductivity.name.c_str());
 
@@ -229,7 +230,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& fluid_thermal_conductivity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__fluid_thermal_conductivity}
-        "fluid_thermal_conductivity", parameters, 1);
+        "fluid_thermal_conductivity", parameters, 1, &mesh);
     DBUG("Use '%s' as fluid thermal conductivity parameter.",
          fluid_thermal_conductivity.name.c_str());
 
@@ -237,7 +238,7 @@ std::unique_ptr<Process> createThermoHydroMechanicsProcess(
     auto& reference_temperature = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_HYDRO_MECHANICS__reference_temperature}
-        "reference_temperature", parameters, 1);
+        "reference_temperature", parameters, 1, &mesh);
     DBUG("Use '%s' as reference temperature parameter.",
          reference_temperature.name.c_str());
 
diff --git a/ProcessLib/ThermoMechanicalPhaseField/CreateThermoMechanicalPhaseFieldProcess.cpp b/ProcessLib/ThermoMechanicalPhaseField/CreateThermoMechanicalPhaseFieldProcess.cpp
index 1d4b4ed53da01f1f46745ccc60a0c3c50361fa02..ab90b585a1c570ab088b82cc54c0eec2a2c2a278 100644
--- a/ProcessLib/ThermoMechanicalPhaseField/CreateThermoMechanicalPhaseFieldProcess.cpp
+++ b/ProcessLib/ThermoMechanicalPhaseField/CreateThermoMechanicalPhaseFieldProcess.cpp
@@ -128,35 +128,35 @@ std::unique_ptr<Process> createThermoMechanicalPhaseFieldProcess(
     auto& residual_stiffness = ParameterLib::findParameter<double>(
         phasefield_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__phasefield_parameters__residual_stiffness}
-        "residual_stiffness", parameters, 1);
+        "residual_stiffness", parameters, 1, &mesh);
     DBUG("Use '%s' as residual stiffness.", residual_stiffness.name.c_str());
 
     // Crack resistance
     auto& crack_resistance = ParameterLib::findParameter<double>(
         phasefield_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__phasefield_parameters__crack_resistance}
-        "crack_resistance", parameters, 1);
+        "crack_resistance", parameters, 1, &mesh);
     DBUG("Use '%s' as crack resistance.", crack_resistance.name.c_str());
 
     // Crack length scale
     auto& crack_length_scale = ParameterLib::findParameter<double>(
         phasefield_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__phasefield_parameters__crack_length_scale}
-        "crack_length_scale", parameters, 1);
+        "crack_length_scale", parameters, 1, &mesh);
     DBUG("Use '%s' as crack length scale.", crack_length_scale.name.c_str());
 
     // Kinetic coefficient
     auto& kinetic_coefficient = ParameterLib::findParameter<double>(
         phasefield_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__phasefield_parameters__kinetic_coefficient}
-        "kinetic_coefficient", parameters, 1);
+        "kinetic_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as kinetic coefficient.", kinetic_coefficient.name.c_str());
 
     // Solid density
     auto& solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__reference_solid_density}
-        "solid_density", parameters, 1);
+        "solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.", solid_density.name.c_str());
 
     // Linear thermal expansion coefficient
@@ -164,7 +164,7 @@ std::unique_ptr<Process> createThermoMechanicalPhaseFieldProcess(
         double>(
         thermal_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__thermal_parameters__linear_thermal_expansion_coefficient}
-        "linear_thermal_expansion_coefficient", parameters, 1);
+        "linear_thermal_expansion_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as linear thermal expansion coefficient.",
          linear_thermal_expansion_coefficient.name.c_str());
 
@@ -172,7 +172,7 @@ std::unique_ptr<Process> createThermoMechanicalPhaseFieldProcess(
     auto& specific_heat_capacity = ParameterLib::findParameter<double>(
         thermal_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__thermal_parameters__specific_heat_capacity}
-        "specific_heat_capacity", parameters, 1);
+        "specific_heat_capacity", parameters, 1, &mesh);
     DBUG("Use '%s' as specific heat capacity.",
          specific_heat_capacity.name.c_str());
 
@@ -180,14 +180,14 @@ std::unique_ptr<Process> createThermoMechanicalPhaseFieldProcess(
     auto& thermal_conductivity = ParameterLib::findParameter<double>(
         thermal_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__thermal_parameters__thermal_conductivity}
-        "thermal_conductivity", parameters, 1);
+        "thermal_conductivity", parameters, 1, &mesh);
     DBUG("Use '%s' as thermal conductivity parameter.",
          thermal_conductivity.name.c_str());
     // Residual thermal conductivity
     auto& residual_thermal_conductivity = ParameterLib::findParameter<double>(
         thermal_parameters_config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICAL_PHASE_FIELD__thermal_parameters__residual_thermal_conductivity}
-        "residual_thermal_conductivity", parameters, 1);
+        "residual_thermal_conductivity", parameters, 1, &mesh);
     DBUG("Use '%s' as residual thermal conductivity parameter.",
          residual_thermal_conductivity.name.c_str());
     // Reference temperature
diff --git a/ProcessLib/ThermoMechanics/CreateThermoMechanicsProcess.cpp b/ProcessLib/ThermoMechanics/CreateThermoMechanicsProcess.cpp
index ca9bcd0565b9c5a02d273a6bb729c00f8b640748..28125527442939c424a92557596461d8a2d51980 100644
--- a/ProcessLib/ThermoMechanics/CreateThermoMechanicsProcess.cpp
+++ b/ProcessLib/ThermoMechanics/CreateThermoMechanicsProcess.cpp
@@ -113,7 +113,7 @@ std::unique_ptr<Process> createThermoMechanicsProcess(
     auto& reference_solid_density = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICS__reference_solid_density}
-        "reference_solid_density", parameters, 1);
+        "reference_solid_density", parameters, 1, &mesh);
     DBUG("Use '%s' as solid density parameter.",
          reference_solid_density.name.c_str());
 
@@ -122,21 +122,21 @@ std::unique_ptr<Process> createThermoMechanicsProcess(
         double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICS__linear_thermal_expansion_coefficient}
-        "linear_thermal_expansion_coefficient", parameters, 1);
+        "linear_thermal_expansion_coefficient", parameters, 1, &mesh);
     DBUG("Use '%s' as linear thermal expansion coefficient.",
          linear_thermal_expansion_coefficient.name.c_str());
     // Specific heat capacity
     auto& specific_heat_capacity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICS__specific_heat_capacity}
-        "specific_heat_capacity", parameters, 1);
+        "specific_heat_capacity", parameters, 1, &mesh);
     DBUG("Use '%s' as specific heat capacity parameter.",
          specific_heat_capacity.name.c_str());
     // Thermal conductivity // TODO To be changed as tensor input.
     auto& thermal_conductivity = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__THERMO_MECHANICS__thermal_conductivity}
-        "thermal_conductivity", parameters, 1);
+        "thermal_conductivity", parameters, 1, &mesh);
     DBUG("Use '%s' as thermal conductivity parameter.",
          thermal_conductivity.name.c_str());
 
diff --git a/ProcessLib/TwoPhaseFlowWithPP/CreateTwoPhaseFlowWithPPProcess.cpp b/ProcessLib/TwoPhaseFlowWithPP/CreateTwoPhaseFlowWithPPProcess.cpp
index 5b9d540019f168c9059e5e91b548dd26801fa24c..bf80c7a301a029ca702b2dfc222476752a4981c3 100644
--- a/ProcessLib/TwoPhaseFlowWithPP/CreateTwoPhaseFlowWithPPProcess.cpp
+++ b/ProcessLib/TwoPhaseFlowWithPP/CreateTwoPhaseFlowWithPPProcess.cpp
@@ -76,7 +76,7 @@ std::unique_ptr<Process> createTwoPhaseFlowWithPPProcess(
     auto& temperature = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_PP__temperature}
-        "temperature", parameters, 1);
+        "temperature", parameters, 1, &mesh);
 
     //! \ogs_file_param{prj__processes__process__TWOPHASE_FLOW_PP__material_property}
     auto const& mat_config = config.getConfigSubtree("material_property");
diff --git a/ProcessLib/TwoPhaseFlowWithPrho/CreateTwoPhaseFlowWithPrhoProcess.cpp b/ProcessLib/TwoPhaseFlowWithPrho/CreateTwoPhaseFlowWithPrhoProcess.cpp
index c738f92a20ebe5508a32957d0136c1d958b22d2b..104272d7d461cf1cfca00b0e3c29d7695854bb01 100644
--- a/ProcessLib/TwoPhaseFlowWithPrho/CreateTwoPhaseFlowWithPrhoProcess.cpp
+++ b/ProcessLib/TwoPhaseFlowWithPrho/CreateTwoPhaseFlowWithPrhoProcess.cpp
@@ -76,15 +76,15 @@ std::unique_ptr<Process> createTwoPhaseFlowWithPrhoProcess(
     auto& diff_coeff_b = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_PRHO__diffusion_coeff_component_b}
-        "diffusion_coeff_component_b", parameters, 1);
+        "diffusion_coeff_component_b", parameters, 1, &mesh);
     auto& diff_coeff_a = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_PRHO__diffusion_coeff_component_a}
-        "diffusion_coeff_component_a", parameters, 1);
+        "diffusion_coeff_component_a", parameters, 1, &mesh);
     auto& temperature = ParameterLib::findParameter<double>(
         config,
         //! \ogs_file_param_special{prj__processes__process__TWOPHASE_FLOW_PRHO__temperature}
-        "temperature", parameters, 1);
+        "temperature", parameters, 1, &mesh);
 
     //! \ogs_file_param{prj__processes__process__TWOPHASE_FLOW_PRHO__material_property}
     auto const& mat_config = config.getConfigSubtree("material_property");