/* 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 */
/*save_colden parse save column density command, or actually do the save lines output */
#include "cddefines.h"
#include "cddrive.h"
#include "input.h"
#include "save.h"
#include "parser.h"
/* this is the limit to the number of lines we can store */
#define	NPUNLM	100L

static char chElement[NPUNLM][5];
static long int nColdenEntered;
static long int ionstage[NPUNLM];

/*save_colden parse save column density command, or actually do the save lines output */
void parse_save_colden(
	Parser &p,
  /* the header for the file, a list of identifications */
  char chHeader[] )
{
	char chTemp[INPUT_LINE_LENGTH];

	long int i;

	DEBUG_ENTRY( "parse_save_colden()" );

	/* very first time this routine is called, chDo is "READ" and we read
	 * in lines from the input stream.  The line labels and wavelengths
	 * are store locally, and output in later calls to this routine */
	
	/* number of lines we will save */
	nColdenEntered = 0;
	
	/* get the next line, and check for eof */
	p.getline();
	if( p.m_lgEOF )
	{
		fprintf( ioQQQ, 
					" Hit EOF while reading line list; use END to end list.\n" );
		cdEXIT(EXIT_FAILURE);
	}
	
	while( p.strcmp( "END" ) != 0 )
	{
		if( nColdenEntered >= NPUNLM )
		{
			fprintf( ioQQQ, 
						" Too many lines have been entered; the %ld limit is.  Increase variable NPUNLM in routine save_colden.\n", 
						NPUNLM );
			cdEXIT(EXIT_FAILURE);
		}
		
		/* order on line is label (col 1-4), ionstage */
		strncpy( chElement[nColdenEntered], p.getCommand(4).c_str() , 4 );
		
		/* null terminate the string*/
		chElement[nColdenEntered][4] = 0;
		
		/* now get ionstage - 1 for atom, 2 for first ion, etc */
		ionstage[nColdenEntered] = (long)p.FFmtRead();
		if( p.lgEOL() )
			p.NoNumb("ion stage");
		
		/* this is total number stored so far */
		++nColdenEntered;
		
		/* get next line and check for eof */
		p.getline();
		if( p.m_lgEOF )
		{
			fprintf( ioQQQ, " Hit EOF while reading line list; use END to end list.\n" );
			cdEXIT(EXIT_FAILURE);
		}
		
	}
	
	/*fprintf( ioPUN, "%li lines were entered, they were;\n", 
	  nColdenEntered );*/
	/* give header line */
	
	sprintf( chHeader , "#colden %s %3li", chElement[0] , ionstage[0] );
	for( i=1; i < nColdenEntered; i++ )
	{
		sprintf( chTemp, "\t%s %3li", chElement[i] , ionstage[i] );
		strcat( chHeader, chTemp );
	}
	strcat( chHeader, "\n" );
}

/*save_colden parse save column density command, or actually do the save lines output */
void save_colden(
  /* the file we will write to */
  FILE * ioPUN )
{
	long int i;

	DEBUG_ENTRY( "save_colden()" );

	/* save some column densities  */
	double colden;
	/* save some column column densities command */
	for( i=0; i < nColdenEntered; i++ )
	{
		if( i )
			fprintf(ioPUN,"\t");
		/* get column density, returns 0 if all ok */
		if( cdColm(
				 /* four char string, null terminated, giving the element name */
				 chElement[i], 
				 /* IonStage is ionization stage */
				 ionstage[i], 
				 /* will be column density */
				 &colden) )
		{
			fprintf( ioQQQ, 
						"\n PROBLEM save_colden could not find a column density for "
						"the species with label %s %li \n\n",
						chElement[i] , ionstage[i] );
			colden = 1.;
		}
		fprintf( ioPUN, "%.4f", log10( MAX2(SMALLFLOAT , colden ) ) );
	}
	fprintf( ioPUN, "\n" );
}