function [g_L,g_G]=gravity(alpha_L,rho_G,data)
%   Function computing the gravitational effects on the flow
%
%   Written by Espen Storkaas, Jan 16th 2003
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
g_L=mean([alpha_L(data.I_i) alpha_L(data.I_i+1)],2)*data.rho_L*data.g.*sin(data.theta);
g_G=mean([(1-alpha_L(data.I_i)).*rho_G(data.I_i) (1-alpha_L(data.I_i+1))...
        .*rho_G(data.I_i+1)],2)*data.g.*sin(data.theta);