% Mass flow equation
function mf = sz_massf(P0,T0,Pv,k,R,A,ws)
if Pv > P0
mf = 0;
return
end
Ff = A * (P0/(R * T0)) * sqrt(1000 * k * R * T0)/ws;
pc = (2/(k+1))^(k/(k-1));
if (Pv/P0) > pc
mf = Ff * ((2/(k-1)) * abs((Pv/P0)^(2/k) - (Pv/P0)^((k+1)/k)))^0.5;
else
mf = Ff * (2/(k+1))^((k+1)/(2*(k-1)));
end
end