D_t pi = phi
D_t phi = c**2 Euc^ij D_ij pi