How to check the full(DerivativeOperators) in MethodOfLines

[DSDS-CMHL]

I am developing my shallow water equation solver, which adds the motion of the entire fluid domain as a boundary condition to the equation system. And a Boussinesq aprouch is used for solving the pde, so maybe it's difficult to write the equation just like a HeatTransfer Eq. Therefore, I only hope to obtain the spatial discretization coefficient matrix of the non-uniform grid for calculating the derivative of the velocity potential in two spatial directions x and y.

https://www.youtube.com/watch?v=okGybBmihOE
I learned from this video that DiffEqOperators can gives the DerivativeOperators(Matrix) , But DiffEqOperators.jl is being abandoned，So how can I output the Matrix using MethodOfLines.