#include "rt_nonfinite.h"

#include "GetAirDensity.h"

#include "GetAirThermalConduct.h"

#include "GetAirViscosity.h"

#include "GetCompCombGas.h"

#include "GetCorrAirComp.h"

#include "GetDensityHCVapor.h"

#include "GetDiffusivityHCVaporToAir.h"

#include "GetEquilGrill.h"

#include "GetEquilOlikara.h"

#include "GetFuelPropertyN_Dodecane.h"

#include "GetHTCoeffHTX.h"

#include "GetIdealNozzleFlow.h"

#include "GetIdealNozzleFlowPTF.h"

#include "GetMEMbZach.h"

#include "GetPTF.h"

#include "GetPTFV1.h"

#include "GetPTx.h"

#include "GetTFromPhF.h"

#include "GetT_atm_p.h"

#include "GetThdynCombGasZach.h"

#include "GetThdynCombGasZachV1.h"

#include "GetThermalPropertyHCVaporReal.h"

#include "GetThermoDynProp.h"

#include "GetThermoDynPropPartial.h"

#include "GetTotalStaticPT.h"

#include "GetViscosityHCVapor.h"

#include "GetCompleteComb.h"

void GetCompleteComb(double F, const double Fc[4], const double Ac[11], double

x[12])

{

double N_L;

double Fa;

double Fb;

double n_H_2;

double n_H_2O;

double n_N_2;

double n_CO;

memset(&x[0], 0, 12U * sizeof(double));

N_L = ((Fc[1] - 2.0 * Fc[2]) + 4.0 * Fc[0]) / ((((((2.0 * Ac[7] - 2.0 * Ac[5])

- Ac[6]) - 2.0 * Ac[2]) + Ac[9]) + 2.0 * Ac[10]) + 4.0);

Fa = -((((((2.0 * Ac[3] - 2.0 * Ac[5]) - Ac[6]) + 2.0 * Ac[7]) + Ac[9]) + 2.0 *

Ac[10]) + 4.0) / ((2.0 * Fc[2] - Fc[1]) - 2.0 * Fc[0]) * N_L;

Fb = -(((((Ac[3] + Ac[4]) + Ac[7]) + Ac[9]) + Ac[10]) + 2.0) / (Fc[2] - Fc[0])

* N_L;

if (F <= 1.0) {

n_H_2 = F * Fc[0] + (Ac[2] + Ac[3]) * N_L;

n_H_2O = F * Fc[1] / 2.0 + (((Ac[4] + Ac[5]) + Ac[6] / 2.0) + Ac[9] / 2.0) *

N_L;

n_N_2 = F * Fc[3] / 2.0 + ((Ac[0] + Ac[8] / 2.0) + Ac[10] / 2.0) * N_L;

Fb = Ac[1] * N_L;

Fa = (1.0 - F) * ((Fc[1] - 2.0 * Fc[2]) + 4.0 * Fc[0]) / 4.0;

N_L = (((n_H_2 + n_H_2O) + n_N_2) + Fb) + Fa;

x[7] = Fa / N_L;

x[9] = n_H_2 / N_L;

x[8] = n_H_2O / N_L;

x[10] = n_N_2 / N_L;

x[11] = Fb / N_L;

} else if (F <= Fa) {

n_H_2 = (Fc[0] + N_L * (Ac[2] + Ac[3])) * ((Fa - F) / (Fa - 1.0));

n_CO = (F - 1.0) * ((Fc[1] - 2.0 * Fc[2]) + 4.0 * Fc[0]) / 2.0;

n_H_2O = F * Fc[1] / 2.0 + (((Ac[4] + Ac[5]) + Ac[6] / 2.0) + Ac[9] / 2.0) *

N_L;

n_N_2 = F * Fc[3] / 2.0 + ((Ac[0] + Ac[8] / 2.0) + Ac[10] / 2.0) * N_L;

Fb = Ac[1] * N_L;

N_L = (((n_H_2 + n_H_2O) + n_CO) + n_N_2) + Fb;

x[9] = n_H_2 / N_L;

x[8] = n_H_2O / N_L;

x[5] = n_CO / N_L;

x[10] = n_N_2 / N_L;

x[11] = Fb / N_L;

} else {

if (F <= Fb) {

n_CO = F * Fc[0] + (Ac[2] + Ac[3]) * N_L;

n_H_2 = (Fb * Fc[1] / 2.0 + (((Ac[4] + Ac[5]) + Ac[6] / 2.0) + Ac[9] / 2.0)

* N_L) / (Fb - Fa) * (F - Fa);

n_H_2O = (Fa * Fc[1] / 2.0 + (((Ac[4] + Ac[5]) + Ac[6] / 2.0) + Ac[9] /

2.0) * N_L) * ((Fb - F) / (Fb - Fa));

n_N_2 = F * Fc[3] / 2.0 + ((Ac[0] + Ac[8] / 2.0) + Ac[10] / 2.0) * N_L;

Fb = Ac[1] * N_L;

N_L = (((n_CO + n_H_2) + n_H_2O) + n_N_2) + Fb;

x[3] = n_H_2 / N_L;

x[8] = n_H_2O / N_L;

x[5] = n_CO / N_L;

x[10] = n_N_2 / N_L;

x[11] = Fb / N_L;

}

}

}