Datos meteorológicos entre 2012 y 2021 para estaciones meteorológicas seleccionadas de España

Datos meteorológicos entre 2012 y 2021, obtenidos mediante la librería climaemet de R. 

La segunda hoja del archivo incluye los metadatos de las variables.

Obtuve el listado completo de estaciones con aemet_stations(verbose = FALSE, return_sf = FALSE); lo transformé a csv con write.csv2(stations, "estaciones.csv"). SIN EMBARGO, hay estaciones en el listado que no tienen datos para el periodo que nos interesa (pero esto no lo sabemos hasta haber descargado los datos) así que descargamos y comprobamos con un JOIN cuáles son las que no tienen datos. Este es el listado final de las estaciones que me interesan

stations <- c('0252D','0076','0200E','0201D','0149X','0229I','0255B','0367','0370B','0372C','0324A','0016A','0002I','1387E','1387','1393','1351','1400','1437O','1473A','1428','1475X','1212E','1283U','1210X','1331A','1207U','1208A','1208H','1208','1183X','1249I','1249X','1221D','1542','1082','1078I','1055B','1057B','1059X','1083L','1159','1109','1110','1111X','1111','1014A','1024E','1050J','1014','1041A','1037Y','1549','1505','1002Y','1631E','1700X','1690B','1690A','1735X','1484','1484C','1495','2444C','2444','2811A','2117D','2331','2298','2737E','2661B','2661','2630X','2462','2400E','2374X','2235U','2401','2491C','2946X','2867','2870','2916A','2150H','2465A','2465','2030','2503X','2539','2422','2755X','2775X','2614','3391','3469','3469A','3526X','3434X','3519X','3463X','3576X','3044X','3094B','3168C','3168A','3168D','3013','3130C','3100B','3110C','3191E','3200','3129','3194U','3196','3195','3266A','3338','3111D','3175','3298X','3365A','3259','3260B','4452','4358X','4244X','4511C','4386B','4410X','4121C','4121','4220X','4103X','4147X','4148','4267X','4090Y','4560Y','4549Y','4642E','4605','4067','4061X','5973','5911A','5960','5910','5972X','5995B','5402','5427X','5390Y','5047E','5530E','5514','5051X','5582A','5860E','5298X','5181D','5038X','5270','5270B','5246','5192','5704B','5641X','5612B','5796','5783','6302A','6277B','6381','6325O','6297','6332X','6367B','6293X','6001','6268X','6106X','6058I','6084X','6155A','6156X','6172O','6032B','6205X','7096B','7247X','7002Y','7228','7119B','7012C','7145D','7209','7178I','7031','7031X','7275C','8175','8178D','8177A','8025','8019','8050X','8500A','8501','8489X','8523X','8096','8368U','8058X','8325X','8309X','8414A','8416Y','8416','8293X','9091O','9087','9091R','9111','9051','9001D','9019B','9563X','9585','9208E','9784P','9898','9201K','9814A','9170','9619','9771','9771C','9990X','9698U','9720X','9263X','9294E','9283X','9263D','9262','9981A','9987P','9573X','9381','9381I','9569A','9576C','9394X','9390','9244X','9434','9434P')

data_monthly <- aemet_monthly_period(station = stations, start = 2012, end = 2021, verbose = FALSE, return_sf = FALSE)

El xls generado tiene tres hojas: data_raw, con los datos tal cual se obtienen; data_work, habiendo extraído año (year) y mes (month) de la variable fecha; variables, con el diccionario de datos.