Need to make the post-extraction smoothing of boundary conditions a function of the coarsest grid length-scale (i.e. max(dx_src,dx_dst) to prevent discontinuities... or the initial nearest neighbour search should be a function of this length scale. e.g. if going from fine to coarse res the number of nearest neighbours is increased to 2*(dx_dst/dx_src)
Need to make the post-extraction smoothing of boundary conditions a function of the coarsest grid length-scale (i.e. max(dx_src,dx_dst) to prevent discontinuities... or the initial nearest neighbour search should be a function of this length scale. e.g. if going from fine to coarse res the number of nearest neighbours is increased to 2*(dx_dst/dx_src)