Changeset 2134

Timestamp:

06/28/08 02:34:32 (5 months ago)

Author:

gary

Message:

radius_next.cpp - bugfix - logic for zoning time dependent sim assumed first zone was zero

Files:

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

trunk/source/radius_next.cpp

@@ -1036 +1036 @@
         }
 
-        /*>>chng 06 apr 03, add this
-         * this is special case of time dependent recombination - conditions become
-         * very homogenious and tend to crash into I fronts - use last iteration's 
-         * radius as guess of current case */
+        /* time dependent recombination - conditions become very homogeneous and 
+         * crash into I fronts - use last iteration's  radius as guess of current case */
         recom_dr_last_iter = BigRadius;
         if( dynamics.lgStatic && dynamics.lgRecom )
         {
+                /* first time through nzone == 1 */
                 static long int nzone_recom = -1;
-                if( nzone==0 )
+                if( nzone<=1 )
                 {
                         /* initialization case */
                         nzone_recom = 0;
                 }
-                else if( radius.depth < struc.depth_last[struc.nzone_last-1] )
-                {
+                else if( radius.depth < struc.depth_last[struc.nzone_last-1] && 
+                        nzone_recom < struc.nzone_last )
+                {
+                        ASSERT( nzone_recom>=0 && nzone_recom<struc.nzone_last );
                         /* case where we are within previous computed structure 
                          * first possibly increase nzone_recom so that it points deeper