Compute the stability of 1-Layer QG model with topography. Assume the flow is periodic in x and a channel in y. Converted from matlab to petsc4py/slepc4py