Changeset 1913 for trunk/source/prt_comment.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_comment.cpp (modified) (2 diffs)

trunk/source/prt_comment.cpp

@@ -834 +834 @@
         }
 
-        if( !iso.lgCritDensLMix[ipHE_LIKE] && dense.lgElmtOn[ipHELIUM] )
-        {
-                sprintf( chLine,
-                        "   Critical density for l-mixing of He I not reached. "
-                        " More resolved levels are needed for accurate He I line ratios.");
-                notein(chLine);
+        /* check for critical density for l-mixing */
+        for( ipISO=ipH_LIKE; ipISO<NISO; ++ipISO )
+        {
+                if( !iso.lgCritDensLMix[ipISO] && dense.lgElmtOn[ipISO] )
+                {
+                        sprintf( chLine,
+                                "   The density is too low to l-mix the lowest %s I collapsed level. "
+                                " More resolved levels are needed for accurate line ratios.",
+                                elementnames.chElementSym[ipISO]);
+                        notein(chLine);
+                }
         }
 
@@ -849 +854 @@
                         if( iso.lgLevelsLowered[ipISO][nelem] && dense.lgElmtOn[nelem] )
                         {
-                                sprintf( chLine, "  !Continuum was lowered into model %s-like %s due to high density.  Max n is %li",
+                                sprintf( chLine, "  !Continuum was lowered into model %s-like %s due to high density.  Highest n is %li",
                                         elementnames.chElementSym[ipISO],
-                                        elementnames.chElementNameShort[nelem],
+                                        elementnames.chElementSym[nelem],
                                         iso.n_HighestResolved_local[ipISO][nelem]+iso.nCollapsed_local[ipISO][nelem]);
                                 bangin(chLine);