/home66/gary/public_html/cloudy/c08_branch/source/parse_backgrd.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 */
/*ParseBackgrd parse options for the BACKGROUND command - this actually enters two continua*/
#include "cddefines.h"
#include "physconst.h"
#include "radius.h"
#include "rfield.h"
#include "parse.h"

void ParseBackgrd(long int *nqh, 
  char *chCard,
  realnum *ar1)
{
        bool lgEOL;
        long int i;
        double a, 
          fac, 
          rlogl, 
          z;
        char chLocal[INPUT_LINE_LENGTH];

        DEBUG_ENTRY( "ParseBackgrd()" );

        /* check that stack of shape and luminosity specifications
         * is parallel, stop if not - this happens is background comes
         * BETWEEN another set of shape and luminosity commands */
        if( rfield.nspec != *nqh )
        {
                fprintf( ioQQQ, " This command has come between a previous ordered pair of continuum shape and luminosity commands.\n Reorder the commands to complete each continuum specification before starting another.\n" );
                fprintf( ioQQQ, " Sorry.\n" );
                cdEXIT(EXIT_FAILURE);
        }

        /* diffuse x-ray background from 
         * >>refer      cosmic  background      Ostriker and Ikeuchi ApJL 268, L63.
         * parameter on cmnd is redshift, followed by optional J21 
         * >>refer      cosmic  background      Ikeuchi, S.; Ostriker, J. P., 1986, ApJ 301, 522 */

        /* >>chng 02 aug 12, do not try to enter bare powerlaw, use table power law default instead */
        strcpy( chLocal, "TABLE POWERLAW  " );
        ParseTable( nqh, chLocal , ar1);

        /* now generate equivalent of luminosity command
         * enter phi(h), the number of h-ionizing photons/cm2 */

        strcpy( rfield.chRSpec[*nqh], "SQCM" );
        strcpy( rfield.chSpNorm[*nqh], "FLUX" );
        /* this is an isotropic radiation field */
        rfield.lgBeamed[*nqh] = false;

        i = 5;
        /* this will be the redshift */
        z = FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
        {
                z = 0.;
        }

        /* optional scale factor */
        fac = FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
        {
                fac = 1.0;
        }

        /* this equation from Ostriker and Ikeuchi Ap.J.L 268, L63.
         * this should be J_nu into 4\pi sr
         * >>chng 96 jul 23, from ostriker to vedel, et al mnras 271, 743. */
        rfield.totpow[*nqh] = (log10(PI4*fac*1.e-21/
                (1.+powi(5./(1.+z),4))));

        /* this is at 1 ryd for H
         * range(nqh,1) = 1.
         *>>chgn 96 dec 18, to H ioniz pot from 1 ryd */
        rfield.range[*nqh][0] = HIONPOT;

        /* >>chng 02 aug 12, call table power law routine rather than try to enter spec here,
         * so no longer increment nspec */
        /* increment continuum indices
        ++rfield.nspec;
        if( rfield.nspec >= LIMSPC )
        {
                fprintf( ioQQQ, " Too many continua entered; increase LIMSPC\n" );
                cdEXIT(EXIT_FAILURE);
        } */
        ++*nqh;
        if( *nqh >= LIMSPC )
        {
                fprintf( ioQQQ, " Too many continua entered; increase LIMSPC\n" );
                cdEXIT(EXIT_FAILURE);
        }

        /* add fireball unless NO FIREBALL or NO CMP in */
        if( !(nMatch("O FI",chCard) || nMatch("O CM",chCard) ) )
        {

                /* put in a black body */
                strcpy( rfield.chSpType[rfield.nspec], "BLACK" );
                /* >>chng 03 may 23, CMB temp from 2.756 to 2.725 */
                rfield.slope[rfield.nspec] = (CMB_TEMP*(1. + z));
                rfield.cutoff[rfield.nspec][0] = 0.;
                rfield.cutoff[rfield.nspec][1] = 0.;
                strcpy( rfield.chSpNorm[*nqh], "LUMI" );
                a = log10(rfield.slope[rfield.nspec]);
                rlogl = log10(4.*STEFAN_BOLTZ) + 4.*a;
                strcpy( rfield.chRSpec[*nqh], "SQCM" );
                rfield.range[*nqh][0] = rfield.emm;
                rfield.range[*nqh][1] = rfield.egamry;
                rfield.totpow[*nqh] = rlogl;

                /* this flag says that CMB has been set */
                rfield.lgCMB_set = true;

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

        /* set radius to very large value if not already set */
        /* >>chng 01 jul 24, from Radius == 0 to this, as per PvH comments */
        if( !radius.lgRadiusKnown )
        {
                *ar1 = (realnum)radius.rdfalt;
                radius.Radius = pow(10.,radius.rdfalt);
        }
        return;
}

---