Diff_Step_MHD_NS

input: Domain *Omega

output:

Description: This function calculates the diffusion terms in MHD N-S equation.