Changeset 1835 for branches/newmole/source/mole_solve.cpp

Timestamp:

03/09/08 10:35:33 (10 months ago)

Author:

rjrw

Message:

Move molecular source and sink copy to mole_effects routine.

Increase tolerance for negative density solutions somewhat.

Files:

• branches/newmole/source/mole_solve.cpp (modified) (2 diffs)

branches/newmole/source/mole_solve.cpp

@@ -39 +39 @@
 double mole_solve()
 {
-        long int nelem;
         int nBad, nPrevBad, i;
         realnum
@@ -121 +120 @@
         
         check_co_ion_converge();
-
-        for (nelem=ipHYDROGEN;nelem<LIMELM;nelem++)
-        {
-                for (long int ion=0;ion<N_MOLE_ION;ion++) 
-                {
-                        if (element_list[nelem]->ipMl[ion] != -1)
-                        {
-                                mole.source[nelem][ion] = mole.zone[element_list[nelem]->ipMl[ion]].src;
-                                mole.sink[nelem][ion] = mole.zone[element_list[nelem]->ipMl[ion]].snk;
-                                if (dense.IonLow[nelem] > ion)
-                                        dense.IonLow[nelem] = ion;
-                                if (dense.IonHigh[nelem] < ion)
-                                        dense.IonHigh[nelem] = ion;                                     
-                        }
-                        else
-                        {
-                                mole.source[nelem][ion] = 0.0;
-                                mole.sink[nelem][ion] = 0.0;
-                        }
-                }
-        }
 
         {