Issue in compare_exp_obs_grids.R
I'm trying to apply a sunset workflow in which I include interpolation + bias adjustment.
Particularly, for the case of the ECMWFS5.1 seasonal predictions and ERA5-Land I obtain an error in the compare_exp_obs_grids.R function. The issue is due to a very small difference in the exp and obs grids (on the order of 2×10⁻⁸ for the conservative interpolation), which causes the following check to fail:
all.equal(as.vector(exp$coords$latitude),
as.vector(obs$coords$latitude)
I modified the code to include a tolerance for these small differences (see attached file), and it now works well. However, I’m not completely sure about the origin of these discrepancies, nor whether they could lead to any inconsistencies in the workflow.
I just wanted to report this issue in case it’s relevant for others. All the best, Vero
All the best,
Vero