Thanks Chermelle, looking forward to a general solution. For the subset of the missing canopy height problem I have a suggested one-line update.
In the file:
~/cylc-run/u-dg767/src/contrib/Apps/CanopyHeights/canopy_heights.py
increase the neighbourhood search radius from 500 km to, say 2500km:
loop_lim_y = index_nearest_neighbour.ydist2index(trees, 2500)
This didn’t seem to slow down the ancil generation process (full app runtime was still only 30 seconds).
Results before and after, with problem cells circled in red.
Before:
After:
This also works for high-resolution domains.
Before:
After:
