G=\frac{g\pi}{6}[D_2^3(\rho_t-\rho_k)-D_1^3\rho_k]