/home66/gary/public_html/cloudy/c08_branch/source/parse_dlaw.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 */
/*ParseDLaw parse parameters on the dlaw command */
#include "cddefines.h"
#include "dense.h"
#include "optimize.h"
#include "abund.h"
#include "input.h"
#include "radius.h"
#include "parse.h"

void ParseDLaw(char *chCard )
{
        char chCAP[5];
        bool lgEOL, 
          lgEnd;
        long int i, 
          j;

        DEBUG_ENTRY( "ParseDLaw()" );

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

        /* call fcn dense_fabden(RADIUS) which uses the ten parameters
         * N.B.; existing version of dense_fabden must be deleted
         * >>chng 96 nov 29, added table option */
        if( nMatch("TABL",chCard) )
        {
                /* when called, read in densities from input stream */
                strcpy( dense.chDenseLaw, "DLW2" );
                if( nMatch("DEPT",chCard) )
                {
                        dense.lgDLWDepth = true;
                }
                else
                {
                        dense.lgDLWDepth = false;
                }

                input_readarray(chCard,&lgEnd);
                i = 1;
                dense.frad[0] = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
                dense.fhden[0] = (realnum)FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
                if( lgEOL )
                {
                        fprintf( ioQQQ, " No pairs entered - can\'t interpolate.\n Sorry.\n" );
                        cdEXIT(EXIT_FAILURE);
                }

                dense.nvals = 2;
                lgEnd = false;

                /* read pairs of numbers until we find line starting with END */
                /* >>chng 04 jan 27, loop to LIMTABDLAW from LIMTABD, as per
                 * var definitions, caught by Will Henney */
                while( !lgEnd && dense.nvals < LIMTABDLAW )
                {
                        input_readarray(chCard,&lgEnd);
                        if( !lgEnd )
                        {
                                cap4(chCAP , chCard);
                                if( strncmp(chCAP , "END" , 3 ) == 0 )
                                        lgEnd = true;
                        }

                        if( !lgEnd )
                        {
                                i = 1;
                                dense.frad[dense.nvals-1] = (realnum)FFmtRead(chCard
                                  ,&i,INPUT_LINE_LENGTH,&lgEOL);
                                dense.fhden[dense.nvals-1] = (realnum)FFmtRead(chCard
                                  ,&i,INPUT_LINE_LENGTH,&lgEOL);
                                dense.nvals += 1;
                        }
                }
                --dense.nvals;

                for( i=1; i < dense.nvals; i++ )
                {
                        /* the radius values are assumed to be strictly increasing */
                        if( dense.frad[i] <= dense.frad[i-1] )
                        {
                                fprintf( ioQQQ, " density.in radii must be in increasing order\n" );
                                cdEXIT(EXIT_FAILURE);
                        }
                }
        }
        else
        {
                /* this is usual case, call dense_fabden to get density */
                i = 4;
                for( j=0; j < 10; j++ )
                {
                        dense.DensityLaw[j] = FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
                }

                /* set flag so we know which law to use later */
                strcpy( dense.chDenseLaw, "DLW1" );

                /* vary option */
                if( optimize.lgVarOn )
                {
                        /* NB - there are 5 = LIMEXT numbers on this line - if LIMEXT ever changes,
                         * chnage this too */
                        strcpy( optimize.chVarFmt[optimize.nparm], "DLAW %f %f %f %f %f " );

                        /* index for where to write */
                        optimize.nvfpnt[optimize.nparm] = input.nRead;
                        for( j=0; j<LIMEXT; ++j )
                        {
                                optimize.vparm[j][optimize.nparm] = (realnum)dense.DensityLaw[j];
                        }
                        optimize.vincr[optimize.nparm] = 0.5;
                        optimize.nvarxt[optimize.nparm] = LIMEXT;
                        ++optimize.nparm;
                }
        }

        /* set fake density to signal that density command was entered */
        /* real density will be set once all input commands have been read */
        /* this is necessary since density may depend on subsequent commands */
        dense.gas_phase[ipHYDROGEN] = 1.f;

        return;
}

---