Cloudy & Associates

Commit 550bf0c0 authored by Ferland, Gary's avatar Ferland, Gary
Browse files

improve diagnostic messages

parent d08d16d5
......@@ -64,13 +64,17 @@ STATIC void RT_line_pumping(
if( (t.Lo()->nelem()-1)!=ipHYDROGEN &&
t.EnergyWN()>=hydro.EnerLyaProf2 &&
t.EnergyWN()<=hydro.EnerLyaProf3)
{
OccNumContTotal += hydro.PhotOccNumLyaCenter;
//fprintf(ioQQQ,"DEBUGGG2 %d %d %d \n",
// (t.Lo()->nelem()-1) , (t.Lo()->IonStg()) , ipIRON);
}
if( trace.lgLyaPump )
{
fprintf(ioQQQ,"DEBUGGG %d %d occNum %.2e p1 %.2e",
(t.Lo()->nelem()-1) , ipHYDROGEN , hydro.PhotOccNumLyaCenter,
(t.Lo()->nelem()-1) , (t.Lo()->IonStg()) , hydro.PhotOccNumLyaCenter,
t.Emis().pump() );
}
......
......@@ -1099,6 +1099,19 @@ STATIC void setXtraRatesFe2(const TransitionProxy& tr, double &xtraExRate, doubl
{
DEBUG_ENTRY( "setXtraRatesFe2()" );
STATIC bool lgPrintFirstTime = true;
if( trace.lgLyaPump && lgPrintFirstTime )
{
fprintf(ioQQQ, "DEBUG setXtraRatesFe2 trace.lgLyaPump %d\n"
"hydro.lgLyaPumpAllLines = %d\n"
"hydro.lgLyaFeIIPumpOn %d\n",
trace.lgLyaPump,
hydro.lgLyaPumpAllLines ,
hydro.lgLyaFeIIPumpOn );
lgPrintFirstTime = false;
}
if( ! hydro.lgLyaFeIIPumpOn )
return;
......@@ -1135,7 +1148,6 @@ STATIC void setXtraRatesFe2(const TransitionProxy& tr, double &xtraExRate, doubl
else
hydro.PhotOccNumLyaCenter = 0.;
/* on first iteration optical depth in line is inward only, on later
* iterations is total optical depth */
double taux = 0.;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment