parse_ratio.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 */
/*ParseRatio derive continuum luminosity of this continuum relative to previous */
#include "cddefines.h"
#include "optimize.h"
#include "input.h"
#include "rfield.h"
#include "parser.h"
/* HP cc cannot compile following except in -O1 mode */
#if defined(__HP_aCC)
#pragma OPT_LEVEL 1
#endif

void ParseRatio(Parser &p)
{
        bool lgAoxOn;
        double aox;

        DEBUG_ENTRY( "ParseRatio()" );

        /* enter a continuum luminosity as a ratio of
         * nuFnu for this continuum relative to a previous continuum
         * format; first number is ratio of second to first continuum
         * second number is energy for this ratio
         * if third numbewr on line, then 2nd number is energy of
         * first continuum, while 3rd number is energy of second continuum */

        if( p.m_nqh == 0 )
        {
                fprintf( ioQQQ, " Can\'t form ratio since this is first continuum.\n" );
                cdEXIT(EXIT_FAILURE);
        }

        /* per sq cm or luminosity is really irrelevant */
        strcpy( rfield.chRSpec[p.m_nqh], "SQCM" );
        strcpy( rfield.chSpNorm[p.m_nqh], "RATI" );

        /* this option is to specify alpha ox */
        if( p.nMatch("ALPH") )
        {
                /* lgAoxOn is flag saying that we will spicify alpha ox */
                lgAoxOn = true;
                /* only one parameter to be recognized, alpha ox */
                aox = p.FFmtRead();

                /* 403.3 is ratio of energies where alpha ox defined,
                 * assumed to be 2500A and 2keV */
                rfield.totpow[p.m_nqh] = pow(403.3,aox);
                rfield.range[p.m_nqh][0] = 0.3645;
                rfield.range[p.m_nqh][1] = 147.;
        }

        else
        {
                /* set flag saying that alpha ox will not be specified */
                lgAoxOn = false;
                /* set this to impossible number since not used, but lint needs a value */
                aox = -DBL_MAX;
                /* specify ratio, two energies */
                rfield.totpow[p.m_nqh] = p.FFmtRead();
                if( p.lgEOL() )
                        p.NoNumb("continuum ratio");

                /* assumed to be a log if negative or zero */
                if( rfield.totpow[p.m_nqh] <= 0. || p.nMatch(" LOG" ) )
                {
                        rfield.totpow[p.m_nqh] = pow(10.,rfield.totpow[p.m_nqh]);
                }

                rfield.range[p.m_nqh][0] = p.FFmtRead();
                if( p.lgEOL() )
                {
                        fprintf( ioQQQ, " There must be at least 2 numbers on this line.\n" );
                        fprintf( ioQQQ, " The ratio, and one or two energies\n" );
                        cdEXIT(EXIT_FAILURE);
                }

                rfield.range[p.m_nqh][1] = p.FFmtRead();
                /* if only one number then assume same for each */
                if( p.lgEOL() )
                        rfield.range[p.m_nqh][1] = rfield.range[p.m_nqh][0];

                if( rfield.range[p.m_nqh][0] < rfield.emm || 
                    rfield.range[p.m_nqh][1] < rfield.emm )
                {
                        fprintf( ioQQQ, " One of the energies is too low, outside the range of the code.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
                if( rfield.range[p.m_nqh][0] > rfield.egamry ||
                    rfield.range[p.m_nqh][1] > rfield.egamry )
                {
                        fprintf( ioQQQ, " One of the energies is too high, outside the range of the code.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
        }

        /* >>chng 06 mar 22, add time option to vary only some continua with time */
        if( p.nMatch( "TIME"  ) )
                rfield.lgTimeVary[p.m_nqh] = true;

        /* vary option */
        if( optimize.lgVarOn )
        {
                /* pointer to where to write */
                optimize.nvfpnt[optimize.nparm] = input.nRead;
                if( lgAoxOn )
                {
                        /* this is the number of parameters to feed onto the input line */
                        optimize.nvarxt[optimize.nparm] = 1;
                        /* vary alpha ox */
                        strcpy( optimize.chVarFmt[optimize.nparm], "RATIO alphox %f" );
                        optimize.lgOptimizeAsLinear[optimize.nparm] = true;
                        /* param is linear scale factor */
                        optimize.vparm[0][optimize.nparm] = (realnum)aox;
                        optimize.vincr[optimize.nparm] = 0.2f;
                }
                else
                {
                        /* this is the number of parameters to feed onto the input line */
                        optimize.nvarxt[optimize.nparm] = 3;
                        strcpy( optimize.chVarFmt[optimize.nparm], "RATIO LOG %f %f %f" );
                        /* param is log of abundance by number relative to hydrogen */
                        optimize.vparm[0][optimize.nparm] = (realnum)log10(rfield.totpow[p.m_nqh]);
                        optimize.vparm[1][optimize.nparm] = (realnum)rfield.range[p.m_nqh][0];
                        optimize.vparm[2][optimize.nparm] = (realnum)rfield.range[p.m_nqh][1];
                        optimize.vincr[optimize.nparm] = 0.2f;
                }
                if( rfield.lgTimeVary[p.m_nqh] )
                        strcat( optimize.chVarFmt[optimize.nparm], " TIME" );
                ++optimize.nparm;
        }

        ++p.m_nqh;
        if( p.m_nqh >= LIMSPC )
        {
                fprintf( ioQQQ, " Too many continua entered; increase LIMSPC\n" );
                cdEXIT(EXIT_FAILURE);
        }

        return;
}
#if defined(__HP_aCC)
#pragma OPTIMIZE OFF
#pragma OPTIMIZE ON
#endif