New function airmass_ext() to compute the extinction in magnitudes due to airmass

* New function airmass_ext() to compute the extinction in magnitudes due to airmass:
  * simple function

  * airmass(): for objects below the horizon, use large but printable value, resulting
    	       in >= 30 magnitudes in extinction