A new terrain-following vertical coordinate formulation for atmospheric prediction models