/* 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 */
/*IonHelium solve ionization balance for helium */
#include "cddefines.h"
#include "dense.h"
#include "trace.h"
#include "conv.h"
#include "iso.h"
#include "ionbal.h"
#include "taulines.h"

void IonHelium( void )
{
	bool lgDebugPrint=false;

	DEBUG_ENTRY( "IonHelium()" );

	/* option to "turn off" helium */
	if( !dense.lgElmtOn[ipHELIUM] )
	{
		dense.xIonDense[ipHELIUM][0] = 0.;
		dense.xIonDense[ipHELIUM][1] = 0.;
		dense.xIonDense[ipHELIUM][2] = 0.;
		return;
	}

	/* populations */
	/* >>chng 01 may 09, add option to set ionization with element name ioniz cmnd */
	if( dense.lgSetIoniz[ipHELIUM] )
	{
		dense.xIonDense[ipHELIUM][2] = dense.SetIoniz[ipHELIUM][2]*dense.gas_phase[ipHELIUM];
		dense.xIonDense[ipHELIUM][1] = dense.SetIoniz[ipHELIUM][1]*dense.gas_phase[ipHELIUM];
		dense.xIonDense[ipHELIUM][0] = dense.SetIoniz[ipHELIUM][0]*dense.gas_phase[ipHELIUM];
	}

	lgDebugPrint = false;
#	if 0
	if( nzone > 197 )
		lgDebugPrint = true;
#	endif

	/* find ionization balance */
	ion_solver( ipHELIUM , lgDebugPrint );

	if( trace.lgHeBug )
	{
		fprintf( ioQQQ, "  %li IonHelium returns;  nzone %ld He0:%.4e He+:%.4e He+2:%.4e\n",
			conv.nTotalIoniz ,
			nzone,
			dense.xIonDense[ipHELIUM][0]/dense.gas_phase[ipHELIUM] ,
			dense.xIonDense[ipHELIUM][1]/dense.gas_phase[ipHELIUM] ,
			dense.xIonDense[ipHELIUM][2]/dense.gas_phase[ipHELIUM] );

		fprintf( ioQQQ, "     He+ /He0:%s smp:%.4e rec:%.4e ion:%.4e rad rec:%.4e 1s Pop:%.4e\n",
			iso_sp[ipHE_LIKE][ipHELIUM].chTypeAtomUsed,
			iso_sp[ipHE_LIKE][ipHELIUM].xIonSimple,
			ionbal.RateRecomTot[ipHELIUM][0],
			ionbal.RateIonizTot(ipHELIUM,0),
			ionbal.RR_rate_coef_used[ipHELIUM][0],
			iso_sp[ipHE_LIKE][ipHELIUM].st[0].Pop() );

		fprintf( ioQQQ, "     He+2/He+:%s smp:%.4e rec:%.4e ion:%.4e rad rec:%.4e\n",
			iso_sp[ipH_LIKE][ipHELIUM].chTypeAtomUsed,
			iso_sp[ipH_LIKE][ipHELIUM].xIonSimple ,
			ionbal.RateRecomTot[ipHELIUM][1],
			ionbal.RateIonizTot(ipHELIUM,1),
			ionbal.RR_rate_coef_used[ipHELIUM][1] );

		fprintf( ioQQQ, "\n" );
	}
	return;
}