pub fn flux_solar_radiation(
    cos_incidences: &DRVector<Float>,
    albedos: &DRVector<Float>,
    sun_distance: Float
) -> DRVector<Float>