Error while using meteo_ogimet() for certain stations in Switzerland about climate

[aglogowski]

Dear Radosław, The station with an error measure only wind speed and this is the specific exception that wasn't considered in the code: you can see this below: Station number 06683 https://www.ogimet.com/cgi-bin/gsynres?lang=en&ind=06683&ndays=32&ano=2020&mes=10&day=01&hora=6&ord=REV&Send=Send Station number 06631 https://www.ogimet.com/cgi-bin/gsynres?lang=en&ind=06631&ndays=32&ano=2020&mes=10&day=01&hora=6&ord=REV&Send=Send I will try to fix this issue as quickly as possible. With regards, Arkadiusz Głogowski pon., 31 maj 2021 o 11:41 Radoslaw Panczak ***@***.***> napisał(a):

…

Hi team, I'm using the package to get data for all stations in Switzerland 🇨🇭. I retrieved all stations in the country using: library(tidyverse) library(climate) nearest_stations_ogimet(country = "Switzerland", date = Sys.Date(), add_map = TRUE, point = c(8.2318, 46.7985), no_of_stations = 200 ) ch_stations_ogimet <- nearest_stations_ogimet(country = "Switzerland", date = Sys.Date(), point = c(8.2318, 46.7985), no_of_stations = 200) %>% arrange(wmo_id) I tried to loop through all of them and some seem to throw an error. Bern works fine: # example of one station bern <- meteo_ogimet( date = c(as.Date("2020-02-01"), Sys.Date() - 1), # date = c(Sys.Date() - 7, Sys.Date() - 1), interval = "daily", coords = FALSE, station = "06631") Whereas here I fail: # example of an error error <- meteo_ogimet( date = c(as.Date("2020-02-01"), Sys.Date() - 1), # date = c(Sys.Date() - 7, Sys.Date() - 1), interval = "daily", coords = FALSE, station = "06683") Error in if ((length(test[2, !is.na(test[2, ])]) == 6 & test[2, 5] == : argument is of length zero Is there any error in my code? Are the data not available for this station of formatted differently? Or is there a problem with the package? Thank you very much for help! 🙇 — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#58>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AJVWKVBACKAP7YKFCUWXLLDTQNKVNANCNFSM452ROY4Q> .

[RPanczak]

Many thanks Arek!

I wasn't aware that some stations report onl y parital data. This must be the case of the failure in retrieval.

I managed to one by one identify stations that caused the trouble and then getting all the data worked

ch_stations_search <- ch_stations_ogimet %>% 
  filter(!wmo_id %in% c("06647", "06656", "06669", "06671", "06683",
                        "06687", "06728", "06730", "06789", "06799"))

And then simply

ch_weather <- purrr::map_df(ch_stations_search$wmo_id, 
                         meteo_ogimet,
                         date = c(Sys.Date() - 7, Sys.Date() - 1), 
                         interval = "daily",
                         coords = FALSE)

It could be probably useful to be able to do it not by trial and error? Is there anyway to get some attribute of station that would allow that perhaps? Or argument to meteo_ogimet to skip those stations?

Best,
Radek

[bczernecki]

This issue has been partly addressed in climate 1.0.5. Unfortunately, CRAN policy don't allow to pass on failure so I believe the easiest solution would be to use tryCatch mechanism.
Closing.