parse_metal.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2017 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*ParseMetal parse parameters on metal command */
#include "cddefines.h"
#include "input.h"
#include "optimize.h"
#include "grainvar.h"
#include "called.h"
#include "abund.h"
#include "parser.h"

void ParseMetal(Parser &p )
{
        bool 
          lgGrains, 
          lgLogOn;
        double dmlog;

        DEBUG_ENTRY( "ParseMetal()" );

        /* parse the metals command */

        /* metal depletion factor, if negative then it is the log */
        abund.lgAbnSolar = false;       
        abund.ScaleMetals = (realnum)p.FFmtRead();
        if( p.lgEOL() )
        {
                if( p.nMatch("DEPL") )
                {
                        /* this option - no numbers on line but keyword depletion is
                         * deplete by set of scale factors */
                        abund.lgDepln = true;
                        for( long int i=0; i < LIMELM; i++ )
                        {
                                abund.depset[i] = abund.Depletion[i];
                        }
                        abund.ScaleMetals = 1.;
                        return;
                }
                else
                {
                        /* no number, so keyword, punch out */
                        if( !called.lgTalk )
                        {
                                p.PrintLine( ioQQQ );
                        }
                        fprintf( ioQQQ, " There must be a number on this line.  Sorry.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
        }

        /* sort out whether log */
        lgLogOn = false;
        if( p.nMatch(" LOG") )
        {
                lgLogOn = true;
        }
        else if( p.nMatch("LINE") )
        {
                lgLogOn = false;
        }

        if( abund.ScaleMetals <= 0. || lgLogOn )
        {
                dmlog = abund.ScaleMetals;
                abund.ScaleMetals = exp10(abund.ScaleMetals);
        }
        else
        {
                dmlog = log10(abund.ScaleMetals);
        }

        /* option to vary grain abundance as well */
        if( p.nMatch("GRAI") )
        {
                lgGrains = true;
                gv.GrainMetal = abund.ScaleMetals;
        }
        else
        {
                lgGrains = false;
                gv.GrainMetal = 1.;
        }

        /* vary option */
        if( optimize.lgVarOn )
        {
                strcpy( optimize.chVarFmt[optimize.nparm], "METALS= %f LOG" );
                if( lgGrains )
                        strcat( optimize.chVarFmt[optimize.nparm], " GRAINS" );
                /* pointer to where to write */
                optimize.nvfpnt[optimize.nparm] = input.nRead;
                optimize.vparm[0][optimize.nparm] = (realnum)dmlog;
                optimize.vincr[optimize.nparm] = 0.5;
                optimize.nvarxt[optimize.nparm] = 1;
                ++optimize.nparm;
        }
        return;
}