Debugging, initial conditions are still incorrectly done-- see Durran pp 343