USE MATH
-- Cálculo de la hora solar real
-- Entradas de la función: hora actual oficia [1-24], minuto [0-60], mes [1-12], ángulo diario [radianes], longitud[grados]
-- Salida de la función: hora solar real [decimal]
FUNCTION REAL horasolar(REAL horaactualoficial, INTEGER minuto, INTEGER mes, REAL angulodiario, REAL longitud)
DECLS
REAL Longitudmeridianocentralhuso=0
REAL diferenciaporposicion
REAL diferenciapormes
REAL Et
REAL diferenciatotal
REAL Ts
BODY
horaactualoficial=horaactualoficial+minuto/60 -- Hora actual oficial en decimal
diferenciaporposicion=(longitud-Longitudmeridianocentralhuso)*4 -- Diferencia por posición en minutos
IF (diferenciaporposicion<0.) THEN
diferenciaporposicion=-1*diferenciaporposicion
END IF
-- Diferencia por horario de verano
IF (mes<4 OR mes>10) THEN
diferenciapormes=60 -- Una hora si no estamos entre abril y octubre
ELSE
diferenciapormes=120 -- Dos horas entre abril y octubre
END IF
-- Ecuación del tiempo en minutos
Et=(0.000075+0.001868*cos(angulodiario)-0.032077*sin(angulodiario)-0.014615*cos(2*angulodiario)-0.04089*sin(2*angulodiario))*229.18
diferenciatotal=diferenciapormes+diferenciaporposicion-Et -- Diferencia total en minutos
Ts=horaactualoficial-(diferenciatotal/60)
RETURN Ts
END FUNCTION
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
