reformulated divergence() and forcing function