12 deut.gas_phase *= factor;
13 deut.xMolecules *= factor;
14 deut.xIonDense[0] *= factor;
15 deut.xIonDense[1] *= factor;
22 deut.xIonDense[1] = 0.;
37 realnum neut = total * xNeutral/( xNeutral + xIonized );
38 realnum ionz = total * xIonized/( xNeutral + xIonized );
42 realnum ion_ratio = xIonized/xNeutral;
45 realnum neut = total / ( ion_ratio + 1.f );
46 realnum ionz = total * ion_ratio / ( ion_ratio + 1.f );
48 if( total > 1e-4 *
deut.gas_phase )
50 deut.xIonDense[0] = neut;
51 deut.xIonDense[1] = ionz;
65 deut.fractionation = frac;
73 deut.gas_phase = Hdensity *
deut.fractionation;
void SetGasPhaseDeuterium(const realnum &Hdensity)
void InitDeuteriumIonization()
void SetDeuteriumFractionation(const realnum &frac)
void SetDeuteriumIonization(const double &xNeutral, const double &xIonized)
void ScaleDensitiesDeuterium(const realnum &factor)
void total_molecule_deut(realnum &total)