Changeset 2045 for trunk/source/rt_diffuse.cpp

Timestamp:

05/12/08 03:16:01 (8 months ago)

Author:

gary

Message:

rt_diffuse.cpp - BUGFIX - in case where extremely large increase in ionization threshold continuum index for ionization edge could extend beyond the number of cells

ion_trim.cpp - when the density falls to VERY low values, too low to compute on 32 bit cpu, error message is printed and code aborts. Improve information in error message.

service.cpp - fix misspelled words, one use in output

init_coreload.cpp - trivial changes in line skips

Files:

• trunk/source/rt_diffuse.cpp (modified) (1 diff)

trunk/source/rt_diffuse.cpp

@@ -384 +384 @@
 
                                 ip = Heavy.ipHeavy[nelem][ion]-1;
+                                ASSERT( ip >= 0 );
 
                                 /* nflux was reset upward in ConvInitSolution to encompass all
                                  * possible line and continuum emission.  this test should not
-                                 * possibly fail.  It could if the ionization were to increase with depth.
-                                 * This is important because the nflux cell in ConInterOut is used to carry out the
-                                 * unit integration, and if it gets clobbered by diffuse emission the code
-                                 * will declare insanity in PrtComment */
-                                ASSERT( ip >= 0 && ip < rfield.nflux );
+                                 * possibly fail.  It could if the ionization were to increase with depth
+                                 * although the continuum mesh is designed to deal with this.
+                                 * This test is important because the nflux cell in ConInterOut 
+                                 * is used to carry out the unit integration, and if it gets 
+                                 * clobbered by diffuse emission the code will declare 
+                                 * insanity in PrtComment */
+                                if( ip >= rfield.nflux )
+                                        continue;
 
                                 /* get shell number, stat weights for this species */