67 fprintf(
ioQQQ,
" lines called\n" );
88 for( nelem=0; nelem<
LIMELM; ++nelem )
127 linadd(0.f,0,
"zero",
'i' ,
"null placeholder");
132 linadd( 1.e-10 , 1 ,
"Unit" ,
'i' ,
"unit integration placeholder");
133 static long int ipOneAng=-1;
136 lindst( 1.e-10 , 1. ,
"UntD" , ipOneAng ,
'i' ,
false,
"unit integration placeholder");
163 fprintf(
ioQQQ,
"ipHi\tipLo\tnu\tlu\tsu\tnl\tll\tsl\tWL\tintens\n" );
166 for(
long ipLo=0; ipLo<ipHi; ipLo++ )
170 double relint, absint;
180 ipLo!=3 && ipLo!=4 && relint >= 0.0009 )
182 fprintf(
ioQQQ,
"%li\t%li\t%li\t%li\t%li\t%li\t%li\t%li\t%e\t%e\n",
214 for (
int ipSpecies=0; ipSpecies <
nSpecies; ++ipSpecies)
217 em !=
dBaseTrans[ipSpecies].Emis().end(); ++em)
220 if( (*em).Tran().ipCont() > 0)
222 PutLine((*em).Tran(),
"lines from third-party databases", (*(*em).Tran().Hi()).chLabel());
247 "Stoy method energy sum ");
252 "recombination lines");
258 for( i=0; i < 471; i++ )
267 LineSave.RecCoefCNO[1][i]+1.01)-1] );
281 if(
dense.eden < 1e8 )
283 nelem = (long)
LineSave.RecCoefCNO[0][i]-1;
284 long int ion = (long)(
LineSave.RecCoefCNO[0][i]-
LineSave.RecCoefCNO[1][i]+2)-1;
286 dense.xIonDense[nelem][ion];
289 f2 = f2*1.99e-8/
LineSave.RecCoefCNO[2][i];
298 "recombination line");
315 if(
TauLine2[i].Coll().cool() > BigstExtra )
317 BigstExtra =
TauLine2[i].Coll().cool();
320 ExtraCool +=
TauLine2[i].Coll().cool();
329 "hyperfine structure lines ");
333 "total cooling all hyperfine structure lines ");
342 "hyperfine structure line");
348 "inner shell lines");
372 fprintf(
ioQQQ,
"DEBUGyeild\t%s\t%.3f\t%.3e\n",
398 static long nLineSave=-1 , ndLineSave=-1;
407 if( nLineSave<= 0 || ndLineSave < 0 )
415 fprintf(
ioQQQ,
"DISASTER number of lines in LineSave.nsum changed between pass 0 and 1 - this is impossible\n" );
416 fprintf(
ioQQQ,
"DISASTER LineSave.nsum is %li and nLineSave is %li\n",
424 fprintf(
ioQQQ,
"DISASTER number of lines in LineSave.nsum changed between pass 0 and 1 - this is impossible\n" );
425 fprintf(
ioQQQ,
"DISASTER LineSave.nsum is %li and ndLineSave is %li\n",
439 fprintf(
ioQQQ,
" lines returns\n" );
460 fprintf(
ioQQQ,
"n\tlab\twl\n");
467 fprintf(
ioQQQ,
" Thanks for checking on the cdLine routine!\n" );
NORETURN void TotalInsanity(void)
#define DEBUG_ENTRY(funcname)
long int cdLine(const char *chLabel, realnum wavelength, double *relint, double *absint)
EmissionProxy::iterator iterator
void rec_lines(double t, realnum r[][471])
realnum yield(long n) const
realnum energy(long n) const
long ipoint(double energy_ryd)
t_elementnames elementnames
t_iso_sp iso_sp[NISO][LIMELM]
void iso_satellite_update(long nelem)
long int StuffComment(const char *chComment)
void lines_molecules(void)
void lines_lv1_li_ne(void)
void lines_lv1_k_zn(void)
void lines_continuum(void)
void lines_lv1_na_ar(void)
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
void lindst(double xInten, realnum wavelength, const char *chLab, long int ipnt, char chInfo, bool lgOutToo, const char *chComment)
static realnum * wavelength
UNUSED const double EN1RYD
UNUSED const double EVRYD
UNUSED const double RYDLAM
STATIC void Drive_cdLine(void)
static bool lgMustPrintHeader
TransitionList TauLine2("TauLine2", &AnonStates)
vector< TransitionList > dBaseTrans
TransitionList HFLines("HFLines", &AnonStates)
void PutExtra(double Extra)
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)