- use **cupy-xarray** instead of xarray - use **cupy** instead of numpy. In particular convolving the slope could be done on GPU, see https://carpentries-incubator.github.io/lesson-gpu-programming/cupy.html - for Millan velocity and ith fields I may use **oggm** - smoothing the elevation can be probably be done using **albumentation** on gpu