Fix

Author: efaulhaber

src/schemes/boundary/open_boundary/rhs.jl

@@ -88,6 +88,7 @@ function interact!(dv, v_particle_system, u_particle_system,
                                        v_particle_system, v_neighbor_system,
                                        particle, neighbor, pos_diff, distance,
                                        m_b, rho_a, rho_b, v_a, v_b, grad_kernel)
+        dv[end, particle] += drho_particle[]
 
         # Open boundary pressure evolution matches the corresponding fluid system:
         # - EDAC: Compute pressure evolution like the fluid system