32 #if defined (__ICC) && defined(__ia64) && __INTEL_COMPILER < 910
33 #pragma optimization_level 1
78 double cse01=-1.,cse12=-1.,cse02 =-1.,
79 csh01=-1.,cshe01=-1.,csh201=-1.,csh12=-1.,cshe12=-1.,csh212=-1.,csh02=-1.,cshe02=-1.,csh202 =-1.,
80 csh2o01=-1.,csh2o02=-1.,csh2o12=-1.,csh2p01=-1.,csh2p02=-1.,csh2p12=-1.,csp01=-1.,csp02=-1.,
82 static double go2[5]={4.,6.,4.,4.,2.};
83 static double exo2[4]={26808.4,21.0,13637.5,1.5};
86 static double oi_cs_tsave=-1. , oi_te_tsave=-1. , oi_dcdt_tsave=-1.;
88 static bool lgFirst=
true;
89 static long int *ipO4Pump=NULL,
91 double rate_OH_dissoc;
112 cs12 =
MAX2(0.1,cs12);
137 double te_scale =
phycon.
te / 6000.;
138 double rate_H0 = (1.74*te_scale + 0.6)*1e-12*
sexp(0.47*te_scale) / sqrt(te_scale ) *
249 oi3Pcs(&cse01,&cse02,&cse12);
292 csh212 = ortho_frac*csh2o12 + (1.-ortho_frac)*csh2p12;
296 csh201 = ortho_frac*csh2o01 + (1.-ortho_frac)*csh2p01;
300 csh202 = ortho_frac*csh2o02 + (1.-ortho_frac)*csh2p02;
388 (cs01-cs)/cs01*
TauLines[ipT63].cool*
399 oi_dcdt_tsave = (cs01-oi_cs_tsave) / (
phycon.
te-oi_te_tsave);
412 oi_dcdt_tsave =
MAX2( 0. , oi_dcdt_tsave);
422 enum{DEBUG_LOC=
false};
426 fprintf(
ioQQQ,
"DEBUG OI\t%.2f\tte\t%.5e\tcool\t%.5e\tcs\t%.4e\told\t%.4e\tnew\t%.4e\n",
493 atom_pop5(go2,exo2,cs21,cs31,cs41,cs51,cs32,cs42,cs52,cs43,cs53,cs54,
494 a21,a31,a41,a51,a32,a42,a52,a43,a53,a54,p,
dense.
xIonDense[
ipOXYGEN][1]);
513 ( (p[3]*(a41+a42+a43) + p[4]*(a51+a52+a53) ) +
514 ( p[3]*(cs41+cs42+cs43)/go2[3] + p[4]*(cs51+cs52+cs53)/go2[4]) *
525 ( (p[1]*a21 + p[2]*a31 ) +
526 ( p[1]*cs21/go2[1] + p[2]*cs31/go2[2]) *
648 o3cs23 = 0.32412181 + 79.051672/
phycon.
sqrte - 4374.7816/
767 Te_bounded =
MIN2(Te_bounded,450000.);
768 Te_bounded_log = log(Te_bounded);
772 cs = -3.0102462 + 109.22973/Te_bounded_log - 18666.357/Te_bounded;
809 cs =
MAX2( cs , 1e-4 );
827 else if( nO4Pump > 0 )
829 ipO4Pump = (
long *)
MALLOC((
unsigned)(nO4Pump)*
sizeof(
long) );
839 ipO4Pump[nO4Pump] = i;
850 for( i=0; i<nO4Pump; ++i )
854 fprintf(
ioQQQ,
"DEBUG C %li %.3e %.3e\n",
870 0.1367 , 0.1560 , 1.1564 , 0.0457 , pump_rate,
"O 4");