partial(C i rho) partial t + nabla* [C_{i}*rho*v] = 0

partial( rho v) partial t + nabla*[ rho*vv - BB +P^ * +G]=- rho nabla phi- nabla* Pi

Ot +V.[(E +P^ * )v-B(B* v)]=- rho v* nabla (Phi + phi) + rho*H - rho ^ 2 * Lambda - nabla* Q- nabla*( Pi* v)

partial B partial t - nabla*[ vB - eta*J - eta_{H}(JB) - eta AD (J* B)* B]=0

G = (gg - (gg)/2) / 4 * pi*G g =- nabla Phi

J = nabla* B

nabla^ 2 Phi = 4pi*G*rho

Q =- kappa 0 nabla T- kappa || hat b hat b * nabla T

Pi=- nu 0 nabla v-3 nu || ( hat b hat b - 1/3 |)( hat b hat b - 1 3 |): nabla v