12 double g1,
double g2,
double g3,
15 double o12,
double o13,
double o23,
18 double a21,
double a31,
double a32,
21 double Tex12,
double Tex23,
55 if( (abund <= 0.) || (ex > 20. && r12<
SMALLFLOAT ) )
63 ASSERT( g1>0. && g2>0. && g3>0. && o12>=0. && o13>=0. && o23>=0. && a21>=0. && a31>=0. && a32>=0. &&
64 Tex12>=0. && Tex23>=0. );
85 alf = a21 + c21 + c23 + gam2;
86 bet = a31 + a32 + c31 + c32;
88 pop3_v = (
atoms.
c13 + *pop2*c23)/bet;
91 fac = abund/(1. + *pop2 + pop3_v);