diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0.prj
index 17e76a809cd9bbb471af054dd8e55e7462fc6bc6..e8285294f1079ab5b9569fa80049f5b51b24b499 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0_neumann.prj
index c5ebb1ddce41140ecba92b0fc3abb1776b5456a1..cffbfbff23f89d9966b35e36e7781921c7cf1d64 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e0_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e1_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e1_neumann.prj
index 0686aa5c821a90b30330530317c11f7f0fb34cbf..2920a2e4d47c1a1cbcaeb925c62dd856931dcb14 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e1_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e1_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e1.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2.prj
index 80ebfa42a62dfa438ece520dc501aed1a593003c..113876fb86bd3f02d5b6850bb26de134d66c830d 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e2.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2_neumann.prj
index 3548da9b4f87b4e06162c4aedacc8c7ab8193a63..19fd6bc3fca31e1fd25e056ffca28b6822bdf36c 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e2_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e2.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3.prj
index db8021bc80447e8fd60021635811f8b1fc4e4425..58d4e42e4733919e0f04b266f30c7da82364450d 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_bottom_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_bottom_neumann.prj
index d12f6f1062ee974df9d79a5cfcc7bb41699687a7..404af09936ef4438d1586133db80241c01ea8c7e 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_bottom_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_bottom_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann.prj
index 45290cf5fae87d39d70e1fea3ac8a63a8b6d7ea4..bbc6d04f78d29e78d1f8559bbdd14d49ac14f506 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_top_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_top_neumann.prj
index b713cddfffb176bac22480d417a902ec0479f49a..c73c5dc419f96f062ff4914302cef7bfd342e47e 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_top_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_top_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4.prj
index 2940cbc29e3ebf3285bd33c1c3d9fe14469ab771..542c50af78cb19313a601c4a9da27e8f8aa367e4 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4_neumann.prj
index 5e978eea2f2034f69b2c128ead7e5629f944d3f4..612992b55f0aabb891afb95fb303af6b76187607 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e4_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5.prj
index 2f3c592586a536eb9f10de47e6c6e11d2d6e680e..33b488d11008458435a0410a1236ee61fe35bb5b 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e5.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5_neumann.prj
index 85b0f5edae3c5b17e1e139c07aaf22cca65c8587..b6ede46aa84e47a12ac2c9d8fb7cbb003224d23d 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e5_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e5.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6.prj
index 2cb2317ad1ac08c62223556b038d2f0208e68bf2..0820b7ba2eaf1286e07f7f72ee3b1281a3c09f9c 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e6.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6_neumann.prj
index 1456da8f7fc04ca62532dbee7cb0eee7a0c279cf..55aa55c7214bb8541df7830c7aba39e945f47f36 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e6_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_1e6.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4.prj
index da10fdeda9ba1373e32c84e1ab432d762a3833a0..aeeb62d48f5888e94fd02a361c36738a960eb4bf 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_2e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4_neumann.prj
index 2fd2bbf27c3d85d9fb85be754249013ebb28954a..123f8eb030c377fb421e33ff79e93d25c31088b3 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e4_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_2e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4.prj
index 99fb415fc0b50ecb98ba2a973ccdd0878a0f8bd2..b1a4cf1b248c212d3c9adc2bc9877cded19a334d 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_3e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4_neumann.prj
index c3e4e789ed5067a7b34f94045ef32a1972c0b5cf..c63c11710128e9a03e3f6b19139944bd7795e773 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_3e4_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_3e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4.prj
index 0a6a0b261076df68826d21b72817c7413acdb8a9..e4800e451d022d418214e78bf285c22b2d419aba 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_4e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4_neumann.prj
index 4a3ee8d05ef61d01387825707f6992b7483de993..5407202b62937cc50d1b759bc5434af3b588e4f6 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_4e4_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_4e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4.prj
index 341d687af5323e5e08140c208fcca1be3a6b3ffd..961d246b7b9d93103f6c87a48e4e7a70aa59a7cc 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_5e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4_neumann.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4_neumann.prj
index 434864ce3e8425a328b2c8f2e59d85b75400c1bd..5122550ee835f25e0f70d7e4e351dd9ee037df4a 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4_neumann.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_5e4_neumann.prj
@@ -3,17 +3,7 @@
     <mesh>cube_1x1x1_hex_5e4.vtu</mesh>
     <geometry>cube_1x1x1.gml</geometry>
     <processes>
-        <process>
-            <name>SteadyStateDiffusion</name>
-            <type>STEADY_STATE_DIFFUSION</type>
-            <integration_order>2</integration_order>
-            <process_variables>
-                <process_variable>pressure</process_variable>
-            </process_variables>
-            <secondary_variables>
-                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
-            </secondary_variables>
-        </process>
+        <include file="SteadyStateDiffusion.xml" />
     </processes>
     <media>
         <medium id="0">