/*negcon sanity check for negative continuum intensities */
#include "cddefines.h"
#include "rfield.h"
#include "zonecnt.h"
#include "phycon.h"
#include "linelabl.h"
#include "negcon.h"

void negcon(void)
{
	int lgFNeg;
	long int i;

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

	/* look for negative continuum points */
	lgFNeg = FALSE;
	for( i=0; i < rfield.nflux; i++ )
	{
		if( rfield.flux[i] < 0. )
		{
			fprintf( ioQQQ, " FLUX negative, value, freq, pointer=%10.2e%10.2e%5ld %4.4s %4.4s\n", 
			  rfield.flux[i], rfield.anu[i], i, LineLabl.chLineLabel[i]
			  , LineLabl.chContLabel[i] );
			lgFNeg = TRUE;
		}
		else if( rfield.otslin[i] < 0. )
		{
			fprintf( ioQQQ, " OTSLIN negative, value, freq, pointer=%10.2e%10.2e%5ld\n", 
			  rfield.otslin[i], rfield.anu[i], i );
			lgFNeg = TRUE;
		}
		else if( rfield.otscon[i] < 0. )
		{
			fprintf( ioQQQ, " OTSCON negative, value, freq, pointer=%10.2e%10.2e%5ld\n", 
			  rfield.otscon[i], rfield.anu[i], i );
			lgFNeg = TRUE;
		}
		else if( rfield.outcon[i] < 0. )
		{
			fprintf( ioQQQ, " OUTCON negative, value, freq, pointer=%10.2e%10.2e%5ld\n", 
			  rfield.outcon[i], rfield.anu[i], i );
			lgFNeg = TRUE;
		}
		else if( rfield.outlin[i] < 0. )
		{
			fprintf( ioQQQ, " OUTLIN negative, value, freq, pointer=%10.2e%10.2e%5ld\n", 
			  rfield.outlin[i], rfield.anu[i], i );
			lgFNeg = TRUE;
		}
	}

	if( !lgFNeg )
	{
		fprintf( ioQQQ, " No parts of the continuum were negative, the electron density was%10.2e te=%10.2e\n", 
		  phycon.eden, phycon.te );
		fprintf( ioQQQ, " This is zone number%4ld\n", ZoneCnt.nzone );
	}

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