/*AgeCheck check various timescales after calculation complete to confirm time steady OK */
#include <string.h>
#include "cddefines.h"
#include "physconst.h"
#include "timesc.h"
#include "struc.h"
#include "zonecnt.h"
#include "printit.h"
#include "warnings.h"
#include "pressure.h"
#include "agecheck.h"

void AgeCheck(void)
{
	char chLine[133];
	long int i, 
	  limit;
	double hold,
	  tlong, 
	  tsound;

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

	/* cloud age of zero means that age command turned off
	 * negative cloud age means was not set */

	/* remember longest timescale */
	tlong = 0.;

	/* thermal equilibrium timescale */
	timesc.ttherm = 0.;
	limit = MAX2(1,ZoneCnt.nzone-1);
	for( i=0; i < limit; i++ )
	{
		timesc.ttherm = 
			MAX2( timesc.ttherm ,
			struc.pdenstr[i]*BOLTZMANN*1.5*struc.testr[i]/struc.coolstr[i]);
		/*>>chng 99 feb 01, had div by heating, changed to cooling so constant
		 * temperature models are more realistic */
	}

	tlong = MAX2(tlong,timesc.ttherm);
	if( printit.lgPrnAges )
	{
		sprintf( chLine, "   AGE: longest thermal timescale=%10.2es.", 
		  timesc.ttherm );
		notein(chLine);
	}

	if( timesc.CloudAgeSet > 0. && timesc.ttherm > timesc.CloudAgeSet )
	{
		sprintf( chLine, " C-AGE: Thermal equilibrium timescale longer than age, =%10.2es", 
		  timesc.ttherm );
		caunin(chLine);
	}

	/* check soundt travel time if constant pressure */
	if( strcmp(pressure.chCPres,"CPRE") == 0 )
	{
		tsound = timesc.sound;
		if( printit.lgPrnAges )
		{
			sprintf( chLine, "   AGE: sound travel time=%10.2es.", 
			  tsound );
			notein(chLine);
		}

		if( timesc.CloudAgeSet > 0. && tsound > timesc.CloudAgeSet )
		{
			sprintf( chLine, " C-AGE: Sound travel time longer than age in constant pressure model =%10.2es", 
			  timesc.ttherm );
			caunin(chLine);
		}
	}

	else
	{
		/* do not check if not constant pressure */
		tsound = 0.;
	}
	tlong = MAX2(tlong,tsound);

	/* molecule formation timescale */
	tlong = MAX2( tlong , timesc.BigH2MoleForm );
	tlong = MAX2( tlong , timesc.BigCOMoleForm );

	/* >>chng 97 jan 3, don't print if zero */
	if( printit.lgPrnAges && timesc.BigH2MoleForm > 0. )
	{
		sprintf( chLine, "   AGE: longest H2 formation timescale=%10.2es.", 
		  timesc.BigH2MoleForm );
		notein(chLine);
	}

	if( timesc.CloudAgeSet > 0. && timesc.BigH2MoleForm > timesc.CloudAgeSet )
	{
		sprintf( chLine, " C-AGE: H2 formation timescale longer than age, =%10.2es", 
		  timesc.BigH2MoleForm );
		caunin(chLine);
	}

	if( printit.lgPrnAges && timesc.BigCOMoleForm > 0. )
	{
		sprintf( chLine, "   AGE: longest CO formation timescale=%10.2es.", 
		  timesc.BigCOMoleForm );
		notein(chLine);
	}

	if( timesc.CloudAgeSet > 0. && timesc.BigCOMoleForm > timesc.CloudAgeSet )
	{
		sprintf( chLine, " C-AGE: CO formation timescale longer than age, =%10.2es", 
		  timesc.BigCOMoleForm );
		caunin(chLine);
	}

	/* hydrogen recombination timescale */
	timesc.threc = 0.;
	limit = MAX2(1,ZoneCnt.nzone-1);
	for( i=0; i < limit; i++ )
	{
		hold = struc.ednstr[i]*2.90e-10*pow(struc.testr[i],-0.77);
		timesc.threc = MAX2(timesc.threc , 1./hold);
	}

	tlong = MAX2(tlong,timesc.threc);
	if( printit.lgPrnAges )
	{
		sprintf( chLine, "   AGE: longest H recombination timescale=%10.2es.", 
		  timesc.threc );
		notein(chLine);
	}

	if( timesc.CloudAgeSet > 0. && timesc.threc > timesc.CloudAgeSet )
	{
		sprintf( chLine, " C-AGE: Hydrogen recombination timescale longer than age, =%10.2es", 
		  timesc.threc );
		caunin(chLine);
	}

	/* give age in various units, depending on order of magnitude */
	if( timesc.CloudAgeSet < 0. )
	{
		/* CloudAgeSet initially set to -1, if still the case then age not set */
		if( tlong < 3600. )
		{
			/* less than one day, give only seconds */
			sprintf( chLine, "  !AGE: Cloud age was not set.  Longest timescale was%10.2es.", 
			  tlong );
			bangin(chLine);
		}

		else if( tlong < 8.64e4 )
		{
			/* less than one day, give seconds and hours */
			sprintf( chLine, "  !AGE: Cloud age was not set.  Longest timescale was %8.2es = %8.2e hours.", 
			  tlong, tlong/3600. );
			bangin(chLine);
		}

		else if( tlong < 3e7/12. )
		{
			/* less than one month, give seconds and days */
			sprintf( chLine, "  !AGE: Cloud age was not set.  Longest timescale was %8.2es = %8.2e days.", 
			  tlong, tlong/86400. );
			bangin(chLine);
		}

		else if( tlong < 3e7 )
		{
			/* less than one year, give seconds and months */
			sprintf( chLine, "  !AGE: Cloud age was not set.  Longest timescale was %8.2es = %8.2e months.", 
			  tlong, tlong/3.15569e7*12. );
			bangin(chLine);
		}

		else
		{
			/* more than one year, give seconds and years */
			sprintf( chLine, "  !AGE: Cloud age was not set.  Longest timescale was %8.2es = %8.2e years.", 
			  tlong, tlong/3.15569e7 );
			bangin(chLine);
		}
	}

	else
	{
		/* age set, and passed tests, still say longest */
		if( tlong < 3e7 )
		{
			/* less than one year, give only seconds */
			sprintf( chLine, "   AGE: Cloud age was%10.2es, Longest timescale was %8.2es.", 
			  timesc.CloudAgeSet, tlong );
			notein(chLine);
		}

		else
		{
			/* more than one year, give seconds and years */
			sprintf( chLine, "  !AGE: Cloud age was %10.2es.  Longest timescale was%10.2es = %8.2e years.", 
			  timesc.CloudAgeSet, tlong, tlong/3.15569e7 );
			bangin(chLine);
		}
	}

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