Sketched out the final eigenvalue problem for the 3D internal wave modes. Need to implement one more minor function, and probably a few cases, then this should run.