/* 
	 * calculate the airmass for some object.
	 * The latitude and longitude of the observer are taken from 
	 * 'observatory.h'
	 */
double tass_airmass(double ra, double dec, double jd,
                         double latitude, double longitude);