/home66/gary/public_html/cloudy/c08_branch/source/parse_globule.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 */
/*ParseGlobule parse parameters off the globule command */
#include "cddefines.h"
#include "radius.h"
#include "dense.h"
#include "optimize.h"
#include "input.h"
#include "parse.h"

void ParseGlobule(char *chCard)
{
        bool lgEOL;
        long int i;

        DEBUG_ENTRY( "ParseGlobule()" );

        if( dense.gas_phase[ipHYDROGEN] > 0. )
        {
                fprintf( ioQQQ, " PROBLEM DISASTER More than one density command was entered.\n" );
                cdEXIT(EXIT_FAILURE);
        }

        /* globule with density increasing inward
         * parameters are outer density, radius of globule, and density power */
        i = 5;
        radius.glbden = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        radius.glbden = lgEOL ? 1.f : (realnum)pow((realnum)10.f,radius.glbden);
        dense.gas_phase[ipHYDROGEN] = radius.glbden;

        if( dense.gas_phase[ipHYDROGEN] <= 0. )
        {
                fprintf( ioQQQ, " PROBLEM DISASTER Hydrogen density must be > 0.\n" );
                cdEXIT(EXIT_FAILURE);
        }

        radius.glbrad = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
        {
                radius.glbrad = 3.086e18f;
        }
        else
        {
                radius.glbrad = (realnum)pow((realnum)10.f,radius.glbrad);
        }

        /* this is largest zone thickness, used to set first zone thickness */
        radius.sdrmax = radius.glbrad/25.;

        /* turn off min dr checking in NEXTDR */
        radius.lgDrMnOn = false;
        radius.glbpow = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
                radius.glbpow = 1.;
        strcpy( dense.chDenseLaw, "GLOB" );

        /* this is distance to globule */
        radius.glbdst = radius.glbrad;

        /* vary option */
        if( optimize.lgVarOn )
        {
                /* pointer to where to write */
                optimize.nvfpnt[optimize.nparm] = input.nRead;

                /* this is the number of parameters to feed onto the input line */
                optimize.nvarxt[optimize.nparm] = 3;
                strcpy( optimize.chVarFmt[optimize.nparm], "GLOBULE %f %f %f" );

                /* param is log of abundance by number relative to hydrogen */
                optimize.vparm[0][optimize.nparm] = (realnum)log10(radius.glbden);
                optimize.vparm[1][optimize.nparm] = (realnum)log10(radius.glbrad);
                optimize.vparm[2][optimize.nparm] = radius.glbpow;
                optimize.vincr[optimize.nparm] = 0.2f;
                ++optimize.nparm;
        }
        return;
}

---