39 double tnurbn[
NRUBIN] = {1.05E-08,1.05E-07,1.05E-06,1.04E-05,1.00E-04,1.00E-03,1.00E-02,3.01E-02,1.00E-01,
40 1.50E-01,2.50E-01,4.01E-01,6.01E-01,9.8E-01,9.96E-01,1.00E+00,1.02445,1.07266,1.12563,1.18411,1.23881,
41 1.29328,1.35881,1.42463,1.48981,1.55326,1.6166,1.68845,1.76698,1.8019,1.808,1.84567,1.9317,2.04891,2.14533,
42 2.19702,2.27941,2.37438,2.43137,2.49798,2.56113,2.59762,2.66597,2.80543,2.95069,3.02911,3.11182,3.22178,
43 3.3155,3.42768,3.50678,3.56157,3.61811,3.75211,3.89643,4.},
44 fnurbn[
NRUBIN] = {1.56E-20,1.55E-17,1.54E-14,1.53E-11,1.35E-08,1.34E-05,1.35E-02,3.62E-01,1.27E+01,
45 3.90E+01,1.48E+02,4.52E+02,1.02E+03,2.27E+03,8.69E+02,8.04E+02,6.58E+02,6.13E+02,6.49E+02,6.06E+02,
46 6.30E+02,5.53E+02,5.55E+02,5.24E+02,4.86E+02,4.49E+02,4.42E+02,3.82E+02,3.91E+02,2.91E+02,2.61E+02,
47 1.32E+02,1.20E+02,1.16E+02,9.56E+01,9.94E+01,9.10E+01,4.85E+01,7.53E+01,9.53E+01,8.52E+01,5.76E+01,
48 6.72E+01,5.20E+01,8.09E+00,3.75E+00,5.57E+00,3.80E+00,2.73E+00,2.22E+00,3.16E-01,1.26E-01,1.39E-01,
49 6.15E-02,3.22E-02,7.98E-03};
71 static double tnuHM96[
NHM96]={-8,-1.722735683,-0.351545683,-0.222905683,-0.133385683,
73 -0.127655683,-0.004575683,0.297544317,0.476753,0.476756,0.588704317,
74 0.661374317,1.500814317,2.245164317};
77 static double fnuHM96[
NHM96]={-32.53342863,-19.9789,-20.4204,-20.4443,-20.5756,-20.7546,
78 -21.2796,-21.6256,-21.8404,-21.4823,-22.2102,-22.9263,-23.32,-24.2865};
95 long int *nhm ,
double redshift )
99 double *table_redshifts = NULL,
104 long int nRedshift , i , n , nz , ipLo , ipHi;
116 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , ioFILE ) != NULL )
120 if( chLine[0] !=
'#')
129 mag_read = (long)
FFmtRead(chLine,&i,(
int)
sizeof(chLine),&lgEOL);
130 if( mag_read != 50808 )
133 " Magic number in Haardt & Madau file is not correct, I expected 50808 and found %li\n",
148 a =
FFmtRead(chLine,&i,(
int)
sizeof(chLine),&lgEOL);
157 table_redshifts = (
double*)
MALLOC( (
size_t)nRedshift*
sizeof(double) );
160 for( n=0; n<nRedshift; ++n )
162 table_redshifts[n] =
FFmtRead(chLine,&i,(
int)
sizeof(chLine),&lgEOL);
172 table_wl = (
double*)
MALLOC( (
size_t)*nhm*
sizeof(double) );
174 table_fn = (
double**)
MALLOC( (
size_t)nRedshift*
sizeof(
double*) );
175 for(n=0; n<nRedshift; ++n )
177 table_fn[n] = (
double*)
MALLOC( (
size_t)(*nhm)*
sizeof(double) );
181 if( fseek( ioFILE , 0 , SEEK_SET ) != 0 )
183 fprintf(
ioQQQ,
" read_hm05 could not rewind HM05 date file.\n");
188 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , ioFILE ) != NULL )
192 if( chLine[0] !=
'#')
200 table_wl[*nhm] =
FFmtRead(chLine,&i,(
int)
sizeof(chLine),&lgEOL);
203 for( nz=0; nz<nRedshift; ++nz )
206 table_fn[nz][*nhm] =
FFmtRead(chLine,&i,(
int)
sizeof(chLine),&lgEOL);
217 enum {DEBUG_LOC=
false};
220 fprintf(
ioQQQ,
"wavelength/z");
221 for(nz=0; nz<nRedshift; ++nz )
222 fprintf(
ioQQQ,
"\t%.3f", table_redshifts[nz] );
224 for( i=0; i<*nhm; ++i )
226 fprintf(
ioQQQ,
"%.3e", table_wl[i] );
227 for( nz=0; nz<nRedshift; ++nz )
228 fprintf(
ioQQQ,
"\t%.3e", table_fn[nz][i] );
235 assert( table_redshifts!=NULL );
238 if( redshift < table_redshifts[0] ||
239 redshift > table_redshifts[nRedshift-1] )
241 fprintf(
ioQQQ,
" The redshift requested on table HM05 is %g but is not within the range of the table, which goes from %g to %g.\n",
242 redshift, table_redshifts[0] , table_redshifts[nRedshift-1] );
243 fprintf(
ioQQQ,
" Sorry.\n");
250 for( nz=0; nz<nRedshift-1; ++nz )
252 if( redshift >= table_redshifts[nz] &&
253 redshift <= table_redshifts[nz+1] )
260 ASSERT( ipLo>=0 && ipHi >=0 );
264 for( i=0; i<*nhm-1; ++i )
266 if( table_wl[i]>=table_wl[i+1] )
268 fprintf(
ioQQQ,
" The wavelengths in the table HM05 data table are not in increasing order. This is required.\n");
269 fprintf(
ioQQQ,
" Sorry.\n");
275 *fnuHM = (
double *)
MALLOC( (
size_t)(*nhm)*
sizeof(double ) );
276 *tnuHM = (
double *)
MALLOC( (
size_t)(*nhm)*
sizeof(double ) );
281 frac_hi = (redshift-table_redshifts[ipLo]) / (table_redshifts[ipHi]-table_redshifts[ipLo]);
282 for( i=0; i<*nhm; ++i )
287 (*tnuHM)[*nhm-1-i] =
RYDLAM / table_wl[i];
289 (*fnuHM)[*nhm-1-i] = table_fn[ipLo][i]*(1.-frac_hi) + table_fn[ipHi][i]*frac_hi;
293 for( n=0; n<nRedshift; ++n )
297 free( table_redshifts );
329 bool lgNoContinuum =
false,
343 fprintf(
ioQQQ,
" %ld is too many spectra entered. Increase LIMSPC\n Sorry.\n",
352 if(
GetQuote( chFile , chCard ,
false ) )
353 lgQuoteFound =
false;
367 if(
nMatch(
" AGN",chCard) )
371 for( i=0; i <
NAGN; i++ )
379 if(
nMatch(
"BREA",chCard) )
387 if(
nMatch(
" NO ",chCard) )
393 fprintf(
ioQQQ,
" There must be a number for the break.\n Sorry.\n" );
400 fprintf(
ioQQQ,
" The break must be greater than 0.2 Ryd.\n Sorry.\n" );
404 if(
nMatch(
"MICR",chCard) )
407 ConBreak = 0.0912/ConBreak;
413 ConBreak = pow(10.,ConBreak);
416 else if( ConBreak == 0. )
418 fprintf(
ioQQQ,
" An energy of 0 is not allowed.\n Sorry.\n" );
424 fprintf(
ioQQQ,
" The energy of the break cannot be greater than%10.2e Ryd.\n Sorry.\n",
431 fprintf(
ioQQQ,
" The energy of the break cannot be less than%10.2e Ryd.\n Sorry.\n",
448 else if(
nMatch(
"AKN1",chCard) )
452 for( i=0; i <
NKN120; i++ )
460 else if(
nMatch(
"CRAB",chCard) )
464 for( i=0; i <
NCRAB; i++ )
472 else if(
nMatch(
"COOL",chCard) )
476 for( i=0; i <
NCFL; i++ )
484 else if( (i=
nMatch(
"HM96",chCard))>0 )
495 for( j=0; j <
NHM96; j++ )
509 scale = log10(scale);
514 fprintf(
ioQQQ,
" %ld is too many continua entered. Increase LIMSPC\n Sorry.\n",
524 fprintf(
ioQQQ,
" This command has come between a previous ordered pair of continuum shape and luminosity commands.\n Reorder the commands to complete each continuum specification before starting another.\n" );
525 fprintf(
ioQQQ,
" Sorry.\n" );
561 else if(
nMatch(
"HM05",chCard) )
563 double *tnuHM , *fnuHM;
572 if(
nMatch(
"QUAS",chCard) )
575 strcpy( chFile ,
"haardt_madau_quasar.dat" );
580 strcpy( chFile ,
"haardt_madau_galaxy.dat" );
593 fprintf(
ioQQQ,
" The redshift MUST be specified on the table HM05 command. I did not find one.\n Sorry.\n");
599 read_hm05( chFile , &tnuHM , &fnuHM , &nhm , redshift );
602 for( j=0; j < nhm-1; j++ )
606 if( tnuHM[j]<=1. && 1. <= tnuHM[j+1] )
619 scale = log10(scale);
631 fprintf(
ioQQQ,
" %ld is too many continua entered. Increase LIMSPC\n Sorry.\n",
641 fprintf(
ioQQQ,
" This command has come between a previous ordered pair of continuum shape and luminosity commands.\n Reorder the commands to complete each continuum specification before starting another.\n" );
642 fprintf(
ioQQQ,
" Sorry.\n" );
654 scale =
SDIV( fnuHM[ipNORM] );
665 for( j=0; j < nhm; j++ )
678 for( j=0; j < nhm-1; j++ )
701 else if(
nMatch(
" ISM",chCard) )
709 for( i=6; i <
NISM; i++ )
724 if( scale > 0. && !
nMatch(
" LOG",chCard) )
725 scale = log10(scale);
730 fprintf(
ioQQQ,
" %4ld is too many continua entered. Increase LIMSPC\n Sorry.\n",
740 fprintf(
ioQQQ,
" This command has come between a previous ordered pair of continuum shape and luminosity commands.\n Reorder the commands to complete each continuum specification before starting another.\n" );
741 fprintf(
ioQQQ,
" Sorry.\n" );
780 else if(
nMatch(
"DRAI",chCard) )
799 if( scale > 0. && !
nMatch(
" LOG",chCard) )
800 scale = log10(scale);
805 fprintf(
ioQQQ,
" %4ld is too many continua entered. Increase LIMSPC\n Sorry.\n",
815 fprintf(
ioQQQ,
" This command has come between a previous ordered pair of continuum shape and luminosity commands.\n Reorder the commands to complete each continuum specification before starting another.\n" );
816 fprintf(
ioQQQ,
" Sorry.\n" );
856 else if(
nMatch(
"LINE",chCard) )
860 static bool lgLines_Malloc =
false;
863 lgNoContinuum =
true;
868 fprintf(
ioQQQ,
" sorry, only one table line per input stream\n");
879 lgLines_Malloc =
true;
888 fprintf(
ioQQQ,
"\n DISASTER PROBLEM ParseTable could not find "
890 fprintf(
ioQQQ,
" Please check the spelling of the file name and that it "
891 "is in either the local or data directory.\n\n");
896 else if(
nMatch(
"POWE",chCard) )
924 else if( brakmm == 0. )
933 else if( brakmm < 0. )
937 brakmm = pow(10.,brakmm);
941 if(
nMatch(
"MICR",chCard) )
942 brakmm =
RYDLAM / (1e4*brakmm);
949 " Check the order of parameters on this table power law command - the low-energy break of %f Ryd seems high to me.\n",
977 else if( brakxr == 0. )
985 brakxr = pow(10.,brakxr);
993 if( brakmm >= brakxr )
995 fprintf(
ioQQQ,
" HELP!! The lower energy for the power law, %f, is greater than the upper energy, %f. This is not possible.\n Sorry.\n",
1016 else if(
nMatch(
"READ",chCard) )
1021 fprintf(
ioQQQ,
" Oops, there are more than one TABLE READ command in this input stream. I can only deal with a single TABLE READ.\n Sorry.\n");
1031 fprintf(
ioQQQ,
" Name of file must appear on TABLE READ.\n");
1053 else if(
nMatch(
"TLUSTY",chCard) && !
nMatch(
"STAR",chCard) )
1056 fprintf(
ioQQQ,
" The TABLE TLUSTY command is no longer supported.\n" );
1057 fprintf(
ioQQQ,
" Please use TABLE STAR TLUSTY instead. See Hazy for details.\n" );
1061 else if(
nMatch(
"RUBI",chCard) )
1064 for( i=0; i <
NRUBIN; i++ )
1074 else if(
nMatch(
"STAR",chCard) )
1076 char *ptr, chMetalicity[5] =
"", chODFNew[8], chVaryFlag[6] =
"";
1077 bool lgPG1159 =
false, lgHCa =
false, lgHNi =
false, lgSolar, lgHalo, lgList;
1079 double Tlow = -1., Thigh = -1.;
1083 if( (ptr = strstr(chCard,
"1-DI")) != NULL )
1085 else if( (ptr = strstr(chCard,
"2-DI")) != NULL )
1087 else if( (ptr = strstr(chCard,
"3-DI")) != NULL )
1089 else if( (ptr = strstr(chCard,
"4-DI")) != NULL )
1095 strncpy(chVaryFlag,ptr,5);
1096 chVaryFlag[5] =
'\0';
1101 if(
nMatch(
"TIME" , chCard ) )
1108 if( (ptr = strstr(chCard,
"Z+1.0")) != NULL )
1109 strncpy( chMetalicity,
"p10",
sizeof(chMetalicity) );
1110 else if( (ptr = strstr(chCard,
"Z+0.75")) != NULL )
1111 strncpy( chMetalicity,
"p075",
sizeof(chMetalicity) );
1112 else if( (ptr = strstr(chCard,
"Z+0.5")) != NULL )
1113 strncpy( chMetalicity,
"p05",
sizeof(chMetalicity) );
1114 else if( (ptr = strstr(chCard,
"Z+0.4")) != NULL )
1115 strncpy( chMetalicity,
"p04",
sizeof(chMetalicity) );
1116 else if( (ptr = strstr(chCard,
"Z+0.3")) != NULL )
1117 strncpy( chMetalicity,
"p03",
sizeof(chMetalicity) );
1118 else if( (ptr = strstr(chCard,
"Z+0.25")) != NULL )
1119 strncpy( chMetalicity,
"p025",
sizeof(chMetalicity) );
1120 else if( (ptr = strstr(chCard,
"Z+0.2")) != NULL )
1121 strncpy( chMetalicity,
"p02",
sizeof(chMetalicity) );
1122 else if( (ptr = strstr(chCard,
"Z+0.1")) != NULL )
1123 strncpy( chMetalicity,
"p01",
sizeof(chMetalicity) );
1124 else if( (ptr = strstr(chCard,
"Z+0.0")) != NULL )
1125 strncpy( chMetalicity,
"p00",
sizeof(chMetalicity) );
1126 else if( (ptr = strstr(chCard,
"Z-0.1")) != NULL )
1127 strncpy( chMetalicity,
"m01",
sizeof(chMetalicity) );
1128 else if( (ptr = strstr(chCard,
"Z-0.2")) != NULL )
1129 strncpy( chMetalicity,
"m02",
sizeof(chMetalicity) );
1130 else if( (ptr = strstr(chCard,
"Z-0.25")) != NULL )
1131 strncpy( chMetalicity,
"m025",
sizeof(chMetalicity) );
1132 else if( (ptr = strstr(chCard,
"Z-0.3")) != NULL )
1133 strncpy( chMetalicity,
"m03",
sizeof(chMetalicity) );
1134 else if( (ptr = strstr(chCard,
"Z-0.4")) != NULL )
1135 strncpy( chMetalicity,
"m04",
sizeof(chMetalicity) );
1136 else if( (ptr = strstr(chCard,
"Z-0.5")) != NULL )
1137 strncpy( chMetalicity,
"m05",
sizeof(chMetalicity) );
1138 else if( (ptr = strstr(chCard,
"Z-0.7")) != NULL )
1139 strncpy( chMetalicity,
"m07",
sizeof(chMetalicity) );
1140 else if( (ptr = strstr(chCard,
"Z-0.75")) != NULL )
1141 strncpy( chMetalicity,
"m075",
sizeof(chMetalicity) );
1142 else if( (ptr = strstr(chCard,
"Z-1.0")) != NULL )
1143 strncpy( chMetalicity,
"m10",
sizeof(chMetalicity) );
1144 else if( (ptr = strstr(chCard,
"Z-1.3")) != NULL )
1145 strncpy( chMetalicity,
"m13",
sizeof(chMetalicity) );
1146 else if( (ptr = strstr(chCard,
"Z-1.5")) != NULL )
1147 strncpy( chMetalicity,
"m15",
sizeof(chMetalicity) );
1148 else if( (ptr = strstr(chCard,
"Z-1.7")) != NULL )
1149 strncpy( chMetalicity,
"m17",
sizeof(chMetalicity) );
1150 else if( (ptr = strstr(chCard,
"Z-2.0")) != NULL )
1151 strncpy( chMetalicity,
"m20",
sizeof(chMetalicity) );
1152 else if( (ptr = strstr(chCard,
"Z-2.5")) != NULL )
1153 strncpy( chMetalicity,
"m25",
sizeof(chMetalicity) );
1154 else if( (ptr = strstr(chCard,
"Z-3.0")) != NULL )
1155 strncpy( chMetalicity,
"m30",
sizeof(chMetalicity) );
1156 else if( (ptr = strstr(chCard,
"Z-3.5")) != NULL )
1157 strncpy( chMetalicity,
"m35",
sizeof(chMetalicity) );
1158 else if( (ptr = strstr(chCard,
"Z-4.0")) != NULL )
1159 strncpy( chMetalicity,
"m40",
sizeof(chMetalicity) );
1160 else if( (ptr = strstr(chCard,
"Z-4.5")) != NULL )
1161 strncpy( chMetalicity,
"m45",
sizeof(chMetalicity) );
1162 else if( (ptr = strstr(chCard,
"Z-5.0")) != NULL )
1163 strncpy( chMetalicity,
"m50",
sizeof(chMetalicity) );
1164 else if( (ptr = strstr(chCard,
"Z-INF")) != NULL )
1165 strncpy( chMetalicity,
"m99",
sizeof(chMetalicity) );
1167 strncpy( chMetalicity,
"p00",
sizeof(chMetalicity) );
1172 strncpy(chVaryFlag,ptr,5);
1173 chVaryFlag[5] =
'\0';
1180 if(
nMatch(
"HALO",chCard) )
1184 lgSolar = ( !lgHalo && strcmp( chMetalicity,
"p00" ) == 0 );
1187 if( (ptr = strstr(chCard,
"PG1159")) != NULL )
1194 if(
nMatch(
"LIST",chCard) )
1199 if(
nMatch(
"AVAI",chCard) )
1207 for( nval=0; nval <
MDIM; nval++ )
1214 if( nval == 0 && !lgList )
1216 fprintf(
ioQQQ,
" The stellar temperature MUST be entered.\n" );
1226 if( ( val[0] <= 10. && !
nMatch(
"LINE",chCard) ) ||
nMatch(
"LOG ",chCard) )
1227 val[0] = pow(10.,val[0]);
1234 else if(
nMatch(
"ATLA",chCard) )
1240 if(
nMatch(
"ODFN",chCard) )
1241 strncpy( chODFNew,
"_odfnew",
sizeof(chODFNew) );
1243 strncpy( chODFNew,
"",
sizeof(chODFNew) );
1246 nstar =
AtlasInterpolate(val,&nval,&ndim,chMetalicity,chODFNew,lgList,&Tlow,&Thigh);
1249 else if(
nMatch(
"COST",chCard) )
1255 if(
nMatch(
"INDE",chCard) )
1266 if( val[1] < 1. || val[1] > 7. )
1268 fprintf(
ioQQQ,
" The second number must be the id sequence number, 1 to 7.\n" );
1269 fprintf(
ioQQQ,
" reset to 1.\n" );
1273 else if(
nMatch(
"ZAMS",chCard) )
1278 fprintf(
ioQQQ,
" A second number, the age of the star, must be present.\n" );
1282 else if(
nMatch(
" AGE",chCard) )
1287 fprintf(
ioQQQ,
" A second number, the ZAMS mass of the star, must be present.\n" );
1309 if( ! ( lgSolar || lgHalo ) )
1311 fprintf(
ioQQQ,
" Please choose SOLAR or HALO abundances.\n" );
1318 else if(
nMatch(
"KURU",chCard) )
1323 else if(
nMatch(
"MIHA",chCard) )
1328 else if(
nMatch(
"RAUC",chCard) )
1330 if( ! ( lgSolar || lgHalo ) )
1332 fprintf(
ioQQQ,
" Please choose SOLAR or HALO abundances.\n" );
1344 else if(
nMatch(
"HYDR",chCard) )
1348 else if(
nMatch(
"HELI",chCard) )
1352 else if(
nMatch(
"H+HE",chCard) )
1367 else if(
nMatch(
"TLUS",chCard) )
1369 if(
nMatch(
"BSTA",chCard) )
1375 else if(
nMatch(
"OSTA",chCard) )
1383 fprintf(
ioQQQ,
" There must be a third key on TABLE STAR TLUSTY;" );
1384 fprintf(
ioQQQ,
" the options are BSTAR, OSTAR.\n" );
1389 else if(
nMatch(
"WERN",chCard) )
1396 else if(
nMatch(
"WMBA",chCard) )
1405 fprintf(
ioQQQ,
" There must be a second key on TABLE STAR;" );
1406 fprintf(
ioQQQ,
" the options are ATLAS, KURUCZ, MIHALAS, RAUCH, WERNER, and WMBASIC.\n" );
1426 for( i=1; i < nval; i++ )
1430 if(
nMatch(
"ATLA",chCard) )
1436 if(
nMatch(
"ODFN",chCard) )
1446 else if(
nMatch(
"COST",chCard) )
1476 else if(
nMatch(
"KURU",chCard) )
1481 "TABLE STAR KURUCZ %f LOG" );
1484 else if(
nMatch(
"MIHA",chCard) )
1489 "TABLE STAR MIHALAS %f LOG" );
1492 else if(
nMatch(
"RAUC",chCard) )
1499 if(
nMatch(
"HYDR",chCard) )
1501 else if(
nMatch(
"HELI",chCard) )
1503 else if(
nMatch(
"H+HE",chCard) )
1510 if( ( lgHCa || lgHNi ) && ndim == 2 )
1522 if(
nMatch(
"H+HE",chCard) )
1529 if(
nMatch(
"TLUS",chCard) )
1534 if(
nMatch(
"BSTA",chCard) )
1536 else if(
nMatch(
"OSTA",chCard) )
1547 else if(
nMatch(
"WERN",chCard) )
1552 "TABLE STAR WERNER %f LOG , LOG(G)= %f" );
1555 else if(
nMatch(
"WMBA",chCard) )
1560 "TABLE STAR WMBASIC %f LOG , LOG(G)= %f , LOG(Z)= %f" );
1570 for( i=1; i < nval; i++ )
1584 fprintf(
ioQQQ,
" There MUST be a keyword on this line. The keys are:"
1585 "AGN, AKN120, CRAB, COOL, DRAINE, HM96, HM05, _ISM, LINE, POWERlaw, "
1586 "READ, RUBIN, and STAR.\n Sorry.\n" );
1610 for( i=0; i < (ncont + 1); i++ )
1620 for( i=0; i < (ncont + 1); i++ )
1628 for( i=0; i < (ncont + 1); i++ )
1633 for( i=0; i < ncont; i++ )
1653 fprintf(
ioQQQ,
" Table for this continuum; TNU, TFAC, TSLOP\n" );
1654 for( i=0; i < (ncont + 1); i++ )
1669 for( i=0; i < (ncont + 1) && !lgHit; i++ )
1678 if( ncont > 0 && lgHit )
1681 for( i=0; i < (ncont + 1); i++ )
1702 const double RESETFACTOR = 0.98;
1711 power = (fluxlog[1] - fluxlog[0] ) / log10( tnu[1]/tnu[0] );
1714 fluxlog[0] = fluxlog[1] + power * log10( tnu[0] /tnu[1] );
1720 power = log10( fluxlog[1]/fluxlog[0]) / log10( tnu[1]/tnu[0] );
1723 fluxlog[0] = log10(fluxlog[1]) + power * log10( tnu[0] /tnu[1] );
1725 fluxlog[0] = pow(10. , fluxlog[0]);
1765 ++i;
tsldrn[i] = -17.7575;
1766 ++i;
tsldrn[i] = -17.7268;
1767 ++i;
tsldrn[i] = -17.7036;
1768 ++i;
tsldrn[i] = -17.6953;
1769 ++i;
tsldrn[i] = -17.7182;
1770 ++i;
tsldrn[i] = -17.7524;
1771 ++i;
tsldrn[i] = -17.8154;
1772 ++i;
tsldrn[i] = -17.9176;
1773 ++i;
tsldrn[i] = -18.1675;
1774 ++i;
tsldrn[i] = -18.1690;
1775 ++i;
tsldrn[i] = -18.2477;
1776 ++i;
tsldrn[i] = -18.4075;
1777 ++i;
tsldrn[i] = -18.5624;
1778 ++i;
tsldrn[i] = -18.7722;
1834 cfle[0] = 0.0000100;
1835 cflf[0] = -0.8046910;
1836 cfle[1] = 0.7354023;
1837 cflf[1] = -0.8046910;
1838 cfle[2] = 1.4708046;
1839 cflf[2] = -0.7436830;
1840 cfle[3] = 2.2062068;
1841 cflf[3] = -0.6818757;
1842 cfle[4] = 2.9416091;
1843 cflf[4] = -0.7168990;
1844 cfle[5] = 3.6770115;
1845 cflf[5] = -0.8068384;
1846 cfle[6] = 4.4124136;
1847 cflf[6] = -0.6722584;
1848 cfle[7] = 5.1478162;
1849 cflf[7] = -0.7626385;
1850 cfle[8] = 5.8832183;
1851 cflf[8] = -1.0396487;
1852 cfle[9] = 6.6186204;
1853 cflf[9] = -0.7972314;
1854 cfle[10] = 7.3540230;
1855 cflf[10] = -0.9883416;
1856 cfle[11] = 14.7080460;
1857 cflf[11] = -1.1675659;
1858 cfle[12] = 22.0620689;
1859 cflf[12] = -1.1985949;
1860 cfle[13] = 29.4160919;
1861 cflf[13] = -1.2263466;
1862 cfle[14] = 36.7701149;
1863 cflf[14] = -1.2918345;
1864 cfle[15] = 44.1241379;
1865 cflf[15] = -1.3510833;
1866 cfle[16] = 51.4781609;
1867 cflf[16] = -1.2715496;
1868 cfle[17] = 58.8321838;
1869 cflf[17] = -1.1098027;
1870 cfle[18] = 66.1862030;
1871 cflf[18] = -1.4315782;
1872 cfle[19] = 73.5402298;
1873 cflf[19] = -1.1327956;
1874 cfle[20] = 147.080459;
1875 cflf[20] = -1.6869649;
1876 cfle[21] = 220.620681;
1877 cflf[21] = -2.0239367;
1878 cfle[22] = 294.160919;
1879 cflf[22] = -2.2130392;
1880 cfle[23] = 367.701141;
1881 cflf[23] = -2.3773901;
1882 cfle[24] = 441.241363;
1883 cflf[24] = -2.5326197;
1884 cfle[25] = 514.7816162;
1885 cflf[25] = -2.5292389;
1886 cfle[26] = 588.3218384;
1887 cflf[26] = -2.8230250;
1888 cfle[27] = 661.8620605;
1889 cflf[27] = -2.9502323;
1890 cfle[28] = 735.4022827;
1891 cflf[28] = -3.0774822;
1892 cfle[29] = 1470.8045654;
1893 cflf[29] = -4.2239799;
1894 cfle[30] = 2206.2067871;
1895 cflf[30] = -5.2547927;
1896 cfle[31] = 2941.6091309;
1897 cflf[31] = -6.2353640;
1898 cfle[32] = 3677.0114746;
1899 cflf[32] = -7.1898708;
1900 cfle[33] = 4412.4135742;
1901 cflf[33] = -8.1292381;
1902 cfle[34] = 5147.8159180;
1903 cflf[34] = -9.0594845;
1904 cfle[35] = 5883.2182617;
1905 cflf[35] = -9.9830370;
1906 cfle[36] = 6618.6206055;
1907 cflf[36] = -10.9028034;
1908 cfle[37] = 7354.0229492;
1909 cflf[37] = -11.8188877;
1910 cfle[38] = 7400.0000000;
1911 cflf[38] = -30.0000000;
1912 cfle[39] = 10000000.0000000;
1913 cflf[39] = -30.0000000;
2009 fprintf(
ioQQQ ,
"lines_table checking lines within data table %s\n",
chLINE_LIST );
2015 double relative , absolute;
2019 fprintf(
ioQQQ,
"lines_table in parse_table.cpp did not find line %4s ",
chLabel[n]);
2021 fprintf(
ioQQQ,
"\n");
2027 fprintf(
ioQQQ ,
" BOTCHED ASSERTS!!! Botched Asserts!!! lines_table could not find a total of %li lines\n\n", miss );
2031 fprintf(
ioQQQ ,
"lines_table found all lines\n\n" );
2036 #if defined(__HP_aCC)
2037 #pragma OPTIMIZE OFF