how to treat source term
 

Hello! I am writing a program to solve chemical nonequilibrium problem.When adding the source term,the density becomes negative,and I think it's uncorrect treating of source term.so,could anyone give me some advise,thank you!