22 long int i, nelem, ipHi, ipLo;
34 fprintf(
ioQQQ,
" lines_hydro called\n" );
39 " start H -like iso sequence ");
42 " total collisional cooling due to all hydrogen lines ");
45 " total collisional heating due to all hydrogen lines ");
53 "collisionally excited La cooling ");
56 " collisionally de-excited La heating ");
59 " cooling due to n>2 Lyman lines ");
62 " heating due to n>2 Lyman lines ");
65 " cooling due to n>3 Balmer lines ");
68 " heating due to n>3 Balmer lines ");
71 " cooling due to higher Paschen lines ");
74 " heating due to higher Paschen lines ");
87 "Lyaa contribution from suprathermal secondaries from ground ");
94 linadd(pump,4861,
"Pump",
'r',
95 "part of Hbeta formed by continuum pumping");
98 "collision ionization cooling of hydrogen ");
101 " this is the heating due to 3-body recombination ");
105 " Stark broadening contribution to line ");
109 " Stark broadening contribution to line ");
113 "Stark broadening contribution to line ");
117 " Stark broadening contribution to line ");
125 "Stark broadening part of line");
133 " portion of line lost due to absorp by background opacity ");
140 "Ha destroyed by background opacity");
147 "portion of line lost due to absorb by background opacity");
154 "portion of line lost due to absorb by background opacity");
161 "portion of line lost due to absorb by background opacity");
165 "Ly-alpha destroyed by overlap with FeII " );
168 " simple high-density case b intensity of Ly-alpha, no two photon ");
191 "Case B H-beta computed from Q(H) and specified covering factor");
204 "Ly-alpha from Q(H), high-dens lim, specified covering factor" );
207 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
213 for( ipLo=0; ipLo < ipHi; ipLo++ )
219 Transitions[ipISO][nelem][ipHi][ipLo].Emis->phots =
226 Transitions[ipISO][nelem][ipHi][ipLo].Emis->xIntensity =
236 for( nelem=0; nelem <
LIMELM; nelem++ )
247 index_of_nHi_P = ipHi;
252 for( ipLo=0; ipLo < ipHi; ipLo++ )
265 if( ipHi == index_of_nHi_P && ipLo == index_of_nLo_S )
291 " this is old case b, had been in code for very long time");
304 " old prediction of He II 1640, case B at low densities");
308 em = 2.52e-20/(pow(
phycon.
te,1.05881));
319 for(nelem=0; nelem<
HS_NZ; ++nelem )
326 for( nelem=0; nelem <
LIMELM; nelem++ )
333 if( nelem <
HS_NZ && (nelem<2 || nelem>4) )
336 for( iCase=0; iCase<2; ++iCase )
338 char chAB[2]={
'A',
'B'};
350 double case_b_Intensity;
351 long int ipCHi , ipCLo;
363 if( case_b_Intensity<=0. )
366 case_b_Intensity = 0.;
371 if( iCase==0 && ipLo==1 )
385 chLab[3] = chAB[iCase];
396 else if( ipCLo == 2 )
401 else if( ipCLo == 1 || ipCLo == 0 )
413 " case a or case b from Hummer & Storey tables" );
431 "two-photon emission");
439 "induced two-photon emission ");
457 char chCommentTemp[23];
459 strcat( chCommentTemp,
" - " );
461 const char* chComment = chCommentTemp;
465 "predicted line, all processes included");
475 fprintf(
ioQQQ,
" lines_helium returns\n" );