[T/LiquidFlow] Use specific_body_force instead of darcy_gravity.

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>

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>