%Slug/Bubble Composition - needs amount inputs
%Author: Tom D Pering - University of Sheffield
%You are free to use and alter with acknowledgement
%inputs aSO2=; aCO2=; aH20=; in kg
function [mmSlug,dSlug]=slug(aSO2,aCO2,aH2O)
load gas_properties.mat;
mmSlug=(((aSO2*1000)/mmSO2)+((aCO2*1000)/mmCO2)+((aH2O*1000)/mmH2O));
total=(aSO2+aCO2+aH2O);
propSO2=aSO2/total;
propCO2=aCO2/total;
propH2O=aH2O/total;
dSlug=(dSO2*propSO2)+(dCO2*propCO2)+(dH2O*propH2O);
clear mmCO mmCO2 mmCl mmF mmH2O mmH2S mmSO2 dCO dCO2 dCl dF dH2O dH2S dSO2 total propSO2 propCO2 propH2O
end