diff --git a/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj b/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
index 531bc745120ca7b4dd58d49d172fab282ad130ab..848ec07f4a83cabc30a5f2e7458440834cff0ea4 100644
--- a/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/BuildupTest/buildup_test.prj b/Tests/Data/Parabolic/LiquidFlow/BuildupTest/buildup_test.prj
index 19bd14f34d17ed2980218d16eb4778a2fe395664..0d42f1f34049deb274d085bd480042b6f1454d5f 100644
--- a/Tests/Data/Parabolic/LiquidFlow/BuildupTest/buildup_test.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/BuildupTest/buildup_test.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/DrainageExcavation/drainage_LiquidFlow.prj b/Tests/Data/Parabolic/LiquidFlow/DrainageExcavation/drainage_LiquidFlow.prj
index db0fc614353d377dfb0b4925f9acc28c0caafcf4..7f3370a3c77cd205fede7e8d85b16448cb28a4ca 100644
--- a/Tests/Data/Parabolic/LiquidFlow/DrainageExcavation/drainage_LiquidFlow.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/DrainageExcavation/drainage_LiquidFlow.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
index f338f4b169b3e95c29ebb14c384c7cc7706093f5..563da1555a10636e2cd89de6e1cf95b0cafad0f8 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>square_1x1_tri_1.8e1_complete_boundary</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
index 21780ed544f8db644cf8d2330c1189d2732205ae..0cb01ecdb015d956d662b649985515cd7303c698 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>square_1x1_quad_1e1_complete_surface</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
index 6e00739b02d31c3d950894093e9bcc9e72aa4bc1..829af844a0e7cf24e15672fedfae490996b5985e 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>cuboid_1x1x1_hex_27_entire_boundary</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
index 3f631010dab581ffc58587d84823aa4568c7dd87..39ca442dcc5da33c2e661601f03f7453356f21b9 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>cuboid_1x1x1_hex_27_entire_boundary</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux-->
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
index b25d23ff7cc53ab27425b9070a12f6304fdc9cea..fe502413ab96cb9fb39903bd55b44a4c9141a366 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>cuboid_1x1x1_pyramid_6000_entire_boundary</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
index 914d43fb9e55570cb20b94494d93294843935dd3..d926872288a92668dd2db32ff0675a9635663cec 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
@@ -11,12 +11,6 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
@@ -27,6 +21,7 @@
<mesh>cube_1x1x1_hex_1e3_complete_surface</mesh>
<property_name>specific_flux</property_name>
</calculatesurfaceflux>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/GravityDriven/gravity_driven.prj b/Tests/Data/Parabolic/LiquidFlow/GravityDriven/gravity_driven.prj
index 6e471d2a40d137980e44f4034c534988917fe488..952e3c9adde588bfee9b98a061e28e5017a2e0eb 100644
--- a/Tests/Data/Parabolic/LiquidFlow/GravityDriven/gravity_driven.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/GravityDriven/gravity_driven.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/anisotropic_gravity_driven3D.prj b/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/anisotropic_gravity_driven3D.prj
index c2f5746b9e42911ba0237079add46bdfb8773a33..8317d6440621c5d913a9022ee5b90fcf6ee3d27d 100644
--- a/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/anisotropic_gravity_driven3D.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/anisotropic_gravity_driven3D.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>2</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/isotropic_gravity_driven3D.prj b/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/isotropic_gravity_driven3D.prj
index be503cd97b121a612951fde108427dd175a5bcb1..2384633e65a4bff9cf89aff61824b13853d8b0bc 100644
--- a/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/isotropic_gravity_driven3D.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/GravityDriven3D/isotropic_gravity_driven3D.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>2</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/LineDirichletNeumannBC/line_dirichlet_neumannBC.prj b/Tests/Data/Parabolic/LiquidFlow/LineDirichletNeumannBC/line_dirichlet_neumannBC.prj
index 1197b2b20abf54472f8bacb098582c4fa1b78500..c7ff6cced8d843fbd535a2628e9492d23c676fd2 100644
--- a/Tests/Data/Parabolic/LiquidFlow/LineDirichletNeumannBC/line_dirichlet_neumannBC.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/LineDirichletNeumannBC/line_dirichlet_neumannBC.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g >= 0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/PressureBCatCornerOfAnisotropicSquare/pressureBC_at_corner_of_anisotropic_square.prj b/Tests/Data/Parabolic/LiquidFlow/PressureBCatCornerOfAnisotropicSquare/pressureBC_at_corner_of_anisotropic_square.prj
index 205fe9f9bb583ac7fd7f0ecf85cc9af95b74a344..7621e960f50154c3549b51443c3509b5190057c0 100644
--- a/Tests/Data/Parabolic/LiquidFlow/PressureBCatCornerOfAnisotropicSquare/pressureBC_at_corner_of_anisotropic_square.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/PressureBCatCornerOfAnisotropicSquare/pressureBC_at_corner_of_anisotropic_square.prj
@@ -7,18 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest.prj
index 9355f400f039fca03e28eae826327c29d22545a4..79801a971b3ed3433cfed94573121b429b1fc089 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest.prj
@@ -10,16 +10,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>1</axis_id>
- <g>0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_1.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_1.prj
index 7c60c031c65960c75cf68c68732465dd91f937cd..4cf18fbd416876784afaffd10471b320f07a776b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_1.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_1.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_2.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_2.prj
index 2eedeb0a90c65dc6ead02e9a673b6e6eae61f03d..07fcf787913d68ef398e684be984308443f658fd 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_2.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_2.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_3.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_3.prj
index 11076532933e6cd2116145488da6ddf4cdbe5fc1..4d4f6f9dcde5df5f903b23db95bc76441d27145c 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_3.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_3.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>0</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>-0. 0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
index 7d7e36f1bf72b246e1b2b78efa9d0f9d50cb3425..2eed785076be2f34dd4dcf08428b187615345019 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
@@ -12,16 +12,13 @@ boundary meshes will be written as well -->
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>1</axis_id>
- <g>0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/constraint_bc_1e3.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/constraint_bc_1e3.prj
index 183b01006574c33fd9e411b7426b4b0835496bb8..9d65144a5269b4b0e7b1e82e4cb716a7fcabf5f3 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/constraint_bc_1e3.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/constraint_bc_1e3.prj
@@ -19,18 +19,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>2</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>p</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousBoundaryConditions/TimeDependentHeterogeneousBoundaryConditions.prj b/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousBoundaryConditions/TimeDependentHeterogeneousBoundaryConditions.prj
index d800527939fedd165a7b3b8eef61d0f77ea441bc..21d940f09d4d71cec2c52737f5962975444155d7 100644
--- a/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousBoundaryConditions/TimeDependentHeterogeneousBoundaryConditions.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousBoundaryConditions/TimeDependentHeterogeneousBoundaryConditions.prj
@@ -10,18 +10,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousSourceTerm/TimeDependentHeterogeneousSourceTerm.prj b/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousSourceTerm/TimeDependentHeterogeneousSourceTerm.prj
index e6c5bcd0f9b0bdc14adbfdda922bd66b8cec8aed..154395e3256eda68f5e01a270e31f4d7f53c36c8 100644
--- a/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousSourceTerm/TimeDependentHeterogeneousSourceTerm.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/TimeDependentHeterogeneousSourceTerm/TimeDependentHeterogeneousSourceTerm.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/TimeIntervalDirichletBC/TimeIntervalDirichletBC.prj b/Tests/Data/Parabolic/LiquidFlow/TimeIntervalDirichletBC/TimeIntervalDirichletBC.prj
index 6a1d63a7a95b994ba1cb0f429dcdd90397694341..eaac7112a2d40e6e071af09387a121b874c124f4 100644
--- a/Tests/Data/Parabolic/LiquidFlow/TimeIntervalDirichletBC/TimeIntervalDirichletBC.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/TimeIntervalDirichletBC/TimeIntervalDirichletBC.prj
@@ -10,18 +10,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC/TestSet_01.prj b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC/TestSet_01.prj
index 0fee8218f8a99589b3c19c9ef9e962900fcfe7d2..595a61945dbeafecc2d948f60474f5e72a172b49 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC/TestSet_01.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC/TestSet_01.prj
@@ -10,18 +10,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE/TestSet_01.prj b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE/TestSet_01.prj
index 8b25a9a5417c3655d15b61edadc6c7968d1581f7..b7f43e9b898a4bfcf0b15bd1389c329504050d9b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE/TestSet_01.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE/TestSet_01.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE2/TestSet_01.prj b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE2/TestSet_01.prj
index 1f37ce38d44dd7d4acb824c8a6542b17a252ac2c..1ebd98e0c29b76698c27a59746c13e5c7ef8186b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE2/TestSet_01.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_RECHARGE2/TestSet_01.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_STORAGE/TestSet_01.prj b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_STORAGE/TestSet_01.prj
index fbdb7a51a4f2917f30219f75764aabe4afe726ee..9525087438ec430bb474dbfbd46d11b62efecb60 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_STORAGE/TestSet_01.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Unconfined_Aquifer/BC_BC_STORAGE/TestSet_01.prj
@@ -11,18 +11,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>0.</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0</specific_body_force>
</process>
</processes>
<time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsource/h1_1Dsource.prj b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsource/h1_1Dsource.prj
index fa63d4a0d8074fd204c6cc0a578e74fcd648108d..7070400d05b5afa1b347322f872820bdab7010b9 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsource/h1_1Dsource.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsource/h1_1Dsource.prj
@@ -7,16 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>2</axis_id>
- <g> 0.0 </g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="q"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsteady/h1_1Dsteady.prj b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsteady/h1_1Dsteady.prj
index 1e93134849ea446b234ea389b55e98ea30826f8a..3fc1b5ebeb6df8f8ae3cb612d40c50c4e2449a1b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsteady/h1_1Dsteady.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_1Dsteady/h1_1Dsteady.prj
@@ -7,16 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>2</axis_id>
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="q"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_3Dhydstat/h1_3Dhydstat.prj b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_3Dhydstat/h1_3Dhydstat.prj
index ebc6e6c09da5315a401baf994db1e8cfb10ee681..77fd57b86aa207a9b9dd657c9577cc5498d0a3e0 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Verification/h1_3Dhydstat/h1_3Dhydstat.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Verification/h1_3Dhydstat/h1_3Dhydstat.prj
@@ -7,16 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>2</axis_id>
- <g>9.81</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="q"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 -9.81</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D1bt/h2_1D1bt.prj b/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D1bt/h2_1D1bt.prj
index 72e76c205d09017801b4fd50b3d62f93f8ea05a4..7f1680f8d953dc398a83e937ab999d1e2f9313dd 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D1bt/h2_1D1bt.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D1bt/h2_1D1bt.prj
@@ -7,16 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>2</axis_id>
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="q"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 0.0</specific_body_force>
</process>
</processes>
<media>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D2bt/h2_1D2bt.prj b/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D2bt/h2_1D2bt.prj
index 1c98959bfa8239478b4ef9ba845fad5cc7af1191..b6f2bb59694438ba7a8311e039dd58580965068b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D2bt/h2_1D2bt.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Verification/h2_1D2bt/h2_1D2bt.prj
@@ -7,16 +7,13 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <darcy_gravity>
- <axis_id>2</axis_id>
- <g>0.0</g>
- </darcy_gravity>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="q"/>
</secondary_variables>
+ <specific_body_force>0.0 0.0 0.0</specific_body_force>
</process>
</processes>
<media>