# Simple example with random data
coarse <- array(rnorm(10*20*30), dim = c(10, 20, 30)) # time x lat x lon
fine <- array(rnorm(10*40*60), dim = c(10, 40, 60)) # time x lat x lon
model <- bcsd(coarse, fine, method = "bilinear", n_quantiles = 100)
# New coarse data for prediction
new_coarse <- array(rnorm(5*20*30), dim = c(5, 20, 30)) # time x lat x lon
predictions <- predict(model, new_coarse)
# Check dimensions of predictions
dim(predictions) # Should be (5, 40, 60) for time x lat x lon
Run the code above in your browser using DataLab