diff --git a/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U.prj b/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U.prj
index 8aafcf08903dcc5d3079ccf115c15710d8de44dc..598f54be615ed70d83ed0639206613d267e5bf19 100644
--- a/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U.prj
+++ b/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U.prj
@@ -145,14 +145,14 @@
                     <type>AqueousLiquid</type>
                     <properties>
                         <property>
-                            <name>specific_heat_capacity</name>
+                            <name>phase_velocity</name>
                             <type>Constant</type>
-                            <value>4068</value>
+                            <value>0 0 0</value>
                         </property>
                         <property>
-                            <name>thermal_conductivity</name>
+                            <name>specific_heat_capacity</name>
                             <type>Constant</type>
-                            <value>0.1</value>
+                            <value>4068</value>
                         </property>
                         <property>
                             <name>density</name>
@@ -169,11 +169,6 @@
                             <type>Constant</type>
                             <value>1780</value>
                         </property>
-                        <property>
-                            <name>thermal_conductivity</name>
-                            <type>Constant</type>
-                            <value>2.4</value>
-                        </property>
                         <property>
                             <name>density</name>
                             <type>Constant</type>
@@ -189,11 +184,6 @@
                             <type>Constant</type>
                             <value>1000</value>
                         </property>
-                        <property>
-                            <name>thermal_conductivity</name>
-                            <type>Constant</type>
-                            <value>3.2</value>
-                        </property>
                         <property>
                             <name>density</name>
                             <type>Constant</type>
@@ -203,6 +193,16 @@
                 </phase>
             </phases>
             <properties>
+                <property>
+                    <name>porosity</name>
+                    <type>Constant</type>
+                    <value>0</value>
+                </property>
+                <property>
+                    <name>thermal_conductivity</name>
+                    <type>Constant</type>
+                    <value>2.4</value>
+                </property>
                 <property>
                     <name>thermal_longitudinal_dispersivity</name>
                     <type>Constant</type>
@@ -268,7 +268,6 @@
                 </variables>
         </output>
     </time_loop>
-
     <parameters>
         <parameter>
             <name>T0</name>
@@ -281,7 +280,6 @@
             <values>305.36 305.13 305.23 305.115</values>
         </parameter>
     </parameters>
-
     <process_variables>
         <process_variable>
             <name>temperature_soil</name>
diff --git a/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U_pcs_0_ts_10_t_7200.000000.vtu b/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U_pcs_0_ts_10_t_7200.000000.vtu
index 884ff5ff382e1859c176ce8947f085d484ab3f07..f46716f74dc471c20c1276c1c7f0f737576d05e8 100644
--- a/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U_pcs_0_ts_10_t_7200.000000.vtu
+++ b/Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U_pcs_0_ts_10_t_7200.000000.vtu
@@ -1,3 +1,3 @@
 version https://git-lfs.github.com/spec/v1
-oid sha256:4312a8c588a91c61251f8a16bf7539e01e6d61b849560ec516350b4924f29924
-size 220400
+oid sha256:b287cb24862ca7a802c67dfb071f26efc252baf1d5ca0bedacb979c28d1f6741
+size 220384
diff --git a/Tests/Data/Parabolic/T/3D_3BHEs_array/bcs_tespy.py b/Tests/Data/Parabolic/T/3D_3BHEs_array/bcs_tespy.py
index 556d3672cbe8cc81f3123d117f31ee3c2b3139ee..3efa453eaaa308bd4b0a13dd61fa62cc1576d9a4 100644
--- a/Tests/Data/Parabolic/T/3D_3BHEs_array/bcs_tespy.py
+++ b/Tests/Data/Parabolic/T/3D_3BHEs_array/bcs_tespy.py
@@ -39,7 +39,7 @@ def timerange(t):
 #month demand
 def consumer_demand(t):#dynamic thermal demand from consumer
     #thermal demand in each month, assumed specific heat extraction rate * lenth of BHE * number of BHE
-    month_demand = [-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3]
+    month_demand = [-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3]
     return month_demand[t -1]
 
 # %% dynamic hydraulic flow rate