input: Element_List **Us, Element_List **Uf
output:
Description: This function calculates the divergence of flux in scalar equations.
