Changeset 2043 for branches/newmole/source/cool_neon.cpp

Timestamp:

05/11/08 14:27:43 (8 months ago)

Author:

rjrw

Message:

Integrate Nick Abel's Ne, Ar additions to the chemistry.

Files:

• branches/newmole/source/cool_neon.cpp (modified) (3 diffs)

branches/newmole/source/cool_neon.cpp

@@ -5 +5 @@
 #include "coolheavy.h"
 #include "taulines.h"
+#include "mole.h"
 #include "dense.h"
 #include "phycon.h"
@@ -24 +25 @@
           cs12,
           cs13,
-          cs23; 
+          cs23,
+          /* >>chng 08 may 06, Add collisional excitation of NeII by H, important in AGN */
+          cs_neII_h; 
         realnum 
           cs, 
@@ -43 +46 @@
         cs=(realnum)(0.132*phycon.te07*phycon.te02 *phycon.te004 *phycon.te0001 );
 
-        PutCS(cs,&TauLines[ipTNe13]);
+        /* >>chng 08 mag 06, Collisional de-excitation rate of fine-structure level of Ne+ by H impact, from *
+         * >>refer      ne2     csh Hollenbach, D. & McKee, C. F. 1989, ApJ, 342, 306 */
+
+        cs_neII_h = 1.3e-9;
+
+        PutCS(cs+cs_neII_h*(dense.xIonDense[ipHYDROGEN][0]+findspecies("H2")->den+findspecies("H2*")->den)/dense.cdsqte,&TauLines[ipTNe13]);
         atom_level2(&TauLines[ipTNe13]);