Changeset 2296 for trunk/source/iso_ionize_recombine.cpp

Timestamp:

08/01/08 13:48:31 (5 months ago)

Author:

gary

Message:

series of small changes resulting from a few hours debugging the negative population bugs now current in trunk. these should have little effect

cont_setintensity.cpp enforce use of EdenSet? when used

eden_sum - tidy up prints, assert EdenSet? used at exit

iso_ionize_recombine.cpp - on very first pass through routine set initial recomb rate

iso_level.cpp - small changes in print statements

Files:

• trunk/source/iso_ionize_recombine.cpp (modified) (2 diffs)

trunk/source/iso_ionize_recombine.cpp

@@ -5 +5 @@
 #include "cddefines.h"
 #include "ionbal.h"
+#include "conv.h"
 #include "atmdat.h"
 #include "dense.h"
@@ -246 +247 @@
         iso.RecomCollisFrac[ipISO][nelem] = Recom3Body* dense.eden / ionbal.RateRecomTot[nelem][nelem-ipISO];
 
+        /* very first pass through here rate RateIonizTot not yet evaluated */
+        if( conv.nTotalIoniz==0 )
+                ionbal.RateIonizTot[ipISO][nelem] = iso.RateLevel2Cont[ipISO][nelem][0];
+
         /* get simple estimate of atom to ion ratio */
         if( ionbal.RateRecomTot[nelem][nelem-ipISO] > 0. )