Changeset 1950 for trunk/source/ion_silic.cpp

Timestamp:

04/18/08 14:22:02 (9 months ago)

Author:

gary

Message:

ion_silic.cpp - comment out code that damped changed in valence shell photo rate - why was it there?

silicon - clean out dead code ion ion_silic.cpp -
sil.h - remove variable used to get pump rate
prt_lines_lv1_na_ar.cpp - use shell photo rate to get same pump rate

Files:

• trunk/source/ion_silic.cpp (modified) (4 diffs)

trunk/source/ion_silic.cpp

@@ -39 +39 @@
         static long nzUsed = -1;
         static double OldRate = 0.;
-        double save_rec;
 
         DEBUG_ENTRY( "IonSilic()" );
 
         /* silicon, atomic number 14 */
-        if( !dense.lgElmtOn[ipSILICON] )
-        {
-                sil.p1895 = 0.;
-                return;
-        }
 
         ion_zero(ipSILICON);
 
         ion_photo(ipSILICON,false);
+
+#       if 0
+        /* 2008 apr 18, this appears to try to damp out changes in the valence
+         * shell photo rate - perhaps due to Lya oscillations in dust free
+         * environments? */
         if( nzone > 1 && OldRate > 0. )
         {
@@ -67 +66 @@
         }
         OldRate = ionbal.PhotoRate_Shell[ipSILICON][0][4][0];
+#       endif
 
         /* find collisional ionization rates */
@@ -78 +78 @@
         /*lint +e64 type mismatch */
 
-        /* pumping to 3p of Si III */
-        sil.p1895 = ionbal.PhotoRate_Shell[ipSILICON][1][2][0];
-
-        /* >>chng 03 sep 29, synch up ion and co solvers.
-         * the co solver will have a different C/C+ balance that
-         * is strongly affected by the chemistryif we are in neutral gas
-         * in this case use co solver's C/C+ balance */
-        save_rec = ionbal.RateRecomTot[ipSILICON][0];
-        /*>>chng 04 apr 27, do not test for ionhigh being 1,
-         * no reason for test on upper stage of ionization,
-         * if codrive called but not evaluted then hevmol is all zero */
-        /* >>chng 04 sep 10, rm check on search phase, no reason for it */
-#       if 0
-        if( dense.IonLow[ipSILICON]==0  && /*!conv.lgSearch &&*/
-                co.hevmol[ipSIP] > SMALLFLOAT && 
-                ionbal.RateIonizTot[ipSILICON][0]*co.hevmol[ipATSI]>0.)
-        {
-                /* >>chng 04 sep 11, use damped ratio, were Si+ oscillations deep in molecular cloud,
-                 * where Si+ is dominant electron donor */
-                ionbal.RateRecomTot[ipSILICON][0] = 
-                        ionbal.RateIonizTot[ipSILICON][0]*
-                        co.hevmol[ipATSI]/co.hevmol[ipSIP];
-                /*fprintf(ioQQQ,"DEBUG sil atm\t%.2f\t%.3e\t%.3e\t%.3e\n",
-                        fnzone , 
-                        co.hevmol[ipATSI],
-                        co.hevmol[ipSIP],
-                        ionbal.RateIonizTot[ipSILICON][0]);*/
-        }
-#       endif
-
         /* solve for ionization balance */
         /*if( nzone>700 ) lgPrtDebug = true;*/
         ion_solver(ipSILICON,lgPrtDebug);
 
-        /* reset the var we just hosed */
-        if( save_rec > 0. )
-                ionbal.RateRecomTot[ipSILICON][0] = save_rec;   
-
-        /* get pmp rate cm-3 s-1 */
-        sil.p1895 *= dense.xIonDense[ipSILICON][1]*0.85;
-
+#       if 0
+        /* 2008 apr 18, this is all dead code from when co network fed back from
+         * the ionization */
         if( trace.lgTrace && trace.lgHeavyBug )
         {
@@ -128 +94 @@
                 fprintf( ioQQQ, "\n" );
         }
-#       if 0
         fprintf( ioQQQ, "DEBUG silicon\t%.2f",fnzone );
         for( i=0; i < 2; i++ )