Resultados para Observer() aparentemente no representa los efectos de elevación en PyEphem

Question

Tengo una consulta sobre los resultados dados por el módulo PyEphem relacionado con las consultas de Observer(), y los efectos de la elevación. Entiendo por un par de fuentes (como http://curious.astro.cornell.edu/question.php?number=388) que la elevación del observador tiene un efecto marcado en la hora del atardecer. Sin embargo, en el siguiente código, veo al lado de ninguna diferencia:

import ephem 

emphemObj = ephem.Observer() 
emphemObj.date = '2011/08/09' 
emphemObj.lat = '53.4167' 
emphemObj.long = '-3' 
emphemObj.elevation = 0 

ephemResult = ephem.Sun() 
ephemResult.compute(emphemObj) 
print "Sunset time @ 0m: " + str(emphemObj.previous_rising(ephemResult)) 

emphemObj.elevation = 10000 
ephemResult.compute(emphemObj) 
print "Sunset time @ 10000m: " + str(emphemObj.previous_rising(ephemResult)) 

me sale la salida:

Sunset time @ 0m: 2011/8/8 04:38:34 
Sunset time @ 10000m: 2011/8/8 04:38:34 

Estoy bastante seguro de que estoy haciendo algo mal y no siendo esto un error, pero habiendo probado varias formas diferentes, me temo que sigo dando los mismos resultados. ¿Alguien sabe lo que estoy haciendo mal aquí?

He publicado esto en https://launchpad.net/pyephem ya, pero no he tenido respuesta. Espero que no haya malinterpretado fundamentalmente el propósito de la función de elevación ...

Answer 1

El elevation de un observador significa la elevación sobre el nivel del mar de su ubicación, como la altitud de Flagstaff, Arizona, por ejemplo. Pero se presume que no solo el observador y su telescopio o binoculares están a esta distancia sobre el nivel del mar; se presume que el suelo, y por lo tanto el horizonte, son también a esta altura. Por lo tanto, un aumento de elevation no le da ninguna ventaja con respecto al horizonte, porque el horizonte se mueve con usted cuando se traslada a una ciudad de mayor altitud.

Al cabo de unos minutos con un lápiz y una libreta de papel amarillo, parece que el ángulo hacia el horizonte hza está relacionado con el radio de la Tierra r y su altura por encima del suelo h de la siguiente manera:

hza = - acos(r/(h + r)) 

Así como continuación de su ejemplo anterior:

import math 
height = 10000 
hza = - math.acos(ephem.earth_radius/(height + ephem.earth_radius)) 
emphemObj.horizon = hza 
print "Sunrise time @ 10000m: " + str(emphemObj.previous_rising(ephemResult)) 

consigo la salida:

Sunrise time @ 10000m: 2011/8/8 04:08:18 

(Tenga en cuenta que "la salida del sol" va con previous_rising() y "puesta de sol" va con next_setting()!)