use constant molecular diffusion coefficient (1.0E-9)
correcting RHS in global residual