Hola Belen,
En el caso del seguimiento de dos ejes, como bien decis, se soluciona imponiendo costita = 1, y por tanto rb = 1 / cosz. Esto es equivalente a setear gamma = gamma_s (azimuth solar) y beta = titaz (ángulo cenital). Evaluar el cálculo de costita con estas dos variables así, debería dar 1 en todo momento. (gamma y beta son los ángulos de posicionamiento de la superficie inclinada)
Para el caso del seguimiento azimutal, beta está fijo en un ángulo dado, y lo único que se mueve es gamma. Por tanto, hay que calcular costita con beta fijo y gamma = gamma_s (que impone el seguimiento del sol en azimuth de la superficie inclinada).
Los totales anuales de irradiación se calculan a partir de los totales diarios de irradiación si, sumando todos los días. Un promedio anual de irradiación diaria es tomar el promedio de los valores día a día. A veces se calculan también los promedios anuales de irradiancia (es común en ingeniería eléctrica, porque matchea mejor con cálculo que siempre ha hecho el sector eléctrico). En ese caso, con una serie completa, se pueden promediar todas las irradiancias incluso en series intradía, pero no hay que olvidar meter en el promedio los ceros de la noche. También se puede promediar la irradiancia de cada día (y se obtiene la irradiancia promedio diaria, incluyendo los ceros nocturnos, claro), y luego promediar estos valores diarios para obtener la irradiancia promedio anual. Ambos caminos deberían llegar al mismo número.
En concreto, el ejercicio pide trabajar con irradiación anual, así que sería total de irradiación anual o promedio de irradiación anual, como prefieran (pasar de uno a otro, debería ser un factor de 365).
Saludos, Rodrigo