/home66/gary/public_html/cloudy/c08_branch/source/parse_absmag.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 */
/*ParseAbsMag parse the absolute magnitude command */
#include "cddefines.h"
#include "rfield.h"
#include "called.h"
#include "parse.h"

void ParseAbsMag(char *chCard, 
  long int *nqh)
{
        bool lgEOL;
        long int i;

        DEBUG_ENTRY( "ParseAbsMag()" );

        /* enter luminosity in absolute magnitudes */
        i = 5;
        strcpy( rfield.chRSpec[*nqh], "4 PI" );
        rfield.totpow[*nqh] = FFmtRead(chCard,&i,INPUT_LINE_LENGTH,&lgEOL);
        if( lgEOL )
        {
                if( called.lgTalk )
                {
                        fprintf( ioQQQ, " There should have been a number on this line.   Sorry.\n" );
                }
                cdEXIT(EXIT_FAILURE);
        }
        if( nMatch("BOLO",chCard) )
        {
                strcpy( rfield.chSpNorm[*nqh], "LUMI" );
                rfield.range[*nqh][0] = rfield.emm;
                rfield.range[*nqh][1] = rfield.egamry;
                /* page 197 allen 76 */
                rfield.totpow[*nqh] = ((4.75 - rfield.totpow[*nqh])/
                  2.5 + 33.5827);
        }
        else if( nMatch("VISU",chCard) )
        {
                strcpy( rfield.chSpNorm[*nqh], "FLUX" );
                /* this is 5550A, the center of the V filter */
                rfield.range[*nqh][0] = 0.164f;
                /* totpow(nqh) = (4.79 - totpow(nqh)) / 2.5  + 18.758
                 * page 197, allen 76, 3rd line from bottom */
                rfield.totpow[*nqh] = (-rfield.totpow[*nqh]/2.5 + 
                  20.65296);
        }
        else
        {
                if( called.lgTalk )
                {
                        fprintf( ioQQQ, " Keyword BOLOmetric or VISUal must appear.\n" );
                }
                cdEXIT(EXIT_FAILURE);
        }

        /* >>chng 06 mar 22, add time option to vary only some continua with time */
        if( nMatch( "TIME" , chCard ) )
                rfield.lgTimeVary[*nqh] = true;

        ++*nqh;
        if( *nqh >= LIMSPC )
        {
                cdEXIT(EXIT_FAILURE);
        }
        return;
}

---