run trajectory with gfs0.25 files

[jordinagili]

Hi everyone, I am trying to run the trajectories with the SPLITR package with GFS files which have 21 hours instead of 23. Then, the hours 22 and 23 are not found in the corresponding day's file because they are on the following day. How could that be fixed? Thank you very much.

[jordinagili]

the error it gives is the following:

Error in dplyr::mutate_at():
! Can't subset columns that don't exist.
✖ Column year doesn't exist.
Run rlang::last_error() to see where the error occurred.
Warning message:
In shell(sys_cmd) :
'(cd "C:/hysplit/exec" && "C:/Users/Jordina/AppData/Local/R/win-library/4.2/splitr/win/hyts_std.exe" > NUL 2>&1)' execution failed with error code 900

rlang::last_error()
<error/vctrs_error_subscript_oob>
Error in dplyr::mutate_at():
! Can't subset columns that don't exist.
✖ Column year doesn't exist.

---

Backtrace:

1. splitr::hysplit_trajectory(...)
2. dplyr::mutate_at(...)
3. dplyr:::manip_at(...)
4. dplyr:::tbl_at_syms(...)
5. dplyr:::tbl_at_vars(tbl, vars, .include_group_vars = .include_group_vars, error_call = error_call)
   Run rlang::last_trace() to see the full context.

[jordinagili]

and the code I used:

trajectory <-
hysplit_trajectory(
lat = 42.411814,
lon = -6.979698,
height = 750,
duration = 120,
days = seq(
lubridate::ymd("2022-08-16"),
lubridate::ymd("2022-08-17"),
by = "1 day"
),
daily_hours = c(22),
direction = "backward",
met_type = "gfs0.25",
extended_met = TRUE,
met_dir = "C:/hysplit/gfs0p25",
exec_dir = "C:/hysplit/exec"
)

[jordinagili]

This just happens when I run the back trajectories from 22 and 23 hours (daily_hours = c(22,23)). For the rest of the hours it works perfectly.

[af-font]

I have the same problem, have you find a solution for that?

[gototophao]

have you solved?

[MRPHarris]

You can try using my fork of the package (https://github.com/MRPHarris/splitr), which forces the supply of an additional day of met data when running the model. Fixes this issue for me.

Once I figure out a more elegant workaround, I'll do a pull request for the main package.

detach(package:splitr, unload = TRUE) # unload splitr if needed
installr::uninstall.packages("splitr") # uninstall current splitr ver
devtools::install_github("MRPHarris/splitr") # install forked ver