65 static bool lgFirst=
true;
66 static long int *ipC2Pump=NULL,
155 csp01 = 1e-9*pow((4.3671821 - 14.39018/log(temp)),2)*(1./3.)*exp(16.4*
T1CM/temp);
156 csp12 = 1e-9*exp(3.2823932 - 60.99754*(log(temp)/temp))*(1./5.)*exp(37.1*
T1CM/temp);
157 csp02 = 1e-9/(0.033932579+ (1503.4042/pow(temp,1.5)))*(3./5.)*exp(43.5*
T1CM/temp);
182 csh2p01 =
MAX2(6.60e-11,
185 csh2o01 =
MAX2(7.10e-11,
193 csh2p01 =
MIN2(8.10e-11,csh2p01);
196 csh2o01 =
MIN2(8.57e-11,csh2o01);
200 csh201 = ortho_frac*csh2o01 + (1.-ortho_frac)*csh2p01;
208 csh2p12 =
MAX2(1.75e-10,
214 csh2p12 =
MIN2(2.61e-10,csh2p12);
220 csh212 = ortho_frac*csh2o12 + (1.-ortho_frac)*csh2p12;
230 csh2p02 =
MAX2(8.40e-11,
236 csh2p02 =
MIN2(1.04e-10,csh2p02);
241 csh202 = ortho_frac*csh2o02 + (1.-ortho_frac)*csh2p02;
335 else if( nC2Pump > 0 )
337 ipC2Pump = (
long *)
MALLOC((
unsigned)(nC2Pump)*
sizeof(
long) );
347 ipC2Pump[nC2Pump] = i;
358 for( i=0; i<nC2Pump; ++i )
362 fprintf(
ioQQQ,
"DEBUG C %li %.3e %.3e\n",
394 cs_c2_h12 = 1e-10*(4.4716028+ 0.69658785*pow(temp, 0.31692387));
399 cs_c2_h12 *= pow(
phycon.
te/2e3, 0.31692387);
429 0.2349 , 0.8237 , 0.8533 , 1.9818 , pump_rate ,
"C 2");
432 enum {DEBUG_LOC=
false};
433 if( DEBUG_LOC &&
nzone > 80 )
435 fprintf(
ioQQQ,
"DEBUG\t%.2f\t%.3e\t%.3e\t%.2e\t%.2e\t%.2e\t%.2e\n",
503 popratio = cs/(cs + a21);