59 static double gS2[5]={4.,4.,6.,2.,4.};
60 static double exS2[4]={14851.9,31.5,9640.8,48.6};
62 static bool lgFirst=
true;
64 static long int *ipS4Pump=NULL,
117 atom_pop5(gS2,exS2,cs12,cs13,cs14,cs15,cs23,cs24,cs25,cs34,cs35,cs45,
118 8.82e-4,2.60e-4,9.06e-2,0.225,3.35e-7,0.163,0.133,0.078,0.179,
303 else if( nS4Pump > 0 )
305 ipS4Pump = (
long *)
MALLOC((
unsigned)(nS4Pump)*
sizeof(
long) );
315 ipS4Pump[nS4Pump] = i;
326 for( i=0; i<nS4Pump; ++i )
330 fprintf(
ioQQQ,
"DEBUG C %li %.3e %.3e\n",
345 1.168 , 3.366 , 2.924 , 7.233 ,
515 telogn1 =
MAX2(3.5,telog);
516 telogn1 =
MIN2(telogn1,5.0);
522 *cs35 =
MIN2(5.82,a+b*telogn1+c*telogn1*telogn1);
523 *cs35 =
MAX2(3.87,*cs35);
529 *cs25 =
MIN2(3.38,a+b*telogn1+c*telogn1*telogn1);
530 *cs25 =
MAX2(2.24,*cs25);
536 *cs34 =
MIN2(2.46,a+b*telogn1+c*telogn1*telogn1);
537 *cs34 =
MAX2(1.64,*cs34);
543 *cs24 =
MIN2(2.14,a+b*telogn1+c*telogn1*telogn1);
544 *cs24 =
MAX2(1.42,*cs24);
550 *cs45 =
MIN2(3.07,a+b/telogn1+c/(telogn1*telogn1));
551 *cs45 =
MAX2(1.85,*cs45);
557 *cs23 =
MIN2(8.01,a+b/telogn1+c/(telogn1*telogn1));
558 *cs23 =
MAX2(4.79,*cs23);
564 *cs15 = a+b*pow(telogn1,c);
565 *cs15 =
MIN2(2.46,*cs15);
566 *cs15 =
MAX2(1.45,*cs15);
575 *cs13 = a+b*pow(telogn1,c);
576 *cs13 =
MIN2(4.77,*cs13);
577 *cs13 =
MAX2(2.54,*cs13);