/home66/gary/public_html/cloudy/c08_branch/source/parse_metal.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2008 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 "parse.h"

void ParseMetal(char *chCard )
{
        bool lgEOL, 
          lgGrains, 
          lgLogOn;
        long int i;
        double dmlog;

        DEBUG_ENTRY( "ParseMetal()" );

        /* parse the metals command */

        /* metal depletion factor, if negative then it is the log */
        i = 5;
        abund.ScaleMetals = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
        {
                if( nMatch("DEPL",chCard) )
                {
                        /* this option - no numbers on line but keyword depletion is
                         * deplete by set of scale factors */
                        abund.lgDepln = true;
                        for( 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 )
                        {
                                fprintf( ioQQQ, " %80.80s\n", chCard );
                        }
                        fprintf( ioQQQ, " There must be a number on this line.  Sorry.\n" );
                        cdEXIT(EXIT_FAILURE);
                }
        }

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

        if( abund.ScaleMetals <= 0. || lgLogOn )
        {
                dmlog = abund.ScaleMetals;
                abund.ScaleMetals = (realnum)pow((realnum)10.f,abund.ScaleMetals);
        }
        else
        {
                dmlog = log10(abund.ScaleMetals);
        }

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

        /* vary option */
        if( optimize.lgVarOn )
        {
                if( lgGrains )
                {
                        strcpy( optimize.chVarFmt[optimize.nparm], "METALS LOG=%f GRAINS" );
                }
                else
                {
                        strcpy( optimize.chVarFmt[optimize.nparm], "METALS LOG=%f" );
                }

                /* 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;
}

---