Hola,
Tengo una consulta sobre la ecuación 2.17, que se implementó para realizar la Figura 2.11 de las notas.
He realizado la función para calcular la cantidad de horas diurnas, usando la misma ecuación, en un script en Python, la chequeé un montón de veces, pero tengo algún error (porque no me da igual que las notas) y no puedo encontrarlo. Estoy comparando la curva para latitud -35°. Este es mi resultado:
Dejo por acá mi script;
# CANTIDAD DE HORAS DIURNAS EN HORAS Y FRACCION (TIEMPO ESTANDAR)
def horas_diurnas(n,latitud):
'''
Calcula la cantidad de horas diurnas en horas y fraccion.
Pre:
Entrada: día ordinal n, latitud
Salida: número de horas diurnas N0, expresado en horas y fracción
'''
res0 = math.tan(latitud) * math.tan(decl_solar(n))
res = (24 / math.pi) * math.acos(-res0)
return res
decl_solar(n) llama a la función que calcula la declinación solar dependiendo el día ordinal, y sé que da bien porque pude chequearlo.
Muchas gracias!
Anabela