We may have to change the interpolation scheme at several places in roms-tools. So far, roms-tools uses linear interpolation because this is the only scheme (besides nearest neighbor interpolation) that xarray supports when interpolating over multiple dimensions.
for the topography fieldhraw, see here; since a smoother is applied after this, the interpolation scheme does arguably not matter that much
for the zonal and meridional wind forcing, as well as the SWR correction, which leads to significant differences in the MATLAB0 vs. python-generated forcing fields, see https://github.com/CWorthy-ocean/roms-tools/pull/41
For the remaining atmospheric forcing fields and tidal forcing fields, the MATLAB-based scripts use linear interpolation.
We may have to change the interpolation scheme at several places in
roms-tools
. So far,roms-tools
uses linear interpolation because this is the only scheme (besides nearest neighbor interpolation) thatxarray
supports when interpolating over multiple dimensions.The MATLAB-based scripts use a modified akima scheme for the following fields:
hraw
, see here; since a smoother is applied after this, the interpolation scheme does arguably not matter that muchFor the remaining atmospheric forcing fields and tidal forcing fields, the MATLAB-based scripts use linear interpolation.