prt_columns.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2022 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*PrtColumns print column densities of all elements */
#include "cddefines.h"
#include "cddrive.h"
#include "dense.h"
#include "elementnames.h"
#include "h2.h"
#include "taulines.h"
#include "molcol.h"
#include "generic_state.h"
#include "prt.h"

void PrtColumns(
        /* this is stream used for io, is stdout when called by final,
         * is save unit when save output generated */
         FILE *ioMEAN  )
{
        DEBUG_ENTRY( "PrtColumns()" );

        /* print molecular element column densities */
        molcol("PRIN" , ioMEAN);
        
        fprintf( ioMEAN, "\n" );
        
        fprintf( ioMEAN, "\n         " );
        for( long i=1; i <= 17; i++ )
        {
                fprintf( ioMEAN, "%7ld", i );
        }
        fprintf( ioMEAN, "\n\n" );
        
        /* ionization column densities  */
        for( long nelem=0; nelem < LIMELM; nelem++ )
        {
                if( dense.lgElmtOn[nelem] )
                {
                        bool lgDONE = false;
                        
                        fprintf( ioMEAN, " %10.10s", elementnames.chElementName[nelem]  );
                        
                        long i = 1;
                        while( !lgDONE )
                        {
                                double aa;
                                if( cdColm(
                                                 /* return value is zero if all ok, 1 if errors happened */
                                                 /* 4-char + eol string that is first
                                                  * 4 char of element name as spelled by cloudy */
                                                 elementnames.chElementNameShort[nelem],
                                                 
                                                 /* integer stage of ionization, 1 for atom, 0 for CO or H2 */
                                                 i,
                                                 
                                                 /* the theoretical column density derived by the code */
                                                 &aa ) )
                                        TotalInsanity();
                                
                                if( aa == 0. )
                                {
                                        aa = -30.;
                                }
                                else if( aa > 0. )
                                {
                                        aa = log10(aa);
                                }
                                
                                if( i == 18 )
                                {
                                        fprintf( ioMEAN, "\n" );
                                }
                                fprintf( ioMEAN, "%7.3f", aa );
                                
                                /* increment counter and check if at upper limit */
                                ++i;
                                /* MAX2 is to include H2 in H array */
                                if( i > MAX2(3,nelem+2) )
                                        lgDONE = true;
                                
                                /* print title for this info if we are done with hydrogen */
                                if( nelem==ipHYDROGEN && lgDONE )
                                        fprintf(ioMEAN," (H2)                Log10 Column density (cm^-2)");
                        }
                        
                        fprintf( ioMEAN, "\n" );
                }
        }
        
        /* only print excited state column densities if level2 lines are included
         * since they populated the upper level by UV pumping.   This process
         * is not included if level2 lines are not considered, by introducing
         * the "no level2" command */
        if( nWindLine>0 )
        {
                const int nEXCIT_COL = 12;
                // use strings instead of char[][] to avoid bogus errors with PGCC bounds checking
                string strExcit_Col[nEXCIT_COL]={
                        "He[2]","C+[2]","C[1]","C[2]","C[3]","O[1]","O[2]","O[3]","Si+[1]","C+2[2]","C+2[3]","C+2[4]"};
                /* print  excited level column densities */
                fprintf(ioMEAN," Exc state ");
                long int nprt = 12;
                for(long i=0; i<nEXCIT_COL; ++i )
                {
                        ASSERT( strExcit_Col[i].size() <= NCHLAB-1 );
                        vector<genericState> v = matchGeneric(strExcit_Col[i].c_str(),false);
                        ASSERT ( v.size() <= 1 );
                        genericState gs= v.size() == 1 ? v[0] : genericState();
                        
                        if( nprt > 120 )
                        {
                                fprintf(ioMEAN,"\n           ");
                                nprt = 12;
                        }
                        fprintf(ioMEAN,"   %*.*s%7.3f",
                                          NCHLAB-1,NCHLAB-1,strExcit_Col[i].c_str(),
                                          log10(SDIV(column(gs)) ));
                        nprt += 10+(NCHLAB-1);
                }
                fprintf(ioMEAN,"\n");
        }
        
        /* print column densities for H2 */
        h2.H2_Prt_column_density(ioMEAN);
        
        fprintf(ioMEAN,"\n");
}