/*coolpr stores coolants before block printed, when printing cooling agents */
#include <limits.h>
#include <string.h>
#include "cddefines.h"
#include "showme.h"
#include "coolpr.h"
#define	NCOLSAV	100

void coolpr(
	/* the line label */
	char *chLabel, 
	/* the line wavelength */
	long int lambda, 
	/* the ratio of cooling to total, negative if a heat source  */
	double ratio, 
	/* the job to do, one of "ZERO", "DOIT", or "DONE"  */
	char *chJOB )
{
	static char chLabsv[NCOLSAV][5];

	static char chSig[NCOLSAV];

	long int i, 
	  ipAr[NCOLSAV], 
	  j, 
	  limit;

	static long int nCoolant, 
	  isav[NCOLSAV];

	float SavMax, 
	  scratch[NCOLSAV];

	static float csav[NCOLSAV];

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

	/* routine is called with two flags, "ZERO" and "DONE" to 
	 * initialize and complete the printout.  Any other label is
	 * interpreted as a line label */
	if( strcmp(chJOB,"ZERO") == 0 )
	{
		/* nCoolant is the counter through the array of coolants,
		 * zero it if new job to do */
		nCoolant = 0;
	}

	else if( strcmp(chJOB,"DOIT") == 0 )
	{
		strcpy( chLabsv[nCoolant], chLabel );

		if( lambda < 10000 )
		{
			isav[nCoolant] = lambda;
		}
		else
		{
			isav[nCoolant] = lambda/10000;
		}

		csav[nCoolant] = (float)ratio;
		/* is this coolant really cooling (+) or a heat source? */
		if( ratio < 0. )
		{
			chSig[nCoolant] = 'n';
		}
		else
		{
			chSig[nCoolant] = ' ';
		}

		/* increment the counter, so this is the number actually in the stack */
		nCoolant += 1;

		/* this is limit to how much we can save */
		if( nCoolant >= NCOLSAV )
		{
			fprintf( ioQQQ, "  coolpr ran out of room, increase NCOLSAV.\n" );
			ShowMe();
			puts( "[Stop in coolpr]" );
			exit(1);
		}
	}

	else if( strcmp(chJOB,"DONE") == 0 )
	{
		/* want to print sorted list of coolants sorted from strongest to faintest */
		for( i=0; i < nCoolant; i++ )
		{
			/* save abs val so we pick up both heating and cooling */
			scratch[i] = (float)fabs(csav[i]);
		}

		for( i=0; i < nCoolant; i++ )
		{
			SavMax = 0.;
			/* following will be reset in following loop */
			ipAr[i] = -LONG_MAX;

			/* find largest of remaining coolants */
			for( j=0; j < nCoolant; j++ )
			{
				if( scratch[j] > SavMax )
				{
					SavMax = scratch[j];
					/* ipAr points to coolant within saved stack */
					ipAr[i] = j;
				}
			}
			/* set it to zero so we can look for next strongest */
			scratch[ipAr[i]] = 0.;
		}

		/* now print this stack in order or strength, seven across a line */
		for( j=0; j < nCoolant; j += 7 )
		{
			limit = MIN2(nCoolant,j+7);
			fprintf( ioQQQ, "     " );
			for( i=j; i < limit; i++ )
			{
				fprintf( ioQQQ, 
					"%4.4s%6ld%c%6.3f", 
					chLabsv[ipAr[i]], 
					isav[ipAr[i]], 
					chSig[ipAr[i]], 
					csav[ipAr[i]] );
			}
			fprintf( ioQQQ, " \n" );
		}
	}

	else 
	{
		fprintf( ioQQQ, "  coolpr called with insane job =%s=\n",chJOB );
		ShowMe();
		puts( "[Stop in coolpr]" );
		exit(1);
	}


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