[spatial_combination] crashes if the mask contained only one cell on overlay mode "aggregate"

fields = c("PRENEI_Q", "PRELIQ_Q", "T_Q")
sf <- sf::read_sf(system.file("extdata/study_area/test.shp", package = "RADIS"))
clim <- get_sim2_daily(
    sf,
    DATE__greater = as.Date("2000-01-01"),
    fields = fields,
    api_format = "ncdf"
)
Error in apply((apply(xd$values, 1, expand)), 1, expand) : 
  dim(X) must have a positive length
In addition: Warning messages:
1: In create_dimension(values = x, is_raster = !sf && i %in% raster_dims,  :
  dimension value(s) non-finite
2: In create_dimension(values = x, is_raster = !sf && i %in% raster_dims,  :
  dimension value(s) non-finite
Called from: apply((apply(xd$values, 1, expand)), 1, expand)

It crashes on poly_data <- sf::st_as_sf(s_data, as_points = FALSE, merge = FALSE) on transforming raster data into vector. Maybe it's a bug in sf?