#ifndef GEOCEN_C
#define GEOCEN_C
/**
* @since 2006-08-28
* @author Richard J. Mathar
*/
#include <stdlib.h>
#include <string.h>
#include <cstdarg>
#include <string>
#include <math.h>

#include <stdio.h>
#include <unistd.h>

using namespace std ;

#include "Geocen.h"

/**
* Constructor and default ctor. Defaults are the Paranal UVW platform zero
* according to VLTSW20050050.
*/
Geocen::Geocen(const double lambda, const double phi)
	: geolon(lambda), geolat(phi)
{
}

/**
* Compute the direction to the NCP in the topocentric coordinates
* Eq (21) of
* \latexonly
* \cite{MatharArxiv0608}.
* \endlatexonly
* \htmlonly
* <a href="http://arxiv.org/abs/astro-ph/0608273">astro-ph/0608273</a>.
* \endhtmlonly
* @return a vector with three components normalized to 1
*   component[0] is the N, component[1] the West and component[2] the zenith direction
*/
vector<double> Geocen::ncp() const
{
	vector<double> dir ;
	/*
	* add the components in x (North), y (West) and up z (up) direction.
	*/
	dir.push_back( cos(geolat)) ;
	dir.push_back( 0.) ;
	dir.push_back( sin(geolat)) ;

	return dir ;
}

// __oOo__
#endif // GEOCEN_C