/*theavy time dependent heavy element ionization balance
 *this routine currently does not do anything */
#include "cddefines.h"
#include "timefc.h"
#include "theavy.h"

void theavy(long int NELEM )
{
	long int i, 
	  l, 
	  n, 
	  n1;
	double ar[LIMELM], 
	  change, 
	  rec;

#	ifdef DEBUG_FUN
	fputs( "<+>theavy()\n", debug_fp );
#	endif

	/* none of this can possibly work */
	n = 0;
	n1 = NELEM - 1;
	rec = 0.;
	for( i=0; i < n1; i++ )
	{
		/* rec = c( i-l+1 , 2 )*dt */
		/* following line only for lint */
		ar[i+1] = 0.;
		change = ar[i+1]*(1. - exp(-rec));
		ar[i] = (ar[i] + change)*timefcCom.timefc;
		ar[i+1] -= change;
	}

	/* following three only to trick lint */
	l = 1;
	n = l;
	n1 = n;


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