Changed to sum (no trapz) for speed. Tested against Motoki's code. Works well, except for issue in forced Rossby Mode.
At this point everything compares favorably with Motoki. Not exact though, still have the Forced Rossby Wave differences.