For all boundary conditions, fluxes into the model are negative. For a mass source as a second kind boundary condition, you have to be careful with the form of the transport equation used: The default 'convective' form will only be OK if you do not have a flow boundary condition on the same nodes!