output:
Description: This function calculates the diffusion terms in MHD N-S equation, including pot source term.
