/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*ParseIonPar parse the ionization parameter command */
#include "cddefines.h"
#include "radius.h"
#include "optimize.h"
#include "rfield.h"
#include "input.h"
#include "parser.h"

void ParseIonParI(Parser &p)
{
	ParseIonPar(p,'I');
}
void ParseIonParX(Parser &p)
{
	ParseIonPar(p,'X');
}

void ParseIonPar(Parser &p,
					  char chType)
{

	DEBUG_ENTRY( "ParseIonPar()" );

	/*  check not too many continua */
	if( p.m_nqh >= LIMSPC )
	{
		/* too many continua were entered */
		fprintf( ioQQQ, " Too many continua entered; increase LIMSPC\n" );
		cdEXIT(EXIT_FAILURE);
	}

	/* this is counter for where to start scanning number on line - different
	 * for XI than for IONIZ */
	/* say that continuum is per unit area, ionization parameter*/
	strcpy( rfield.chRSpec[p.m_nqh], "SQCM" );
	if( chType == 'I' )
	{
		/* this is the usual ionization parameter, U */
		strcpy( rfield.chSpNorm[p.m_nqh], "IONI" );
	}
	else if( chType == 'X' )
	{
		/* the X-Ray ionization parameter, xi */
		strcpy( rfield.chSpNorm[p.m_nqh], "IONX" );
		/* >>chng 06 jan 20, had used 5 as in ionization parameter so
		 * possible to miss beginning of number */
	}
	else
	{
		fprintf(ioQQQ," ParseIonPar hit chCard insanity.\n");
		cdEXIT(EXIT_FAILURE);
	}

	/* get the ionization parameter*/
	rfield.totpow[p.m_nqh] = p.FFmtRead();
	if( p.lgEOL() )
		p.NoNumb("ionization parameter");

	/* check for linear option, if present take log since rfield.totpow[p.m_nqh]
	 * being log ionization parameter is the default */
	if( p.nMatch( "LINE"  ) )
		rfield.totpow[p.m_nqh] = log10(rfield.totpow[p.m_nqh]);

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

	/* vary option */
	if( optimize.lgVarOn )
	{
		if( chType == 'I' )
		{
			/* this is the usual ionization parameter, U */
			strcpy( optimize.chVarFmt[optimize.nparm], "IONIZATION PARAMETER= %f LOG" );
		}
		else if( chType == 'X' )
		{
			/* the X-Ray ionization parameter, xi */
			strcpy( optimize.chVarFmt[optimize.nparm], "XI= %f LOG" );
		}
		else
		{
			fprintf( ioQQQ, " Insanity in detecting which ionization parameter.\n" );
			cdEXIT(EXIT_FAILURE);
		}
		if( rfield.lgTimeVary[p.m_nqh] )
			strcat( optimize.chVarFmt[optimize.nparm], " TIME" );
		/* pointer to where to write */
		optimize.nvfpnt[optimize.nparm] = input.nRead;
		optimize.vparm[0][optimize.nparm] = (realnum)rfield.totpow[p.m_nqh];
		optimize.vincr[optimize.nparm] = 0.5;
		optimize.nvarxt[optimize.nparm] = 1;
		++optimize.nparm;
	}

	/* set R to large value if U specified. */
	/* 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 )
	{
		radius.Radius = pow(10.,radius.rdfalt);
	}
	/* increment nmber of specifications of continuum intensities, */
	++p.m_nqh;
	return;
}