Bug in Intbc when bc_method = "qm"

I am trying to use the Intbc method with quantile mapping but I'm running into an error message.

Running the example script inst/examples/example_Intbc.R at L101 I change only bc_method = "qm" and get error

Error in approx(x = fobj$par$modq[, 1], y = fobj$par$fitq[, 1], xout = x[wet], : need at least two non-NA values to interpolate

The same example works fine with bc_method = "evmos" and the following example of downscaling to points with qm also works (down_2), so it seems maybe the issue is with qm to a grid.

Let me know if you need any more info. Thanks!

Emily