RD diffusion on sphere with RM or MOL
bug fix in interp3_matrix and add unit tests