9 # pragma warning( disable : 4100 )
39 double (*fun)(
struct COmole_rate_s *rate),
double a,
double b,
double c);
61 if( lgReactInitialized )
64 lgReactInitialized =
true;
1375 int i,j,prod,exists;
1383 p->
p = (
void *) rate;
1399 for(i=0;!i || label[i-1]!=
'\0';i++)
1401 if(label[i] ==
',' || label[i] ==
'=' || label[i] ==
'\0')
1415 fprintf(stderr,
"Mole_co_etc: Too many reactants in %s, only %d allowed\n",label,
MAXREACTANTS);
1423 fprintf(stderr,
"Mole_co_etc: Too many rate species in %s, only %d allowed\n",label,
MAXREACTANTS);
1436 fprintf(stderr,
"Mole_co_etc: Too many products in %s, only %d allowed\n",label,
MAXPRODUCTS);
1447 fprintf(stderr,
"Could not find %s\n",buf);
1448 if(strncmp(buf,
"PHOTON",6) == 0) {
1501 return pow(te/300.,rate->
b)*exp(-rate->
c/te);
1531 for(n=0;n<nreact;n++)
1611 for( ns=0; ns<2; ++ns )
1667 double binding_energy, exponent, vib_freq, number_of_sites ;
1674 binding_energy = rate->
b;
1680 for( nd=0; nd <
gv.
nBin; nd++ )
1682 exponent += exp(-binding_energy/
gv.
bin[nd]->
tedust);
1685 number_of_sites = 1.5e15;
1705 double den_times_area;
1709 den_times_area = 0.0;
1721 for( nd=0; nd <
gv.
nBin; nd++ )