Body Forces Bug Fix (#920)

Author: wilfonba

examples/2D_rayleigh_taylor/README.md

@@ -1,5 +1,4 @@
 # Rayleigh-Taylor Instability (2D)
 
 ## Final Condition and Linear Theory
-<img src='final_condition.png' height='MAX_HEIGHT'/>
-<img src='linear_theory.png' height='MAX_HEIGHT'/>
+<img src='result.png' height='MAX_HEIGHT'/>

examples/2D_rayleigh_taylor/final_condition.png

@@ -98,7 +98,7 @@ contains
     !! so the system can be advanced in time
     !! @param q_cons_vf Conservative variables
     !! @param q_prim_vf Primitive variables
-    subroutine s_compute_body_forces_rhs(q_cons_vf, q_prim_vf, rhs_vf)
+    subroutine s_compute_body_forces_rhs(q_prim_vf, q_cons_vf, rhs_vf)
 
         type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
         type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf