parse_rangeoption.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 */
/*ParseRangeOption parse the range option on the luminosity command */
#include "cddefines.h"
#include "rfield.h"
#include "parser.h"

void ParseRangeOption(
        /* the input line image */
        Parser &p)
{
        bool lgP1Absent, lgP2Absent;
        double p1, 
          p2;

        DEBUG_ENTRY( "ParseRangeOption()" );

        if( p.nMatch("TOTA") )
        {
                rfield.range[p.m_nqh][0] = rfield.emm();
                rfield.range[p.m_nqh][1] = rfield.egamry();
        }
        else if( p.nMatch("RANG") )
        {
                p.set_point(0);
                /* first call return the luminosity on the card, ignore it */
                (void) p.FFmtRead();
                (void) p.lgEOL();

                /* read lower limit and upper limit*/
                p1 = p.FFmtRead();
                lgP1Absent = p.lgEOL();
                p2 = p.FFmtRead();
                lgP2Absent = p.lgEOL();

                /* option to enter log if first energy is neg */
                if( p1 < 0. || p.nMatch(" LOG") )
                {
                        p1 = exp10(p1);
                        p2 = exp10(p2);
                }

                if( lgP1Absent )
                        p1 = rfield.emm();
                if( lgP2Absent )
                        p2 = rfield.egamry();

                /* make sure that energies are within array bounds */
                rfield.range[p.m_nqh][0] = MAX2(p1,rfield.emm());
                rfield.range[p.m_nqh][1] = MIN2(p2,rfield.egamry());
                if( rfield.range[p.m_nqh][0] >= rfield.range[p.m_nqh][1] )
                {
                        fprintf( ioQQQ, " Range MUST be in increasing order - sorry.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
        }
        else
        {
                /* range not specified, use default - total ionizing field
                 * >>chng 96 dec 18, from 1.001 to 1 Ryd for H mass nuc */
                rfield.range[p.m_nqh][0] = HIONPOT;
                rfield.range[p.m_nqh][1] = rfield.egamry();
        }
        return;
}