I also implemented the porous media model through a source term rather than using the build in porous media model. Many authors went this route. I only used the built-in model for stability and simplicity of insertion ( as well as told to by my advisor). 


What would you recommend?