Add external flux term to small-scale shallow water model

same as in the Kinematic Wave solver: add support for source term modeling rain/evaporation/another source of water.