Calculates the solar beam energy on a horizontal surface at the top of the Martian atmosphere (before atmospheric attenuation) between specified start and end times. Implements Equation 11 from Appelbaum & Flood (1990).
I_obh(Ls, phi, Ts_start, Ts_end)Beam insolation at top of atmosphere over specified time period [Wh/m²]
Areocentric longitude [deg]
Planetary latitude [deg]
Start of integration period [h]
End of integration period [h]