Changeset 1599

Timestamp:

12/01/07 11:54:54 (7 months ago)

Author:

rjrw

Message:

Diagnose poor ion convergence

Comment out redundant set of xIonDense

Files:

• branches/newmole/source/conv_base.cpp (modified) (5 diffs)
• branches/newmole/source/mole_co_solve.cpp (modified) (1 diff)

branches/newmole/source/conv_base.cpp

@@ -378 +378 @@
-        {
-
+                fprintf(ioQQQ,"Start of loop %g\n",dense.xIonDense[15][0]);
-                check_conserve(__FILE__,__LINE__);
-        
@@ -412 +413 @@
-                /*fprintf(ioQQQ,"DEBUG h2\t%.2f\t%.3e\t%.3e", fnzone,hmi.H2_total,findspecies("H2")->hevmol);*/
-                iso_solve( ipH_LIKE );
+                fprintf(ioQQQ,"After H_LIKE %g\n",dense.xIonDense[15][0]);
-
-                /* now do level populations for H2 */
@@ -432 +434 @@
-                /* do all he-like species */
-                iso_solve( ipHE_LIKE );
+                fprintf(ioQQQ,"After He_like %g\n",dense.xIonDense[15][0]);
-
-                /*>>chng 04 may 09, add option to abort here, inspired by H2 pop failures
@@ -506 +509 @@
-                IonOxyge();
-                IonNitro();
+                fprintf(ioQQQ,"Before Si %g\n",dense.xIonDense[15][0]);
-                IonSilic();
+                fprintf(ioQQQ,"Before S %g\n",dense.xIonDense[15][0]);
-                IonSulph();
+                fprintf(ioQQQ,"After ion %g\n",dense.xIonDense[15][0]);
-                IonChlor();
-                check_conserve(__FILE__,__LINE__);
@@ -513 +519 @@
-                /* do carbon monoxide after oxygen */
-                CO_drive();
+                fprintf(ioQQQ,"After mole %g\n",dense.xIonDense[15][0]);
-
-                check_conserve(__FILE__,__LINE__);

branches/newmole/source/mole_co_solve.cpp

@@ -73 +73 @@
-                                mole.source[nelem][ion] = mole.src[element_list[nelem]->ipMl[ion]];
-                                mole.sink[nelem][ion] = mole.snk[element_list[nelem]->ipMl[ion]];
-
+// 
-                        }
-                        else