Added function refract() to compute atmospheric refraction

* Added function refract() to compute atmospheric refraction:
  * simple routine, based on altitude
  * optionally correct for pressure and temperature