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