Changeset 1797 for branches/newmole/source

Timestamp:

02/13/08 14:15:39 (9 months ago)

Author:

rjrw

Message:

Add molecular transfer rates to iso_sequence

Multiply rates by abundance in molecular network balance debug print,
to make interpretation clearer.

Location:

branches/newmole/source

Files:

• iso_ionize_recombine.cpp (modified) (1 diff)
• iso_level.cpp (modified) (2 diffs)
• mole_eval_balance.cpp (modified) (1 diff)

branches/newmole/source/iso_ionize_recombine.cpp

@@ -183 +183 @@
         iso.RateCont2Level[ipISO][nelem][level] += gv.GrainChTrRate[nelem][nelem+1-ipISO][nelem-ipISO];
 
+        for(ion=0; ion<nelem+1; ++ion)
+                if ( ion < nelem-ipISO )
+                        iso.RateLevel2Cont[ipISO][nelem][level] += mole.xMoleChTrRate[nelem][nelem-ipISO][ion];
+                else if ( ion > nelem-ipISO )
+                        iso.RateCont2Level[ipISO][nelem][level] += mole.xMoleChTrRate[nelem][ion][nelem-ipISO];
+
         /* now charge transfer - all into/from ground, two cases, H and not H */
         if( nelem==ipHYDROGEN )

branches/newmole/source/iso_level.cpp

@@ -406 +406 @@
                 }
 
-#if     0
                 /* add in source and sink terms from molecular network. */
-                CodeReview(); /* Check... */
                 if( nelem == ipISO && conv.nTotalIoniz )
                 {
@@ -420 +418 @@
                         }
                 }
-#endif
 
                 /* >>chng 04 nov 30, atom XX-like collisions off turns this off too */

branches/newmole/source/mole_eval_balance.cpp

@@ -205 +205 @@
                 {
                         fprintf(stdout,"%20.20s snk %13.7g of %13.7g [",
-                                                        ratesnk->label,snkx,mole.zone[debug_species->index].snk);
+                                                        ratesnk->label,snkx*mole.list[debug_species->index]->den,
+                                                        mole.zone[debug_species->index].snk*mole.list[debug_species->index]->den);
                         for (j=0;j<ratesnk->nreactants;j++) 
                         {