40 char chL[21]={
'S',
'P',
'D',
'F',
'G',
'H',
'I',
'K',
'L',
'M',
'N',
'O',
'Q',
'R',
'T',
'U',
'V',
'W',
'X',
'Y',
'Z'};
49 static int nCalled = 0;
88 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
106 EnergyRyd = HIonPoten/
POW2((
double)
N_(ipHi));
124 for( ipLo=0; ipLo < ipHi; ipLo++ )
135 EnergyWN = -1.0 * EnergyWN;
165 for( ipHi=2; ipHi <
iso.
nLyman[ipISO]; ipHi++ )
204 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
234 Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul = Aul;
257 Transitions[ipISO][nelem][ipHi][ipLo].EnergyWN <= 0.)
259 Transitions[ipISO][nelem][ipHi][ipLo].Emis->gf = 0.;
260 Transitions[ipISO][nelem][ipHi][ipLo].Emis->opacity = 0.;
271 Transitions[ipISO][nelem][ipHi][ipLo].Emis->opacity =
289 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
296 for( ipLo=
ipHe1s1S; ipLo<ipHi; ipLo++ )
321 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
327 StatesElem[ipISO][nelem][0].lifetime = -FLT_MAX;
333 for( ipLo=0; ipLo < ipHi; ipLo++ )
344 for( ipLo=0; ipLo < ipHi; ipLo++ )
346 if(
Transitions[ipISO][nelem][ipHi][ipLo].EnergyWN <= 0. )
353 (1.f/
StatesElem[ipISO][nelem][ipHi].lifetime)/
369 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
378 Transitions[ipISO][nelem][1+ipISO][0].Emis->TauTot = 0.;
382 Transitions[ipISO][nelem][1+ipISO][0].Emis->opacity /= 1e4f;
397 for( nelem = ipISO; nelem <
LIMELM; nelem++ )
418 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
430 iso.
pestrk[ipISO][nelem][ipLo][ipHi] = 0.;
461 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
471 iso.
gamnc[ipISO][nelem][ipHi] = -DBL_MAX;
496 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
508 for( ipLo=0; ipLo < ipHi; ipLo++ )
514 for( ipHi=2; ipHi <
iso.
nLyman[ipISO]; ipHi++ )
533 for( i=0; i <
nUTA; i++ )
602 for( nelem=ipISO; nelem <
LIMELM; ++nelem )
640 for( i1 = 0; i1 < i; ++i1 )
724 for( nelem=ipISO; nelem <
LIMELM; ++nelem )
748 for( nelem=ipISO; nelem <
LIMELM; ++nelem )
765 for( ipLo=0; ipLo < ipHi; ipLo++ )
775 for( ipHi=2; ipHi <
iso.
nLyman[ipISO]; ipHi++ )
790 static bool lgNotSetYet =
true;
805 for( nelem=ipISO; nelem <
LIMELM; ++nelem )
816 for( j = 0; j<i+1; j++ )
859 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
873 for( il = 0L; il < in; ++il )
888 StatesElem[ipISO][nelem][level].S = -LONG_MAX;
889 StatesElem[ipISO][nelem][level].l = -LONG_MAX;
892 for( il = 0; il < in; ++il )
909 for( il = 0L; il < in; ++il )
937 for( il = 0L; il < in; ++il )
939 for( is = 3L; is >= 1L; is -= 2 )
944 if( (il == 1L) && (is == 1L) )
947 if( (in == 1L) && (is == 3L) )
962 else if( (in == 2) && (il == 1) && (is == 3) )
1004 StatesElem[ipISO][nelem][level].S = -LONG_MAX;
1005 StatesElem[ipISO][nelem][level].l = -LONG_MAX;
1008 for( il = 0; il < in; ++il )
1010 for( is = 1; is <= 3; is += 2 )
1028 for( il = 0L; il < in; ++il )
1030 for( is = 3L; is >= 1; is -= 2 )
1033 if( (in == 1L) && (is == 3L) )
1052 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
1059 StatesElem[ipISO][nelem][ipLo].nelem = (int)(nelem+1);
1060 StatesElem[ipISO][nelem][ipLo].IonStg = (int)(nelem+1-ipISO);
1066 else if(
StatesElem[ipISO][nelem][ipLo].l >= 0 )
1083 char chConfiguration[11] =
" ";
1084 long nCharactersWritten = 0;
1089 nCharactersWritten = sprintf( chConfiguration,
"n=%3li",
1092 else if(
StatesElem[ipISO][nelem][ipLo].j > 0 )
1094 nCharactersWritten = sprintf( chConfiguration,
"%3li^%2li%c_%li",
1102 nCharactersWritten = sprintf( chConfiguration,
"%3li^%2li%c",
1108 ASSERT( nCharactersWritten <= 10 );
1109 chConfiguration[10] =
'\0';
1111 strncpy(
StatesElem[ipISO][nelem][ipLo].chConfig, chConfiguration, 10 );
1119 #if defined(__ICC) && defined(__i386)
1120 #pragma optimization_level 1
1129 t->
Hi->
nelem = (int)(nelem+1);
1130 t->
Hi->
IonStg = (int)(nelem+1-ipISO);
1153 Aul = (1.508e10) / pow((
double)nHi,2.975);
1161 Aul = 1.375E10 * pow((
double)nelem, 3.9) / pow((
double)nHi,3.1);
1186 enum {DEBUG_LOC=
false};
1189 fprintf(
ioQQQ,
"%li\t%li\t%.2e\t%.2e\n",
1205 double tau, t0, eps2;
1210 double mu = (m*M)/(M+m);
1212 long Z = nelem + 1 - ipISO;
1216 eps2 = 1. - ( l*l + l + 8./47. - (l+1.)/69./n ) /
POW2( (
double)n );
1218 t0 = 3. *
H_BAR * pow( (
double)n, 5.) /
1220 POW2( (m + M)/(Z*m + z*
M) );
1222 tau = t0 * ( 1. - eps2 ) /
1223 ( 1. + 19./88.*( (1./eps2 - 1.) * log( 1. - eps2 ) + 1. -
1224 0.5 * eps2 - 0.025 * eps2 * eps2 ) );
1231 tau *= 1.1722 * pow( (
double)nelem, 0.1 );
1248 long int i, j, ipLo, ipHi;
1290 for( ipLo=ipLoStart; ipLo<ipHi; ipLo++ )
1292 SumAs +=
Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul;
1295 for( ipLo=ipLoStart; ipLo<ipHi; ipLo++ )
1306 Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul / SumAs;
1310 SumAPerN[
N_(ipHi)] +=
Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul;
1324 (
double)
Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul, 2. );
1335 for( ipLo=0; ipLo<ipHi; ipLo++ )
1337 double SigmaA, SigmaCul = 0.;
1339 for( i=ipLo; i<ipHi; i++ )
1358 SigmaCul = sqrt(SigmaCul);
1368 enum {DEBUG_LOC=
false};
1392 fprintf(
ioQQQ,
"Bm(n,%ld,%ld;%ld)\n",hi_l,hi_s,ipLo);
1393 fprintf(
ioQQQ,
"m\t2\t\t3\t\t4\t\t5\t\t6\n");
1403 for( i = ipLo; i<=ipHi; i++)
1420 fprintf(
ioQQQ,
"%2.4e\t",Bm);
1426 fprintf(
ioQQQ,
"%2.4e\t",Bm);
1433 fprintf(
ioQQQ,
"\n\n");
1448 enum {DEBUG_LOC=
false};
1453 fprintf(
ioQQQ,
"n %ld\t lifetime %.4e\n", i, 1./SumAPerN[i]);
1468 long i, ipISO, nelem;
1474 for( nelem = ipISO; nelem <
LIMELM; nelem++ )
1527 double ConBoltz, LTE_pop=
SMALLFLOAT+FLT_EPSILON, factor1, ConvLTEPOP;
1528 long i, ipISO, nelem;
1535 for( nelem = ipISO; nelem <
LIMELM; nelem++ )
1574 LTE_pop =
SatelliteLines[ipISO][nelem][i].Hi->g * ConBoltz * ConvLTEPOP;
1577 LTE_pop =
max( LTE_pop, 1e-30f );
1602 double lifetime = 1e-10;
1618 long tot_num_levels;
1625 tot_num_levels = (long)( nmaxResolved * 0.5 *( nmaxResolved + 1 ) ) + numCollapsed;
1629 tot_num_levels = nmaxResolved*nmaxResolved + nmaxResolved + 1 + numCollapsed;
1634 return tot_num_levels;
1649 fprintf(
ioQQQ,
"The number of levels for ipISO %li, nelem %li, has been increased since the initial coreload.\n",
1651 fprintf(
ioQQQ,
"This cannot be done.\n" );
1676 double bnl_array[4][3][4][10] = {
1680 {6.13E-01, 2.56E-01, 1.51E-01, 2.74E-01, 3.98E-01, 4.98E-01, 5.71E-01, 6.33E-01, 7.28E-01, 9.59E-01},
1681 {1.31E+00, 5.17E-01, 2.76E-01, 4.47E-01, 5.87E-01, 6.82E-01, 7.44E-01, 8.05E-01, 9.30E-01, 1.27E+00},
1682 {1.94E+00, 7.32E-01, 3.63E-01, 5.48E-01, 6.83E-01, 7.66E-01, 8.19E-01, 8.80E-01, 1.02E+00, 1.43E+00},
1683 {2.53E+00, 9.15E-01, 4.28E-01, 6.16E-01, 7.42E-01, 8.13E-01, 8.60E-01, 9.22E-01, 1.08E+00, 1.56E+00}
1686 {5.63E-01, 2.65E-01, 1.55E-01, 2.76E-01, 3.91E-01, 4.75E-01, 5.24E-01, 5.45E-01, 5.51E-01, 5.53E-01},
1687 {1.21E+00, 5.30E-01, 2.81E-01, 4.48E-01, 5.80E-01, 6.62E-01, 7.05E-01, 7.24E-01, 7.36E-01, 7.46E-01},
1688 {1.81E+00, 7.46E-01, 3.68E-01, 5.49E-01, 6.78E-01, 7.51E-01, 7.88E-01, 8.09E-01, 8.26E-01, 8.43E-01},
1689 {2.38E+00, 9.27E-01, 4.33E-01, 6.17E-01, 7.38E-01, 8.05E-01, 8.40E-01, 8.65E-01, 8.92E-01, 9.22E-01}
1692 {2.97E-01, 2.76E-01, 2.41E-01, 3.04E-01, 3.66E-01, 4.10E-01, 4.35E-01, 4.48E-01, 4.52E-01, 4.53E-01},
1693 {5.63E-01, 5.04E-01, 3.92E-01, 4.67E-01, 5.39E-01, 5.85E-01, 6.10E-01, 6.20E-01, 6.23E-01, 6.23E-01},
1694 {7.93E-01, 6.90E-01, 4.94E-01, 5.65E-01, 6.36E-01, 6.79E-01, 7.00E-01, 7.09E-01, 7.11E-01, 7.11E-01},
1695 {1.04E+00, 8.66E-01, 5.62E-01, 6.31E-01, 7.01E-01, 7.43E-01, 7.63E-01, 7.71E-01, 7.73E-01, 7.73E-01}
1701 {6.70E-02, 2.93E-02, 1.94E-02, 4.20E-02, 7.40E-02, 1.12E-01, 1.51E-01, 1.86E-01, 2.26E-01, 3.84E-01},
1702 {2.39E-01, 1.03E-01, 6.52E-02, 1.31E-01, 2.11E-01, 2.91E-01, 3.61E-01, 4.17E-01, 4.85E-01, 8.00E-01},
1703 {4.26E-01, 1.80E-01, 1.10E-01, 2.09E-01, 3.18E-01, 4.15E-01, 4.93E-01, 5.54E-01, 6.34E-01, 1.04E+00},
1704 {6.11E-01, 2.55E-01, 1.51E-01, 2.74E-01, 3.99E-01, 5.02E-01, 5.80E-01, 6.41E-01, 7.30E-01, 1.21E+00}
1707 {6.79E-02, 3.00E-02, 2.00E-02, 4.30E-02, 7.48E-02, 1.11E-01, 1.44E-01, 1.70E-01, 1.87E-01, 1.96E-01},
1708 {2.40E-01, 1.04E-01, 6.62E-02, 1.32E-01, 2.11E-01, 2.87E-01, 3.51E-01, 3.98E-01, 4.32E-01, 4.58E-01},
1709 {4.26E-01, 1.81E-01, 1.11E-01, 2.10E-01, 3.17E-01, 4.12E-01, 4.89E-01, 5.53E-01, 6.14E-01, 6.84E-01},
1710 {6.12E-01, 2.55E-01, 1.51E-01, 2.73E-01, 3.97E-01, 4.98E-01, 5.77E-01, 6.51E-01, 7.82E-01, 1.18E+00}
1713 {4.98E-02, 3.47E-02, 2.31E-02, 4.54E-02, 7.14E-02, 9.37E-02, 1.08E-01, 1.13E-01, 1.13E-01, 1.11E-01},
1714 {1.75E-01, 1.16E-01, 7.36E-02, 1.36E-01, 2.01E-01, 2.50E-01, 2.76E-01, 2.84E-01, 2.81E-01, 2.77E-01},
1715 {3.38E-01, 1.97E-01, 1.18E-01, 2.13E-01, 3.06E-01, 3.72E-01, 4.06E-01, 4.15E-01, 4.10E-01, 4.04E-01},
1716 {6.01E-01, 2.60E-01, 1.53E-01, 2.76E-01, 3.95E-01, 4.87E-01, 5.45E-01, 5.76E-01, 5.93E-01, 6.05E-01}
1722 {1.77E-01, 3.59E-01, 1.54E-01, 2.75E-01, 3.98E-01, 4.94E-01, 5.51E-01, 5.68E-01, 5.46E-01, 4.97E-01},
1723 {4.09E-01, 7.23E-01, 2.83E-01, 4.48E-01, 5.89E-01, 6.78E-01, 7.22E-01, 7.30E-01, 7.07E-01, 6.65E-01},
1724 {6.40E-01, 1.02E+00, 3.74E-01, 5.49E-01, 6.85E-01, 7.63E-01, 7.98E-01, 8.03E-01, 7.84E-01, 7.53E-01},
1725 {8.70E-01, 1.28E+00, 4.42E-01, 6.17E-01, 7.44E-01, 8.13E-01, 8.42E-01, 8.46E-01, 8.34E-01, 8.13E-01}
1728 {1.78E-01, 3.62E-01, 1.55E-01, 2.73E-01, 3.91E-01, 4.73E-01, 5.10E-01, 5.04E-01, 4.70E-01, 4.32E-01},
1729 {4.08E-01, 7.26E-01, 2.83E-01, 4.45E-01, 5.79E-01, 6.54E-01, 6.78E-01, 6.64E-01, 6.30E-01, 5.98E-01},
1730 {6.37E-01, 1.03E+00, 3.73E-01, 5.46E-01, 6.75E-01, 7.40E-01, 7.57E-01, 7.43E-01, 7.15E-01, 6.92E-01},
1731 {8.65E-01, 1.28E+00, 4.41E-01, 6.14E-01, 7.35E-01, 7.92E-01, 8.05E-01, 7.95E-01, 7.74E-01, 7.59E-01}
1734 {2.07E-01, 3.73E-01, 1.73E-01, 2.85E-01, 4.03E-01, 4.76E-01, 5.06E-01, 5.03E-01, 4.84E-01, 4.63E-01},
1735 {4.32E-01, 7.13E-01, 3.06E-01, 4.54E-01, 5.81E-01, 6.44E-01, 6.59E-01, 6.49E-01, 6.28E-01, 6.11E-01},
1736 {6.40E-01, 9.85E-01, 3.98E-01, 5.53E-01, 6.74E-01, 7.27E-01, 7.36E-01, 7.26E-01, 7.10E-01, 6.98E-01},
1737 {8.38E-01, 1.21E+00, 4.67E-01, 6.20E-01, 7.34E-01, 7.79E-01, 7.87E-01, 7.79E-01, 7.69E-01, 7.63E-01}
1743 {9.31E-02, 3.96E-01, 1.36E-01, 2.74E-01, 3.99E-01, 4.95E-01, 5.52E-01, 5.70E-01, 5.48E-01, 4.96E-01},
1744 {2.25E-01, 8.46E-01, 2.49E-01, 4.46E-01, 5.89E-01, 6.79E-01, 7.23E-01, 7.31E-01, 7.08E-01, 6.64E-01},
1745 {3.59E-01, 1.24E+00, 3.30E-01, 5.47E-01, 6.85E-01, 7.63E-01, 7.98E-01, 8.04E-01, 7.85E-01, 7.53E-01},
1746 {4.93E-01, 1.60E+00, 3.91E-01, 6.15E-01, 7.44E-01, 8.13E-01, 8.42E-01, 8.47E-01, 8.35E-01, 8.12E-01}
1749 {9.32E-02, 3.99E-01, 1.35E-01, 2.72E-01, 3.91E-01, 4.75E-01, 5.14E-01, 5.09E-01, 4.73E-01, 4.31E-01},
1750 {2.25E-01, 8.49E-01, 2.49E-01, 4.44E-01, 5.79E-01, 6.56E-01, 6.81E-01, 6.68E-01, 6.31E-01, 5.96E-01},
1751 {3.58E-01, 1.25E+00, 3.29E-01, 5.44E-01, 6.76E-01, 7.42E-01, 7.60E-01, 7.46E-01, 7.16E-01, 6.91E-01},
1752 {4.92E-01, 1.60E+00, 3.90E-01, 6.12E-01, 7.36E-01, 7.93E-01, 8.07E-01, 7.97E-01, 7.74E-01, 7.58E-01}
1755 {1.13E-01, 4.21E-01, 1.47E-01, 2.83E-01, 4.04E-01, 4.80E-01, 5.13E-01, 5.12E-01, 4.93E-01, 4.71E-01},
1756 {2.52E-01, 8.56E-01, 2.61E-01, 4.50E-01, 5.82E-01, 6.48E-01, 6.66E-01, 6.56E-01, 6.35E-01, 6.16E-01},
1757 {3.85E-01, 1.23E+00, 3.41E-01, 5.49E-01, 6.75E-01, 7.30E-01, 7.41E-01, 7.31E-01, 7.15E-01, 7.02E-01},
1758 {5.14E-01, 1.56E+00, 4.01E-01, 6.15E-01, 7.34E-01, 7.82E-01, 7.90E-01, 7.83E-01, 7.72E-01, 7.65E-01}
1763 double temps[4] = {5000., 10000., 15000., 20000. };
1764 double log_dens[3] = {2., 4., 6.};
1773 ASSERT( (ipTe >=0) && (ipTe < 3) );
1774 ASSERT( (ipDens >=0) && (ipDens < 2) );
1778 for(
long lHi=0; lHi<nHi; lHi++ )
1780 for(
long sHi=1; sHi<4; sHi++ )
1782 if( ipISO ==
ipH_LIKE && sHi != 2 )
1784 else if( ipISO ==
ipHE_LIKE && sHi != 1 && sHi != 3 )
1787 double bnl_at_lo_den, bnl_at_hi_den, bnl;
1788 double bnl_max, bnl_min, temp, dens;
1790 long ipL =
MIN2(9,lHi);
1815 temp =
MIN2( temps[3], temp );
1818 dens =
MIN2( log_dens[2], dens );
1824 if( temp < temps[0] && dens < log_dens[0] )
1825 bnl = bnl_array[ip1][0][0][ipL];
1826 else if( temp < temps[0] && dens >= log_dens[2] )
1827 bnl = bnl_array[ip1][2][0][ipL];
1828 else if( temp >= temps[3] && dens < log_dens[0] )
1829 bnl = bnl_array[ip1][0][3][ipL];
1830 else if( temp >= temps[3] && dens >= log_dens[2] )
1831 bnl = bnl_array[ip1][2][3][ipL];
1834 bnl_at_lo_den = ( temp - temps[ipTe]) / (temps[ipTe+1] - temps[ipTe]) *
1835 (bnl_array[ip1][ipDens][ipTe+1][ipL] - bnl_array[ip1][ipDens][ipTe][ipL]) + bnl_array[ip1][ipDens][ipTe][ipL];
1837 bnl_at_hi_den = ( temp - temps[ipTe]) / (temps[ipTe+1] - temps[ipTe]) *
1838 (bnl_array[ip1][ipDens+1][ipTe+1][ipL] - bnl_array[ip1][ipDens+1][ipTe][ipL]) + bnl_array[ip1][ipDens+1][ipTe][ipL];
1840 bnl = ( dens - log_dens[ipDens]) / (log_dens[ipDens+1] - log_dens[ipDens]) *
1841 (bnl_at_hi_den - bnl_at_lo_den) + bnl_at_lo_den;
1846 bnl_max =
MAX4( bnl_array[ip1][ipDens][ipTe+1][ipL], bnl_array[ip1][ipDens+1][ipTe+1][ipL],
1847 bnl_array[ip1][ipDens][ipTe][ipL], bnl_array[ip1][ipDens+1][ipTe][ipL] );
1848 ASSERT( bnl <= bnl_max );
1850 bnl_min =
MIN4( bnl_array[ip1][ipDens][ipTe+1][ipL], bnl_array[ip1][ipDens+1][ipTe+1][ipL],
1851 bnl_array[ip1][ipDens][ipTe][ipL], bnl_array[ip1][ipDens+1][ipTe][ipL] );
1852 ASSERT( bnl >= bnl_min );
1870 for(
long is = 1; is<=3; ++is)
1874 else if( ipISO ==
ipHE_LIKE && is != 1 && is != 3 )
1877 char chSpin[3][9]= {
"singlets",
"doublets",
"triplets"};
1880 fprintf(
ioQQQ,
" %s %s %s bnl\n",
1886 fprintf(
ioQQQ,
" n\\l=> ");
1889 fprintf(
ioQQQ,
"%2ld ",i);
1891 fprintf(
ioQQQ,
"\n");
1896 if( is==3 && in==1 )
1899 fprintf(
ioQQQ,
" %2ld ",in);
1901 for(
long il = 0; il < in; ++il)
1905 fprintf(
ioQQQ,
"\n");
1918 for(
long ipLo=0; ipLo<ipFirstCollapsed; ipLo++ )
1920 long spin =
StatesElem[ipISO][nelem][ipLo].S;
1941 EffectiveAul /= (2.f*nHi*nHi);
1943 EffectiveAul /= (4.f*nHi*nHi);
1950 Transitions[ipISO][nelem][ipHi][ipLo].Emis->Aul = EffectiveAul;
1967 for(
long ipLo=0; ipLo < ipHi; ipLo++ )
1978 for(
long ipLo=0; ipLo < ipHi; ipLo++ )
1980 if(
Transitions[ipISO][nelem][ipHi][ipLo].EnergyWN <= 0. )
1987 (1.f/
StatesElem[ipISO][nelem][ipHi].lifetime)/
1998 STATIC void Prt_AGN_table(
void )
2007 sprintf( &chLevel.front(ipLo) ,
"%li %li%c",
N_(ipLo),
S_(ipLo),
chL[
MIN2(20,
L_(ipLo))] );
2018 double te[NTEMP]={6000.,8000.,10000.,15000.,20000.,25000. };
2019 double telog[NTEMP] ,
2023 fprintf(
ioQQQ,
"trans");
2024 for(
long i=0; i < NTEMP; ++i )
2026 telog[i] = log10( te[i] );
2027 fprintf(
ioQQQ,
"\t%.3e",te[i]);
2029 for(
long i=0; i < NTEMP; ++i )
2031 fprintf(
ioQQQ,
"\t%.3e",te[i]);
2033 fprintf(
ioQQQ,
"\n");
2037 for(
long ipLo=
ipHe1s1S; ipLo < ipHi; ++ipLo )
2042 if(
N_(ipHi) ==
N_(ipLo) )
2046 fprintf(
ioQQQ,
"%s - %s",
2047 &chLevel.front(ipLo) , &chLevel.front(ipHi) );
2050 for(
long i=0; i < NTEMP; ++i )
2055 fprintf(
ioQQQ,
"\t%.2e", cs );
2059 for(
long i=0; i < NTEMP; ++i )
2068 fprintf(
ioQQQ,
"\t%.2e", ratecoef );
2070 fprintf(
ioQQQ,
"\n");