Changeset 1210

Timestamp:

06/24/07 10:03:57 (1 year ago)

Author:

rjrw

Message:

Include molecular terms in helike balance.

Files:

• branches/newmole/source/helike_level.cpp (modified) (1 diff)
• branches/newmole/source/hydrolevelpop.cpp (modified) (1 diff)

branches/newmole/source/helike_level.cpp

@@ -393 +393 @@
-                        /*iso.RateLevel2Cont[ipHE_LIKE][nelem][ipHe1s1S] += dynamics.Rate*dynamics.lgISO[ipHE_LIKE];*/
-                }
+                
+                CodeReview(); /* Check... */
+                if( nelem == ipHELIUM && (conv.nTotalIoniz > 1 || iteration > 1 ) )
+                {
+                        
+                        /* these are the external source and sink terms */
+                        /* source first */
+                        creation[ipHe1s1S] += mole.source[nelem][nelem-1]/SDIV(dense.xIonDense[nelem][nelem]);
+                        
+                        for( ipLo=ipHe1s1S; ipLo<iso.numLevels_local[ipHE_LIKE][nelem]; ++ipLo )
+                        {
+                                z[ipLo][ipLo] += mole.sink[nelem][nelem-1];
+                        }
+                }
-
-                /* error in how Hx12 is allocated, sec dim has only LIMELM elements,

branches/newmole/source/hydrolevelpop.cpp

@@ -333 +333 @@
-        }
-
- signs
-conv.nTotalIoniz > 1 || iteration > 1
+...
+nelem == ipHYDROGEN && (conv.nTotalIoniz > 1 || iteration > 1 )
-        {
-                
-                /* these are the external source and sink terms */
-                /* source first */
-i
+nelem
-                
-                for( ipLo=ipH1s; ipLo<iso.numLevels_local[ipH_LIKE][nelem]; ++ipLo )
-                {
-ipH_LIKE
+nelem
-                }
-        }