/*radacl radiative acceleration due to incident continuum, calls forlin for lines*/
#include "cddefines.h"
#include "physconst.h"
#include "rfield.h"
#include "opac.h"
#include "wind.h"
#include "forlin.h"
#include "densty.h"
#include "radacl.h"

void radacl(void)
{
	long int i;
	double rforce;

#	ifdef DEBUG_FUN
	fputs( "<+>radacl()\n", debug_fp );
#	endif
	/*radiative acceleration due to incident continuum, calls forlin for lines
	 *used for wind solution derives AccelLine and AccelCont */

	rforce = 0.;
	for( i=0; i < rfield.nflux; i++ )
	{
		rforce += (rfield.flux[i] + rfield.outlin[i] + rfield.outcon[i])*
		  rfield.anu[i]*(opac.opac[i] + opac.scatop[i]);
	}

	/* radiative acceleration; xMassDensity is gm per cc, eval when PTOT called */
	wind.AccelLine = (float)(forlin()/SPEEDLIGHT/densty.xMassDensity);
	wind.AccelCont = (float)(rforce*EN1RYD/SPEEDLIGHT/densty.xMassDensity);

#	ifdef DEBUG_FUN
	fputs( " <->radacl()\n", debug_fp );
#	endif
	return;
}