An analytical stretching function that combines the best attributes of geopotential and terrain-following vertical coordinates