Why not convert the limit to cos(limit) to avoid computing bth cos_zen an Zen_angles ? (zen_angles uses cos_zen...)