Changeset 1913 for trunk/source/prt_alltau.cpp

Timestamp:

04/03/08 12:44:50 (9 months ago)

Author:

gary

Message:

these changes affect the predicted iso sequence lines. there are three major parts:

cont_createmesh.cpp the fine continuum did not extend down to the wavelength of HeI 10830 so line overlap was not being treated. Extended low energy limit to below this line. this caused some asserted 10830 lines to change.

do not print H-like or He-like lines that are produced from collapsed levels. Their emission has problems which will be documented in a trac ticket.

the data/LineList* files had H and He like lines from collapsed levels commented out. The code now predicts far fewer lines from H and He like species.

iso_level - now keep track of whether the l-mixing collisions out of p are fast enough to mix the levels. if not a comment is printed saying that the H, He lines are not accurate

prt_alltau.cpp - do not print He I line optical depths twice

parse_atom_iso.cpp - add print levels option to atom H-like levels print or atom He-like levels print to print sizes of all ions of H and He sequence. the command does not change the number of levels.

Files:

• trunk/source/prt_alltau.cpp (modified) (2 diffs)

trunk/source/prt_alltau.cpp

@@ -116 +116 @@
         PrintE93( ioQQQ , opac.TauTotalGeo[0][iso.ipIsoLevNIonCon[ipH_LIKE][1][ipH1s]-1]/fcon);
         fprintf( ioQQQ, "\n");
-
-        if( dense.lgElmtOn[ipHELIUM] )
-        {
-                /* print optical depths of previous iteration with old atom */
-
-                /* these results are for new HeI as part of he iso sequence */
-                fprintf( ioQQQ, "  Line Optical Depths: 10830:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", (Transitions[ipHE_LIKE][ipHELIUM][ipHe2p3P0][ipHe2s3S].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe2p3P1][ipHe2s3S].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe2p3P2][ipHe2s3S].Emis->TauIn)/flin));
-                fprintf( ioQQQ, "  3889:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", Transitions[ipHE_LIKE][ipHELIUM][ipHe3p3P][ipHe2s3S].Emis->TauIn/flin));
-                fprintf( ioQQQ, "    5876:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", (Transitions[ipHE_LIKE][ipHELIUM][ipHe3d3D][ipHe2p3P0].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe3d3D][ipHe2p3P1].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe3d3D][ipHe2p3P2].Emis->TauIn)/flin));
-                fprintf( ioQQQ, " 7065:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", (Transitions[ipHE_LIKE][ipHELIUM][ipHe3s3S][ipHe2p3P0].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe3s3S][ipHe2p3P1].Emis->TauIn+
-                        Transitions[ipHE_LIKE][ipHELIUM][ipHe3s3S][ipHe2p3P2].Emis->TauIn)/flin));
-                fprintf( ioQQQ, "  2.06m:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", Transitions[ipHE_LIKE][ipHELIUM][ipHe2p1P][ipHe2s1S].Emis->TauIn/flin));
-                fprintf( ioQQQ, "  21c:");
-                fprintf( ioQQQ,PrintEfmt("%9.2e", HFLines[0].Emis->TauIn/flin));
-                /* the HeII 3 cm equivalent of H 21 cm */
-                fprintf( ioQQQ, "\n" );
-        }
 
         /* print optical depths of some metal lines */
@@ -251 +224 @@
                 /* ================================================================================*/
 
-
                 fprintf( ioQQQ, "          Old He Is Lines:" );
                 fprintf( ioQQQ, " %4d",584 );