numerical simulation of gravity wave responses in a stratified fluid
This version is a clean-up version for the previous geostrophic adjustment code
Boundary values should be updated each time inside Runge-Kutta loop. It is not needed if you choose Adams_bashforth method;