bug fix and golden generation

Author: wilfonba

src/simulation/m_body_forces.fpp

@@ -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