Los objetivos de este trabajo son presentar un análisis de datos de radiación solar de calidad y desarrollar varios modelos horarios de radiación fotosintéticamente activa (PAR) utilizando combinaciones de variables radiométricas como la irradiancia horizontal global (GHI), la irradiancia horizontal difusa (DHI) y la irradiancia normal directa (DNI) a partir de sus índices adimensionales: índice de claridad atmosférica (kt), fracción difusa horizontal (kd) y fracción directa normal (kb), junto con el ángulo de elevación solar (alpha). Se utilizaron datos de GHI, DHI y DNI con frecuencias de 1 minuto en el período de 2016 a 2021 de CEDER-CIEMAT, en una meseta del norte, y PSA-CIEMAT en el sureste de la Península Ibérica, para comparar dos ubicaciones con climas muy diferentes según la clasificación de Köppen-Geiger. Se ajustaron y validaron un total de 15 modelos multilineales (con datos de entrenamiento y validación independientes) utilizando primero el conjunto de datos completo y luego por intervalos de kt. En la mayoría de los casos, los modelos que incluían el índice de claridad mostraron un mejor rendimiento, y entre ellos, los modelos que también utilizan el ángulo de elevación solar como variable obtuvieron resultados notables. Además, según la validación estadística, estos modelos presentaron buenos resultados cuando se compararon con modelos de la bibliografía. Finalmente, las estadísticas de validación del modelo indican un mejor rendimiento de los modelos por intervalos que de los modelos completos.