Greetings, I am using the oges solver to solve for a two component grid function. For the type of scheme I am using the interpolation coefficients I am getting in the coefficient matrix are not compatible with this scheme, in the sense that they use interpolee points other than the ones I would like to use and they lead to wrong results. (my grids consist of orthogonal grids of different levels of refinement) At which part the code sets the interpolation coefficients in the coefficient matrix so that I can see how it picks the interpolee and interpolation points? If I want to customize the interpolation coefficients, you've got a hint of how should I go about in order to zero out the default interpolation coefficients at the coefficient matrix and replace them with alternative ones at the proper coefficient matrix locations? Thanks, Giorgos