/* 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 */
/*ParseCosmology parse cosmological parameters and options */
#include "cddefines.h"
#include "radius.h"
#include "rfield.h"
#include "parse.h"
#include "physconst.h"
#include "cosmology.h"
#include "parser.h"

void ParseCosmology( Parser &p )
{
	realnum param;

	DEBUG_ENTRY( "ParseCosmology()" );

	cosmology.lgDo = true;
	
	/* set an omega parameter */
	if( p.nMatch("OMEG") )
	{
		param = (realnum)p.FFmtRead();

		if( p.nMatch( "BARY") )
		{
			cosmology.omega_baryon = param;
		}
		else if( p.nMatch( "RADI") )
		{
			cosmology.omega_rad = param;
		}
		else if( p.nMatch( "MATT") )
		{
			cosmology.omega_matter = param;
		}
		else if( p.nMatch( "LAMB") )
		{
			cosmology.omega_lambda = param;
		}
		else if( p.nMatch( " K ") || p.nMatch( "CURV") )
		{
			cosmology.omega_k = param;
		}
		else
		{
			// must terminate clause safely
			fixit();// what should it do?
			TotalInsanity();
		}
	}
	else if( p.nMatch("HUBB") )
	{
		param = (realnum)p.FFmtRead();

		cosmology.h = param;

		if( param <= 0 || param > 1.1 )
		{
			fprintf( ioQQQ," This sets the variable h, which has units 100 km/s/Mpc, and is typically 0.71\n" );
		}
	}

	return;
}