prt_linesum.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*PrtLineSum parse print line sum command to enter set of lines into sum  */
#include "cddefines.h"
#include "cddrive.h"
#include "radius.h"
#include "lines.h"
#include "input.h"
#include "parser.h"
/* this is the limit to the number of lines we can save */
#define NRDSUM  300L
#include "prt.h"
 
static char **chSMLab;
static long int *ipLine;
static bool lgFirst=true;
static long nlsum;
static realnum *wavelength;
 
void ParsePrtLineSum(Parser &p)
{
        /*static char chSMLab[NRDSUM][5];*/
        bool lgEND;
        long int i;
 
        /* remember whether we have been called before */
 
        DEBUG_ENTRY( "ParsePrtLineSum()" );
 
#               if 0
        if( !lgFirst )
        {
                /* error - more than one read in input stream */
                fprintf(ioQQQ," more than one print line sum has appeared - only first one is used.\n");
                fprintf(ioQQQ," Sorry.\n");
                cdEXIT(EXIT_FAILURE);
        }
        else
#               endif
                /* >>chng 03 jan 23, if not first call, do not allocate space, 
                 * had aborted, which was bad in optized runs, or in a grid. 
                 * Bug caught by Melekh Bohdan */
                if( lgFirst )
                {
                        /* do not malloc space again */
                        lgFirst = false;
                        wavelength = ((realnum *)MALLOC( sizeof(realnum )*NRDSUM ));
                        
                        /* create space for the array of array indices for lines*/
                        ipLine = ((long int *)MALLOC(NRDSUM*sizeof(long)));
                        
                        /* create space for the array of labels*/
                        chSMLab = ((char **)MALLOC(NRDSUM*sizeof(char *)));
                        
                        for( i=0; i<NRDSUM; ++i )
                        {
                                chSMLab[i] = ((char *)MALLOC(5*sizeof(char )));
                        }
                }
        
        /* now read in lines */
        nlsum = 0;
        lgEND = false;
        while( !lgEND )
        {
                p.getline();
                if( p.m_lgEOF )
                {
                        fprintf( ioQQQ, " Hit EOF while reading line list; use END to end list.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
                
                if( p.strcmp("END" ) != 0 )
                {
                        if( nlsum >= NRDSUM )
                        {
                                fprintf( ioQQQ, 
                                                        " Too many lines have been entered; the limit is %li.  Increase NRDSUM in PrtLineSum.\n", 
                                                        NRDSUM );
                                cdEXIT(EXIT_FAILURE);
                        }
                        
                        /*  order on line is label (col 1-4), wavelength */
                        strncpy( chSMLab[nlsum], p.getCommand(4).c_str() , 4 );
                        chSMLab[nlsum][4] = 0;
                        wavelength[nlsum] = (realnum)p.getWaveOpt();
                        ++nlsum;
                }
                else
                {
                        lgEND = true;
                }
        }
}
double PrtLineSum(void)
{
        long int i;
 
        /* remember whether we have been called before */
 
        double absint, 
          relint ,
          sum=-1.;
 
        DEBUG_ENTRY( "PrtLineSum()" );
 
        sum = 0.;
        /* this can be called during setup mode, in which case we do nothing */
        if( LineSave.ipass <= 0 )
        { 
                return sum;
        }
        
        if( nzone == 1 )
        {
                for( i=0; i < nlsum; i++ )
                {
                        /* save the array index for each line */
                        if( (ipLine[i] = cdLine((char*)chSMLab[i],wavelength[i],&relint,&absint) ) <=0 )
                        {
                                fprintf( ioQQQ, " PrtLineSum could not fine line %4.4s %5f\n", 
                                                        chSMLab[i], wavelength[i] );
                                cdEXIT(EXIT_FAILURE);
                        }
                }
        }
        
        /* now sum the line */
        for( i=0; i < nlsum; i++ )
        {
                /* this version of chLine uses index, does not search*/
                cdLine_ip(ipLine[i],&relint,&absint);
                absint = pow(10.,absint - radius.Conv2PrtInten);
                sum += absint;
        }
        return sum;
}