Skip to content

Commit acf692c

Browse files
committed
Update the conservative regridding benchmark
1 parent 845286d commit acf692c

File tree

4 files changed

+72
-71
lines changed

4 files changed

+72
-71
lines changed

benchmarks/benchmark_utils.py

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -2,13 +2,13 @@
22
from mpl_toolkits.axes_grid1 import ImageGrid
33

44

5-
def plot_comparison(data_regrid, data_esmf, data_cdo, vmin, vmax):
5+
def plot_comparison(data_regrid, data_esmf, data_cdo, vmin, vmax, varname):
66
def relative_error(a, b):
77
return (a - b) / a
88

9-
esmf_v_regrid = relative_error(data_esmf, data_regrid).isel(time=0)["d2m"]
10-
cdo_v_regrid = relative_error(data_cdo, data_regrid).isel(time=0)["d2m"]
11-
esmf_v_cdo = relative_error(data_esmf, data_cdo).isel(time=0)["d2m"]
9+
esmf_v_regrid = relative_error(data_esmf, data_regrid).isel(time=0)[varname]
10+
cdo_v_regrid = relative_error(data_cdo, data_regrid).isel(time=0)[varname]
11+
esmf_v_cdo = relative_error(data_esmf, data_cdo).isel(time=0)[varname]
1212

1313
# Set up figure and image
1414
fig = plt.figure(1, (11, 5), dpi=250)
@@ -39,8 +39,8 @@ def relative_error(a, b):
3939
ax.set_ylabel("")
4040
axes[0].set_ylabel("latitude")
4141

42-
axes[0].set_title("xESMF vs. xarray interp")
43-
axes[1].set_title("CDO vs. xarray interp")
42+
axes[0].set_title("xESMF vs. xarray-regrid")
43+
axes[1].set_title("CDO vs. xarray-regrid")
4444
axes[2].set_title("xESMF vs. CDO")
4545

4646
plt.show()

0 commit comments

Comments
 (0)