Abstract: A three-dimensional elastic nonhydrostatic mesoscale(β-γ) model with nested-rid is developed. It uses a set of full equations in terrain-following coordinates as its basic dynamic frame, which is solved with a time-spliting algorithm for acoustic and gravity waves. The model physical parameterization includes a K-theory subgrid eddy mixing for cloud and free atmosphere, a bulk planetary boundary layer parameterization, and three sofisticated cloud microphysics schemes with double parameters for hail cloud, warm cloud, and snow cloud, respectively. The model is designed to be used flexibly for simulation of a variety of mesoand small-scale atmospheric processes, and can be improved as a regional and local operational NWP system in future.