The first main point of this project is to be able to calculate the variogram divergence between the simulation and observation variorums. This is explained technically in https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2023WR036710, Equation 66.
So the thing to do would be to load in the data for simulations and observations at a given data, then do the variogram and kriging for each, then calculate the variogram divergence score, with code that represents that equation 66.